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Notation and list of symbols 

Throughout this work we consider a superlattice, which is grown in the z direction. Vectors 
within the (x, y) plane parallel to the interfaces are denoted by bold face letters k, r, while 
vectors in 3 dimensional space are f,k,.... All sums and integrals extend from — cx) to cx) if not 
stated otherwise. 

The following relations are frequently used in this work and are given here for easy reference: 
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cross section 

spectral function 

electron annihilation and creation operators 

phonon annihilation and creation operators 

period of the superlattice structure 

integration and differentiation symbol 

energy 

center of energy for miniband u 

= h'^k'^/2mc kinetic energy in the direction parallel to the layers 

= 2.718 .. . base of natural logarithm 

charge of the electron (e < 0) 

electric field in the superlattice direction 

semiclassical distribution function 

Hamilton operator 

imaginary unit 

= AJ electric current. In section 5 there is an additional prefactor sgn(e) 

so that the direction is identical with the electron flow. 

current density in the superlattice direction 

Bessel function of first kind and order I 

wavevector in (x,y)-plane [i.e., plane || to superlattice interfaces] 

Boltzmann constant 

length in superlattice direction 

well indices 

effective mass of conduction band 

electron mass 9.11 x 10^^^ kg. 

number of wells 

doping density per period and area (unit [cm^^]) 

_ ^^E/ksT _ -^^-1 gQgg distribution function 

_ ^^E/kBT _|_ j^j-i pgj.j^j distribution function 
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electron density per period and area (unit [cm~^]) in well m 

Bloch vector in superlattice direction 

coupling between Wannier-states of miniband u separated by h barriers 

temperature 

bias applied to the superlattice 

= eFg^cd/h^ argument of Bessel function for irradiation 

= 2Ti/eFd argument of Bessel functions for Wannier-Stark states 

= mc/vr/i^ free-particle density of states for the 2D electron gas 

density operator 

one-particle density matrix 

indices of energy bands/levels 

chemical potential in well m, measured with respect to the bottom of the well 

frequency of the radiation field 

electrical potential 

Bloch function of band v 

Wannier function of band u localized in well m 

Wannier-Stark function of band u centered around well j 

scattering time 

Heavyside function G(x) = for x < and Q{x) = 1 for x > 

imaginary part 

real part 

principal value 

order of x" 

= ab — ba commutator 

= ab + ba anticommutator 









Fig. 1. Sketch of the spatial variation for the conduction band edge Ec{z), together with minibands i^ = a,b 
(shaded areas) for a semiconductor superlattice. 

1. Introduction 



In this review, the transport properties of semiconductor superlattices are studied. These 
nanostructures consist of two different semiconductor materials (exhibiting similar lattice con- 
stants, e.g., GaAs and AlAs), which are deposited alternately on each other to form a periodic 
structure in the growth direction. The technical development of growth techniques allows one to 
control the thicknesses of these layers with a high precision, so that the interfaces are well de- 
fined within one atomic monolayer. In this way it is possible to tailor artificial periodic structures 
which show similar features to conventional crystals. 

Crystal structures exhibit a periodic arrangement of the atoms with a lattice period a. This 
has strong implications for the energy spectrum of the electronic states: Energy bands [1] appear 
instead of discrete levels, which are characteristic for atoms and molecules. The corresponding 
extended states are called Block states and are characterized by the band index i' and the Bloch 
vector k. Their energy is given by the dispersion relation E'^{k). If an electric field F is applied, 
the Bloch states are no longer eigenstates of the Hamiltonian, but the Bloch vector k becomes 
time dependent according to the acceleration theorem 



h 



dk 
"dt 



eF. 



(1) 



where e < is the charge of the electron. Since the Bloch vectors are restricted to the Brillouin 
zone, which has a size ~ 27r/a, a special feature arises when the acceleration of the state lasts 
for a time of reioch ~ 2TTh/{eFa): If interband transitions are neglected, the initial state k is 
then reached again, and the electron performs a periodic motion both in the Brillouin zone and 
in real space [2], which is conventionally referred to as a Bloch oscillation. For typical materials 
and electric fields, TBioch is much larger than the scattering time, and thus this surprising effect 
has not been observed yet in standard crystals. 

In 1970, Esaki and Tsu suggested that superlattice structures with an artificial period d can 
be realized by the periodically repeated deposition of alternate layers from different materials 
[3]. This leads to spatial variations in the conduction and valence band of the material with 
period d implying the formation of energy bands as sketched in Fig. 1. Both the energy width 
A of these bands, as well as the extension 2TT/d of the Brillouin zone, are much smaller than 
the corresponding values for conventional conduction bands. Thus, the energy bands originating 
from the superlattice structure are called minibands. As d can be significantly larger than the 
period a of the crystal, reioch can become smaller than the scattering time for available structures 

6 



and applicable electric fields. 

It is crucial to note that the picture of Bloch-oscillations is not the only possibility to un- 
derstand the behavior of semiconductor superlattices in an electric field. The combination of a 
constant electric field and a periodic structure causes the formation of a Wannier-Stark ladder 
[4], a periodic sequence of energy levels separated by eFd in energy space. This concept is com- 
plementary to the Bloch-oscillation picture, where the frequency WBioch = eFd/h corresponds to 
the energy difference between the Wannier-Stark levels. 

Both the occurrence of Bloch oscillations and the nature of the Wannier-Stark states predict 
an increasing localization of the electrons with increasing electric field. This causes a significant 
drop of the conductivity at moderate fields, associated with the occurrence of negative differ- 
ential conductivity [3]. Similar to the Gunn diode, this effect is likely to cause the formation of 
inhomogeneous field distributions. These provide various kinds of interesting nonlinear behavior, 
but make it difficult to observe the Bloch oscillations. 

The presence of a strong alternating electric field (with frequency 0, in the THz range) along 
the superlattice structure provides further interesting features. Both photon-assisted resonances 
(shifted by hi^ from the original resonance) and negative dynamical conductance have been 
predicted on the basis of a simple analysis [5]. For specific ratios between the field strength and 
its frequency, dynamical localization [6] occurs, i.e. the dc-conductance becomes zero, which can 
be attributed to the collapse of the miniband [7] . 

Some fundamental aspects of superlattice physics have already been reviewed in Ref. [8]. 
Refs. [9,10] consider the electronic structure in detail and the review article [11] focuses on 
infrared spectroscopy. Much information regarding the growth processes as well as transport 
measurements can be found in Ref. [12]. The relation between Bloch oscillations and Wannier- 
Stark states has been analyzed in [13]. Ref. [14] provides an early review on high-frequency 
phenomena. In addition to superlattices consisting of different semiconductor materials, it is 
possible to achieve similar properties by a periodic sequence of n- and p-type doped layers [15]. 



1.1. Experimental summary 

A large variety of superlattice structures has been studied since the original proposal of Esaki 
and Tsu in 1970. These investigations can be divided into four different areas: the nonlinear 
current-field relation and its implications, the Wannier-Stark ladder, the search for Bloch oscil- 
lations, and the interaction with THz-fields. 

The simple model by Esaki and Tsu [3] predicts a nonlinear current-field relation exhibiting a 
maximum for field strengths of eFd = h/r, where r denotes the scattering time. For higher fields, 
the current drops with increasing field, yielding a region of negative differential conductivity. 
Such behavior was first observed in the experiment by Esaki and Chang [16] in 1974, where the 
conductance exhibited a sequence of dips, reaching negative values, as a function of bias voltage. 
This complicated behavior was attributed to the formation of domains with two different field 
values in the superlattice. With improving sample quality the sawtooth structure of the current- 
voltage characteristic due to domain formation could be resolved [17] more than a decade later. 
Traveling field domains (already proposed in 1977 [18]) were observed as well [19]. They cause 
self-generated current oscillations with frequencies up to 150 GHz [20]. Domain formation effects 



typically hinder the direct observation of negative differential conductance, as there is no simple 
proportionality between the measured bias and local electric field in the sample. In Ref. [21], 
the local relation between current and field could be extracted from an analysis of the global 
current-voltage characteristic. A direct observation of the Esaki-Tsu shape was possible from 
time-of-flight-measurements [22], and the analysis of the frequency response [23]. 

The concept of the Wannier-Stark ladder could be corroborated by the observation of the 
typical spacing eFd in the optical excitation spectrum of superlattice structures [24,25]. More 
recent studies refer to the transition between the Franz-Keldysh oscillations and the Wannier- 
Stark ladder [26], and the influence of higher valleys in the band structure [27]. 

The dynamical nature of Block oscillations with period TBioch was observed by transient four- 
wave mixing [28] and by a direct observation of the THz emission at WBioch [29] . Under stationary 
conditions, the phases for the oscillation cycles of individual electrons are randomized by scatter- 
ing processes, and the global signal averages out. Therefore, decaying signals have been observed 
in these experiments after a short pulse excitation, which synchronizes the dynamics in the 
very beginning. More recently, the spatial extension of the Bloch-oscillation was resolved by 
measuring its dipole field [30]. 

With the development of strong THz sources, the interaction of THz fields with transport 
through superlattices have been studied in the last few years. In particular, dynamical localiza- 
tion, photon-assisted tunneling, and absolute negative conductance were observed under irradi- 
ation by a free-electron laser [31]. Recent work aims at applying these effects to the detection of 
THz signals [32]. 

Further experiments will be discussed in the subsequent sections in direct comparison with 
the theory. 



1.2. Outline of this work 

In this work the theory of electrical transport in semiconductor superlattices is reviewed with 
a strong emphasis towards nonlinear electric transport. Here two different issues arise, which will 
be treated thoroughly: 

How can the electric transport in semiconductor structures be described quantitatively? This 
is not a straightforward issue as different energy scales like the miniband width, the scattering 
induced broadening, and the potential drop per period are typically of the same order of magni- 
tude in semiconductor superlattices. For this reason, standard concepts from bulk transport (like 
the semiclassical Boltzmann equation), which rely on the large band width, become question- 
able. Therefore, different approaches, such as miniband transport [3], Wannier-Stark hopping 
[33], or sequential tunneling [34,35], have been suggested to study the transport properties of 
semiconductor superlattices. These standard approaches imply different approximations, and 
their relation to each other was only recently resolved within a quantum transport theory [36]. 
As a result, one finds that all standard approaches are likely to fail if the miniband width, the 
scattering induced broadening, and the potential drop per period take similar values. In this 
case, one has to apply a full quantum transport calculation. For the linear response, the quan- 
tum aspects of the problem can be treated within the Kubo formula [37], which is evaluated in 
thermal equilibrium. For the nonlinear transport discussed here, this is not sufficient and a more 



involved treatment of nonequilibriuni quantum transport is necessary. An overview regarding 
different aspects of quantum transport in mesoscopic systems can be found in recent textbooks 
[38-42]. 

What is the implication of a strongly nonlinear relation between the current- density and the lo- 
cal field? As long as the field and current distribution remain (approximately) homogeneous, the 
ratio between this local relation and the global current-voltage characteristic is given by geomet- 
rical factors. In the region of negative differential conductivity, the stationary homogeneous field 
distribution becomes unstable, which may lead to complex spatio-temporal behavior. Typically, 
one observes complex scenarios, where current filaments or electric field domains form, which 
may yield both stationary or oscillating behavior. In some cases, chaotic behavior is observed as 
well. Such effects can be treated within standard concepts of nonlinear dynamics for a variety 
of different semiconductor systems [43-47] . 

The nature of quantum transport as well as pronounced nonlinearities are characteristic prob- 
lems of high-field transport in semiconductor nanostructures. In such structures the electric 
transport is determined by various quantum phenomena such as resonant tunneling (e.g. the 
resonant tunneling diode [48]), or transmission through funnel injectors (e.g. in the quantum 
cascade laser [49]). In such cases, neither standard semiclassical bulk transport models nor 
linear-response theories apply, and more advanced simulation techniques are required. The ex- 
cellent possibilities for tailoring different structures with specific superlattice periods, miniband 
widths, or doping densities make semiconductor superlattices an ideal testing ground for non- 
linear quantum transport. 

This review is organized as follows: Section 2 introduces the basis state functions, such as 
Bloch, Wannier, or Wannier-Stark states, which will be used in the subsequent sections. In 
particular, it is shown how the miniband widths and coupling parameters can be calculated 
from the material parameters on the basis of the envelope function theory. Section 3 reviews 
the three standard approaches of superlattice transport: miniband transport, Wannier-Stark 
hopping, and sequential tunneling. Each of these approaches is valid in a certain parameter 
range and allows for a quantitative determination of the current density. It is a common feature 
of these approaches that they display negative differential conductivity in qualitative agreement 
with the simple Esaki-Tsu result. These approaches can be viewed as limiting cases of a quantum 
transport theory, which is derived in section 4 on the basis of nonequilibriuni Green functions. 
The occurrence of stationary and traveling field domains in long superlattices structures is 
discussed in section 5. Here, specific criteria are presented, which allow the prediction of the 
global behavior on the basis of the current-field relation and the contact conditions. Finally, 
transport under irradiation by a THz field is addressed in section 6. Some technical matters are 
presented in the appendices. 



2. States in super lattices 

In order to perform any quantum calculation one has to define a basis set of states to be 
used. While in principle the exact result of any calculation must not depend on the choice of 
basis states, this does not hold if approximations are made, which is necessary for almost any 
realistic problem. Now different sets of basis states suggest different kinds of approximations 
and therefore a good choice of basis states is a crucial question. For practical purposes the basis 
set is usually chosen as the set of eigenstates of a soluble part Hq of the total Hamiltonian. If 
the remaining part H — Hq is small, it can be treated in lowest order of perturbation theory 
(e.g., Fermi's golden rule for transition rates) which allows for a significant simplification. This 
provides an indication for the practicability of a set of states. 

In the following discussion we restrict ourselves to the states arising from the conduction band 
of the superlattice, which is assumed to be a single band with spin degeneracy. All wave functions 
employed in the following have to be considered as envelope functions /c(^) with respect to this 
conduction band, which are determined by the Schrodinger-like equation 
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Ec{r) - V- — V + e 
Znic 



fc{r) = Efc 



where rric denotes the effective mass in the conduction band. Assuming ideal interfaces, the 
structure is translational invariant within the x and y direction perpendicular to the growth 
direction. Therefore the (x, y) dependence can be taken in the form of plane waves e where 
k and r are vectors within the two-dimensional (x, y) plane. The crucial point in this section is 
the choice of the z dependence of the basis states. This reflects the current direction considered 
here and has therefore strong implications on the description of transport. 

Semiconductor superlattices are designed as periodic structures with period d in the growth 
direction. Thus their eigenstates can be chosen as Bloch-states ^^„{z) (where q G [— vr/d, vr/ci] 
denotes the Bloch-vector and v is the band index) which extend over the whole structure. The 
corresponding eigenvalues E'^{q) of the Hamiltonian form the miniband (subsection 2.1). This 
provides the exact solution for a perfect superlattice without applied electric field. An alternative 
set of basis functions can be constructed by employing localized wave functions which resemble 
eigenstates of single quantum wells labeled by the index n. Here we use the Wannier-states 
^^{z — nd), which can be constructed separately for each miniband u (subsection 2.3). At third 
we may consider the Hamiltonian of the superlattice in the presence of a finite electric field F. 
Then the energy levels take the form E^ = Eq — neFd and one obtains the Wannier-Stark states 
^'^{z — nd), where we neglect the field-dependent coupling between the subbands (subsection 
2.4). The spatial extension of these states is inversely proportional to the electric field. In the 
subsequent subsections the different basis sets will be derived and their properties will be studied 
in detail. 

2.1. Minibands 

The periodicity of the superlattice structure within the z direction implies that the eigenstates 
of the Hamiltonian can be written as Bloch states (Pq{z), where q G [—TT/d,Tr/d] denotes the 
Bloch vector. The construction of these eigenstates can be performed straightforward within 
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the transfer matrix formulation, see e.g. [50,51]. Within a region Zj < z < Zj+i of constant 
potential and constant material composition the envelope function can be written as f{z) = 
^^.eifej(^)(^-^j)+5^.e-i'=^(^)(^-^j). Then the connection rules [52] (see also Refs. [10,51,53] for a 
detailed discussion) 



/c(r) 
'/c(f) 



m. 



■cj 



dz 



^j+i 



Zj + l 



+0- 



7c(r)._ 
1 



dfcif) 



+0- m-cj+i dz 



Zj + l 



+0+ 



(3) 

(4) 



apply|3]where rric^j is the effective mass in region Zj < z < Zj^i. 
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Mj{E) 



\ 



1 + -jj- — 

kjUlcJ + l 
kj + lTTlcJ 



^ikj{zj+i-Zj) 



1 

1 + 



kj+iiricj 



-ikj{z 



'j+ 



i-zA \ 



-ikj{zj+i-Zj) 



(6) 



/ 



If a single period of the superlattice consists of M regions with constant material composition, 
the Bloch-condition (pq{z + d) = e^'^'^ ipq{z) implies 




M 



A, 







(7) 



For standard superlattices (M = 2) the solutions resemble those of the Kronig-Penney model 
[61], except for the use of effective masses associated with the connection rule (4). Within the 
transfer formalism the extension to superlattices with a basis [62] (larger M) is straightforward. 
For given q, Eq. (7) is only solvable for selected values of E which define the miniband-structure 
E^{q). An example is shown in Fig. 2. Next to the envelope function approximation discussed 
here, different approaches can be used to calculate the superlattice band structure [9,63]. 
For a given miniband v the following quantities may be defined: 

/■n /d 
dqE^iq) (8) 

-n/d 

miniband width: A'' = MaXg{S'"(g)} - Minj£;''(g)} (9) 



^ Throughout this work we use the energy dispersion E{k) = Ec + \^ with the band edge Ec = 0.8a: meV 
and the effective mass m^ = (0.067 + 0.083a;)mo for the conduction band of Al^jGai-iiAs with x < 0.45 [54]. For 
GaAs/AlAs structures, nonparabohcity effects are included using the energy-dependent effective mass mc(E) = 
mc{E — Ev) / {Ec — Ev) with the parameters rric = 0.067mo, Ec = 0, Ev — —1.52 meV for GaAs and iric — 0.152mo, 
Ec — 1.06 meV, Ev = —2.07 meV for AlAs [55,56], where £„ denotes the edge of the valence band. X,L-related 
effects are neglected for simplicity. They become relevant in some transport studies [57]. Some approaches to the 
theoretical study of tunneling via these minima can be found in [58-60]. 



11 



U.8- 




0.6- 


^^^\ /-^ 


? 0-4- 


'^-^^^^^^^-^-"""^ 






§0.2- 


—^^ ^^^ 


0.0- 


— — a 

1 1 1 1 1 1 1 1 



-1.0 



-0.5 0.0 

q [7t/d] 



0.5 



1.0 



Fig. 2. Calculated miniband structure for the GaAs-Alo.3Gao.7As superlattice used in [64] with well width 6.5 nm 
and barrier width 2.5 nm. 

which characterize the miniband structure. For the miniband structure shown in Fig. 2 the values 
E"- = 54.5 meV, A" = 23.6 meV, E^ = 220 meV, A'' = 98 meV, E" = 491 meV, and A'^ = 233 
meV are found, where the band indices are labeled hy i' = a, b,c. .. . The increase of A with i' 
can be easily understood in terms of the increasing transparency of the barrier with the electron 
energy. As E'^{q) is restricted to —ir/d < q < ir/d and E'^{q) = E'^{—q) by Kramers degeneracy, 
the function E^{q) can be expanded as follows: 



E'^iq) =E'' + Y^ 2Tj^ cos{hdq) 



h=l 



(10) 



Typically, the terms T^ for h > 2 are much smaller than the T^ term. E.g., for the bandstructure 

0.48 meV for the lowest band, T^ = 23.7 meV, 



of Fig. 2 one obtains Tf 



-5.84 meV 
2.1 meV for the second band, and Tf 



rpa 
^2 



-53.3 meV, 



^2 



1 

5.3 meV for the third band. 

4|r{'|. 



This demonstrates that the band structure is essentially of cosine-shape and thus A'' 

The dispersion E'^{q) = E^ + 2T^ cos{dq) can be viewed as the result of a standard tight-binding 

calculation with next-neighbor coupling T^ . 

In order to perform many-particle calculations the formalism of second quantization, see, e.g., 
[65], is appropriate. Let aq and a'^ be the creation and annihilation operator for electrons in the 
Bloch-state of band u with Bloch-vector q. Then the Hamiltonian reads 






[ID 



which is diagonal in the Bloch-states, as the Bloch-states are eigenstates of the unperturbed 
superlattice. 

2.2. Bloch-states of the three-dimensional superlattice 

In the preceding subsection only the z direction of the superlattice was taken into account. For 
an ideal superlattice Eq. (2) does not exhibit an (x, y)-dependence and thus a complete set of 
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i?^(<Z,k)«i?'^(g,0) + (v.,^,ol:T--73Ko> (13) 



eigenstates states can be constructed by products of plane waves e'*''" /(27r) and a z-dependent 
function fk{z) which satisfies the eigenvalue equation 

As Ec{z) and mc{z) are periodic functions with the superlattice period d, the eigenstates are 
Bloch state fkiz) = ^^g\JyZ) with energy E'^{q, k).[^ Within first order perturbation theory in k^ 
one obtains the energy 

' 2mc(2;) ' 

Now (/9^q(z) will exhibit a larger probability in the well, so that it seems reasonable to replace 
the second term by 

Ek = ,^ 14 

where m^ is the effective mass of the quantum well. In analogy to Eq. (11) the full Hamiltonian 
reads 

^SL = Y.r d^ / d'^ [^'(9) + ^'^] <(k)<(k) (15) 

^ J—ir/d J 

where flq (k) and aq(k) are the creation and annihilation operator for electrons in the Bloch- 
state of band v with Bloch-vector q and wave vector k in [x, y) plane. In order to evaluate matrix 
elements for scattering processes the zeroth order envelope wave-functions (p'^ q{z) e / {2tt) are 
applied in subsequent sections. The treatment is completely analogous for Wannier and Wannier- 
Stark states discussed in the subsequent subsections. 

2.3. Wannier functions 

By definition the Bloch-functions are delocalized over the whole superlattice structure. This 
may provide difficulties if electric fields are applied or effects due to the finite length of the 
superlattice are considered. Therefore it is often helpful to use different sets of basis states which 
are better localized. A tempting choice would be the use of eigenstates of single quantum wells, 
see, e.g., [66,67]. Nevertheless such a choice has a severe shortcoming: The corresponding states 
are solutions of two different Hamiltonians, each neglecting the presence of the other well. Thus 
these states are not orthogonal which provides complications. Typically, the coupling is estimated 
by the transfer Hamiltonian [68] within this approach. 

For these reasons it is more convenient to use the set of Wannier functions [69] 

^"{z - nd) = \ -^ r dq e-'"^'^' <f''(z) (16) 



^ Eq. (12) shows that the effective Hamiltonian is not exactly separable in a z and r-dependent part, as the 
z-dependent effective mass affects the k-dependence, describing the behavior in the {x, y) plane. Nevertheless, this 
subtlety is not taken into account here, as discussed below. 
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Fig. 3. Wannier functions for the superlattice from Fig. 2. Full line: 'i'°{z), dashed line '^''{z), The thin line 
indicates the conduction band edge profile. 

which are constructed from the Bloch functions with the normahzation J dz[ip'^,{z)]*(pq{z) = 
^vu'^{q — q')- Here some care has to be taken: The Bloch functions are only defined up to a 
complex phase which can be chosen arbitrarily for each value q. The functions ^(z) depend 
strongly on the choice of these phases. In [70] it has been shown that the Wannier functions 
are maximally localized if the phase is chosen in the following way at a symmetry point Zsym 
of the superlattice: If V'o(zsym) / choose ipg{zsyra) G ^ (i)- Otherwise choose iipq{zsyra) G ^ 
(ii). In both cases the phase is chosen such, that (Pq{z) is an analytic function in q. (The latter 
requirement defines the phase when (Pq{zsym) becomes zero and prevents from arbitrary sign 



changes of (pq{z< 



sym)-) For such a choice the Wannier functions ^^(z) are real and symmetric 
Now there are two symmetry points, one in the center of 



(i) or antisymmetric (ii) around Zsym 
the well (z^^) and one in the barrier 



U, 



barrier '^ 
sym 



for a typical superlattice. If the energy of the 
miniband is below the conduction band of the barrier, the Wannier functions seem to be strongly 
localized for z^^, while z^^^^^'^ may be suited as well for larger energies, where the minibands 
are above the barrier. This point has also been addressed in [71]. In Fig. 3 the Wannier functions 
for the first two bands (using z^^) are displayed. One finds that both functions are essentially 
localized to the central quantum well where they resemble the bound states. Outside the well 
they exhibit a decaying oscillatory behavior which ensures the orthonormality relationP] 



dz ^"{z - nd)^'^{z - md) 



^n,m"u,ij. 



(17) 



Within second quantization the creation a' 
ciated with the Wannier functions ^'^( 



^\ 



and annihilation a^ operators of the states asso- 



nd) are defined via 



^ Note that the orthonormality is not strictly fulfilled for different bands u, fi if an energy dependent effective 
mass is used. In this case, energy-dependent Hamiltonians (2) are used for the envelope functions, and therefore 
the orthonormality of eigenfunctions belonging to different energies is not guaranteed. In principle this problem 
could be cured by reconstructing the full wave functions from the envelope functions under consideration of the 
admixtures from different bands. 
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2tt 



•^-'"'^'^ al 



(18) 



Inserting into Eq. (11) and using Eq. (10) one obtains the Hamiltonian within the Wannier basis 



^SL = E 



E^a'^^^al + Y, n (<i,< + <Ih< 



h=l 



(19) 



As the Wannier functions are hnear combinations of Bloch functions with different energies, 
they do not represent stationary states. Neglecting terms with h > 1 the time evolution of the 
annihilation operators in the Heisenberg representation is given by 



For the initial condition a5^(t = 0) = SnflO'Q this set of equations has the solution 



<(t) = i-V„ 



2T^ 



t e" 






(20) 



(21) 



where Jn{x) is the Bessel function of first kind [72]. This shows that the initially occupied 
Wannier state decays on a time scale of 

h 



TWa 



2rf 



(22) 



.^^1 



At this time (oq (t)aQ(t)) = Jo(l) ^ 0.586, thus rwannicr may be viewed as a kind of half-life 
period, although there is no exponential decay. 

If an electric field F is applied to the superlattice, the additional potential (/>(z) = —Fz has 
to be taken into account. Within the Wannier basis the corresponding terms of the Hamiltonian 
can be evaluated directly by the corresponding matrix elements J dz ^'^(z — md)eF z'^^ {z — nd). 
Including the parallel degrees of freedom k, the total Hamiltonian H = Hq + Hi + H2 reads: 



Ho = Y f d^KE" + Ek- eFRl" - eFdn)a^^(k)a^(k) 



(23) 



^^=Y.Y. ^^K^h [a:;Vjk)<(k) + a^t(k)a5;^,(k) 

n,u,^ih=l 



K,u 



eFK 



-^tni^KX^) + <\^)<+h{^)]} 



^2 = E / d^fe(-eFi?r)<^(k)a:;(k) 

n 1/ (I "^ 



(24) 
(25) 






with the couplings i?^^ = J dz ^^{z — hd)z^'^{z). If the superlattice exhibits inversion symmetry 
the coefficients R'j^'^ vanish for h > 0. Finally note, that the expression of H is still exact, except 
for the separation of z and (x, y) direction. 
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The term Hq describes the energy of the states in the superlattice neglecting any couphngs 
to different bands or different wehs. Hi gives the couphng between different wells. Finally H2 
describes the field-dependent mixing of the levels inside a given well. In particular it is responsible 
for the Stark shift. 

The term Hq + H2 for n = can be diagonalized [73] by constructing the new basis 

|/i) = j;C/^^|z.) (26) 

satisfying {Hq + il'2)\li) = E^\jl) where the columns of V^^ are the components of |/i) with 
respect to the basis \iy). This shows, that the level separation becomes field dependent, which 
has been recently observed in superlattice transport [74] under irradiation. In the new basis the 
Hamiltonian is given hy H = Hq^^ + Hf^^^ with 

^ren ^Y^ f ^2j^ ^^u + Ek- eFdn)a^1"(k)a^(k) (27) 

n,i/ ■^ 

00 „ 

Hr =J2J2 r^'^ ^r [<U^)<i^) + a:;t(k)G^^,(k)] (28) 

and the matrix elements 

^f = E i^""^) * (^h^-'>^ - ^FRh"W . (29) 



fJ,,U 



It will turn out later that in the limit of sequential tunneling it is more appropriate to use Hl*^^ 
as a perturbation instead of Hi. 

2.4. Wannier-Stark ladder 

If an electric field F is applied to the superlattice structure the Hamiltonian exhibits an 
additional scalar potential e(f>{z) = —eFz which destroys the translational invariance. In this 
case we can easily see: If there exists an eigenstate with wavefunction <I>o('2) and energy Eq, 
then the set of states corresponding to wavefunctions ^j{z) = ^o{z — jd) are eigenstates of the 
Hamiltonian with energies Ej = Eo—jeFd as well. These states are equally spaced both in energy 
and real space and form the so-called Wannier-Stark ladder [4] . This feature has to be considered 
with some care, as the potential e(j){z) is not bounded for the infinite crystal, which implies 
a continuous energy spectrum [75]. Nevertheless, the characteristic energy spectrum of these 
Wannier-Stark ladders could be resolved experimentally [24,25] in semiconductor superlattices. 
For a more detailed discussion of this subject see [13,76,77]. 

If one restricts the Hamiltonian H^h in Eq. (11) to a given miniband v, an analytical solution 
for the eigenstates of HsL — eFz exists [78] : 
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Fig. 4. Wannier-Stark states calculated from Eq. (30) for the superlattice from Fig. 2. The thin line indicates the 
conduction band edge profile, (a) Full line: $"(2), dashed line $''(z), (b) $"(2) for different fields. 

where E'^ = E^ — jeFd is the ladder of energies corresponding to the u miniband with av- 
erage energy E'^ similar to the discussion of the Wannier statesp]. The field-induced coupling 
to different bands induces a finite lifetime of these single band Wannier-Stark states due to 
Zener tunneling [78]. Thus these Wannier-Stark states can be viewed as resonant states (an 
explicit calculation of these resonances has been performed in [79]). For a cosine-shaped band 
E^[q) = E^ + 2T^ cos{qd) the Wannier-Stark states from Eq. (30) can be expanded in Wannier 
states 1^;^,) [80] 






E 



(7rj 



"^ \eFd 



^,^) for next neighbor coupling 



(31) 



where the definition (16) has been used. This relation can be obtained directly by diagonalizing 
Eqs. (23,24) within the restriction to a single band, nearest-neighbor coupling Ti, and Rq'^ = 0, 



I.e. z. 



sym 



0. In Fig. 4 examples for the Wannier-Stark states are shown. It can be clearly seen 
that the localization of these states increases with the electric field. They exhibit an oscillatory 
structure within a region of approximately ISy jeFd periods and a strong decay outside this 

= x'^12 [Eq. (8.536(2) 



2 72/ 



region. This magnitude can be estimated via Eq. (31). As ^Y^=-oo '^ J-, 

of [81]] we can conclude that Jn I -^^ ) deviates from zero essentially in the range —2\T^\leFd 



\eFdj 

n < 2\T^\/eFd which, together with /V 



< 



4|Tf I, provides the result given above. 



* This representation depends crucially on the relative choice of phases in the Bloch functions ifiqiz). The situation 
resembles the construction of Wannier states (16) and it is suggestive to use the same choice of phase although 
I am not aware of a proof. For consistency the origin of 2 has to be chosen such that Zsym — holds for the 
symmetry point of the superlattice. 
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3. The standard approaches for superlattice transport 

If an external bias is applied to a conductor, such as a metal or a semiconductor, typically an 
electrical current is generated. The magnitude of this current is determined by the band structure 
of the material, scattering processes, the applied field strength, as well as the equilibrium carrier 
distribution of the conductor. In this section the question is addressed, how the special design of 
a semiconductor superlattice, which allows to vary the band-structure in a wide range, influences 
the transport behavior. Throughout this section we assume that a homogeneous electrical field 
F is applied in the direction of the superlattice (the z direction) and consider the current parallel 
to this field. Due to symmetry reasons the transverse current parallel to the layers should vanish. 

A very elementary solution to the problem has been provided by Esaki and Tsu in their pio- 
neering paper [3]. Consider the lowest miniband of the superlattice labeled by the superscript a. 
The eigenstates are the Bloch-states ^q{z) with the Bloch- vector q and the dispersion is approx- 
imately given by E^'lq) ~ E"" — 2\T^\ cos{qd) (see Sec. 2.1 for details) as depicted in Fig. 5a. At 
low temperatures the states close to the minimum at (7 ~ are occupied in thermal equilibrium. 
If an electric field is applied (in z-direction) the Bloch-states are no longer eigenstates of the full 
Hamiltonian but change in time. According to the acceleration theorem [1] 

the states remain Bloch states in time, but the Bloch-vector becomes time dependent and we 
find q{t) = eFt/h if the electron starts in the minimum of the band at t = 0. For t = ■Kh/{eFd) 
the boundary of the Brillouin zone {q = vr/d) is reached. This point is equivalent with the 
point at q = — vr/d, so that the trajectory continues there which is often called Bragg-reflection. 
Finally, at t = TBioch = 2TTh/{eFd) the origin is reached again. Neglecting transitions to diff'erent 
bands (Zener transitions, whose probability is extremely small for low fields) the state remains 
in the given band and thus the same state is reached after TBioch resulting in a periodical motion 
of the state through the Brillouin zone [2]. This oscillation is called Bloch-oscillation and is 
quite general for arbitrary bandstructures. It could be observed in superlattices [28,29]. The 
Bloch-states q travel with the velocity 

1 dE-{q) ^ 2d\Tf\ . 
.(g) = -^^«^— sm(,d) (33) 

Thus, we find v{t) = VmSm{eFdt/h) and the position of a wave packet z(t) = zq + Zm{^ — 
cos{eFdt/h)} with Vm = 2d\Ti\/h and Zm = 2|Tf |/eF. In [82] this behavior has been nicely 
demonstrated by an explicit solution of the Schrodinger equation. The spatial amplitude of this 
oscillation has been resolved recently [30]. 

Scattering processes will interrupt this oscillatory behavior. As scattering processes are likely 
to restore thermal equilibrium it makes sense to assume that the scattered electron will be 
found close to q = 0, the initial point used before. As long as the average scattering time r is 
much smaller than Trh/{2eFd) the electrons will remain in the range < g < t: /{2d) where the 
velocity increases with q and thus an increase of F will generate larger average drift velocities. 
Thus for eFd <^ 7r/i/(2r) a linear increase of Wdrift(-^) is expected. In contrast, if r > ■7rh/{eFd) 
the electrons reach the region —n/d < q < with negative velocities and thus the average drift 
velocity can be expected to drop with the field for eFd > Trh/r. For high fields, eFd ^ 27r/i/r, the 
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Fig. 5. (a) Dispersion E°'{q) and velocity w"(g) for the lowest miniband. (b) Velocity-field relation according to 
Esaki-Tsu [3]. 

electrons perform many periods of the Bloch-oscillation before they are scattered and thus the 
average drift velocity tends to zero for F ^ oo. K detailed analysis for a constant (momentum- 
independent) scattering time r gives the Esaki-Tsu relation [3]: 

eFdh/T 



VAriit{F) = Ve.t{F) 



(34) 



""{eFdf + ih/Tf 

This result will be derived in Sec. 3.1.1 as well. The drift velocity exhibits a linear increase with 
F for low fields, a maximum at eFd = h/r and negative differential conductivity for eFd > ?i/r, 
see Fig. 5b. This general behavior could be observed experimentally [83,22]. 

The rather simple argument given above neglects the plane wave states in (x, y)-direction, 
the thermal distribution of carriers, and treats scattering processes in an extremely simplified 
manner. In Section 3.1 a more realistic treatment is given within the miniband transport model 
where the electrons occupy Bloch-states and the dynamical evolution of the single states is 
described by the acceleration theorem (32). 

A complementary approach to miniband transport is the use of Wannier-Stark states, which 
are the 'real' eigenstates of the superlattice in an electric field [see Sec. 2.4 for a discussion of 
the problems involved with these states]. Scattering processes cause transitions between these 
states yielding a net current in the direction of the electric field [33]. This approach is called 
Wannier-Stark hopping and will be described in detail in Sec. 3.2. 

For superlattices with thick barriers (i.e. narrow minibands) it seems more appropriate to 
view the structure as a series of weakly-coupled quantum wells with localized eigenstates. Due 
to the residual coupling between the wells tunneling processes through the barriers are possible 
and the electrical transport results from sequential tunneling from well to well, which will be 
discussed in Sec. 3.3. Generally, the lowest states in adjacent wells are energetically aligned for 
zero potential drop. Thus, the energy conservation in Fermi's golden rule would forbid tunneling 
transitions for finite fields. Here it is essential to include the scattering induced broadening of 
the states which allows for such transitions. 
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Fig. 6. Overview of the difFerent standard approaches for superlattice transport. 

These three complementary approaches are schematically depicted in Fig. 6. They treat the 
basic ingredients to transport, band structure (coupling Ti), field strength (eFd), and scattering 
(F = h/r), in completely different ways. In section 4 and appendix B these approaches will 
be compared with a quantum transport model, which will determine their respective range of 
validity as sketched in Fig. 7. 

It is an intriguing feature that all three approaches provide a velocity-field relation in quali- 
tative agreement with Fig. 5b, except that the linear increase is missing in the Wannier-Stark 
hopping model. Therefore the qualitative features from the Esaki-Tsu model persist but details 
as well as the magnitude of the current may be strongly altered. These points will be discussed 
in detail in the subsequent subsections. 

3.1. Miniband transport 

Conventionally, the electrical transport in semiconductors or metals is described within a 
semiclassical approach. Due to the periodicity of the crystal, a basis of eigenstates from Bloch- 
functions can be constructed. For the superlattice structure considered here it is convenient 
to treat the Bloch-vector q in superlattice direction and the Bloch-vector k in the direction 
parallel to the layers separately. The eigenfunctions are c/?J^(z) e'*''" /27r and the corresponding 
energies are -E(g, k, i/) = E^{q) + E^ with the superlattice dispersion E^{q) and the in-plane 
energy Ej. = /i^k^/2m,c, for details see Sec. 2.1. The occupation of these states is given by the 
distribution function /(g, k, i/, t) describing the probability that the state (g, k, i/) is occupied 
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Fig. 7. Ranges of validity for the different standard approaches for superlattice transport. Miniband transport 
holds for eFd <^ 2|Ti| and F <g 2|ri|; Wannier-Stark hopping holds for F <^ eFd; Sequential tunneling holds for 
2!ri| < eFd or 2|ri| < F, from [36]. 

and /(g, k, I/, t) d(7d^A;/(27r)^ is the particle density within the the volume element dqd'^k in 
momentum space. 

An electric field breaks the translational invariance of the system and the Bloch states are no 
longer eigenstates. Within the semiclassical theory the temporal evolution of the distribution 
function is given by the Boltzmann equation[_] 

df 
dt 



df{q,k,i^,t) ^ eFdf 
dt h dq 



(35) 



which is derived in most textbooks on solid state physics. (The inclusion of magnetic fields 
is straightforward and results for superlattice structures are given in [84,85]; see also [86] for 
corresponding experimental results. Recently, a mechanism for spontaneous current generation 
due to the presence of hot electron in a magnetic field perpendicular to the superlattice direction 
was proposed [87].) The right side describes the change of the distribution function due to 
scattering. For impurity or phonon scattering the scattering term reads 

/ dfiq,k,u,ty 
V dt 

Y, f dq' J d\' {P{q', k', z.' ^ q, k, u)fiq', k' , v\ t)[l - f{q, k, u, t)] (35) 



scatt 



P{q, k,i^^q',k', u')fiq, k, v, t) [1 - f{q', k',i^', t)] } . 



^ If spatially inhomogeneous distributions /(r,g, k, ;/, t) are considered, the convection term v'^il)^ + h]<./mc- ^ 
has to be added on the left hand side. 
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P{q,\i.,i' -^ g',k',z^') denotes the scattering probability from state (g, k, zv) to state {q',k.',u') 
which can be calculated by Fermi's golden rule using appropriate matrix elements for the different 
scattering processes. Details regarding these scattering processes (for bulk systems) can be found 
textbooks, such as [88,89]. Specific calculations for superlattice structures can be found in [90- 
92]. Once the distribution function is known, the current density Jmbt for miniband transport 



in z direction can be evaluated directly by 

Jmbt = ^^^^^ E f^] ^^ / ^'^ ^ K -, tX{q) (37) 



fTv/d 

-n/d 

and the electron density per superlattice period (in units [1/cm^]) is given by 



2(for Spin)d ^-^ /"^/'^ , f . 



(27r)S 



VT dq f d^kf{q,k,u,t). (38) 

^J-TT/d J 



This approach for the electric transport is called miniband transport. An earlier review has been 
given in [93] where several experimental details are provided. 

One has to be aware that Boltzmann's equation holds for classical particles under the assump- 
tion of independent scattering events. The only quantum mechanical ingredient is the use of the 
dispersion E^{q, k), thus the term semiclassical approach is often used. Therefore deviations may 
result from various quantum features, such as scattering induced broadening of the states, the 
intracollisional field effect, or correlations between scattering effects leading, e.g., to weak local- 
ization. While these features are notoriously difficult to describe, operational solution methods 
like Monte-Carlo methods [94,95] exist for the Boltzmann equation explaining the popularity of 
the semiclassical approach. 

3.1.1. Relaxation time approximation 

Boltzmann's equation can be solved easily if the scattering term is approximated by 

( df{q,k,t) \ ^ nF{E{q,k)-fi)-f{q,k,t) ^^^^ 

V 9t y^^^.j rscatt(g,k) 

with the relaxation time Tsca.ttiq,k). The relaxation time approximation is correct in the linear 
response regime for a variety of scattering processes. Here it is applied to nonlinear transport, 
in order to obtain some insight into the general features. The underlying assumption is that any 
scattering process restores the thermal equilibrium described by the Fermi function npiE) = 
[eyip{E/kBT) + 1]~^ and the chemical potential fi. (The discussion is restricted to the lowest 
miniband here and thus the miniband index i' is neglected.) Then the stationary Boltzmann 
equation reads: 

eFTsci,tt{q,k)df{q,k) , „. , ^ rjpf i\ \ fAn\ 
^ g + f{q,^)=nF{E{q,k)- n). (40) 

This is an inhomogeneous linear partial differential equation which, together with the boundary 
condition /(— Tr/d, k, t) = /(vr/d, k, t), can be integrated directly and one finds: 
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/(g,k) 



■n/d 
-■Kid 



e-FTscatt(go,k) 
&{q - go) + 



■ exp 



dq' 



h 



90 



ei^rscatt(g',k) 



X < 



exp 






T/d 



171^ 



1 



(41) 



Assuming a constant scattering rate Tscatt(<?;k) = r, the g' integrals become trivial. For the 
simplified miniband structure E{q) = E"" — 2|T"| cos{qd) one obtains the electron density from 
Eq. (38) 

2d f^^^ f 

n = neq{fJ,,T) = -j—-^ / dqo / d knF{E{qo,k) - fi) 

J-Tv/d J 



(2vr)= 



(42) 



and the current density from Eq. (37) 



with 



Cc<l{^J^,T) 



Jmbt 



2d 



,^K1 I 



eFdh/T 



h {h/Tf + {eFdy 



rCcq(At,T) 



Tx/d 



dqo / d^k cos{qod)nF{E{qo, k) - ^) 



(43) 



(44) 



(27r)3 J_^/^ 

Note that the field dependence of the current density is identical with the simple Esaki-Tsu 
result (34), but the prefactor has a complicated form, which will be analyzed in the following. 
The /e-integration can be performed analytically in neq(/u,T) and Ccq{p,T). 



d knF{E{q,'k) - ix) = po I dEk- 





1 



pQksT log 



( ^i-E(q) \ 



+1 



(45) 



l^^i ^ ^0 e ^bt +1 

Here po = mc/vr/l^ is the density of states of the two-dimensional electron gas parallel to the 
layers including spin degeneracy. 

At first consider the degenerate case T = which holds for low temperatures ksT <^ 2|Tf |. If 
p > E" + 2\T^\ one obtains [96] n = po{p - E"") and 



JUBT 



2/>o|Tf 



■a 1 2 



eFdh/T 



h {h/Tf + (eFdy 



for 



keT « 2|Tf I 
n > 2|rf Ipo 



(46) 



This expression is independent from the carrier density. (For superlattices with very thin barriers 
a different behavior has been reported [97] which was explained by the one-dimensional character 
of the tranport due to inhomogeneities.) A second instructive result is obtained for very low 
densities < (^ - £^")/2|rf | -M <Cl when 

3 



neq « Ccq « Po\Ti\—- 

6tt 



holds and thus 



-'^MBT = en ^ 



eFdh/T 



h {h/Tf + {eFdf 



arccos 



for 



E,,, 



^ 



9 T'* 



ksT « 2|rf I 
n < 2|rf Ipo 



(47) 



(48) 
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which gives J = e{n/d)vET{F) providing the Esaki-Tsu result (34) mentioned above. 
In the non- degenerate case ^u <C -E" — 2|Tf | one obtains [8] 

2\T?\ 



n = pokBTe\''^^Jloi-l^] (49) 

. , ^ (^) ^ f2\Tf\\ 2|rf| eFdh/r 

with the modified Bessel functions Io,Ii, see Eq. (9.6.19) of [72]. For low temperatures ksT <C 
2|Tf I the argument x of the Bessel functions becomes large and loix) ~ Ii{x). Then Eq. (48) is 
recovered again. For high temperatures kBT ^ 2|T"| the Bessel functions behave as Io{x) ~ 1 
and Ii{x) ~ x/2, and 



2|rf|2 eFdh/T 



n <^ poksT 



Such a 1/T dependence of the current density has been observed experimentally in [98,99] albeit 
the superlattices considered there exhibit a rather small miniband width and the justification of 
the miniband transport approach is not straightforward. 

3.1.2. Two scattering times 

A severe problem of the relaxation-time model is the fact that all scattering processes restore 
thermal equilibrium. While this may be correct for phonon scattering, where energy can be 
transferred to the phonon systems, this assumption is clearly wrong for impurity scattering, 
which does not change the energy of the particle. The significance of this distinction can be 
studied by applying the following scattering term [5,100] 

( df{q,k,t) \ ^ nF{E{q,k)-p)-f{q,k,t) ^ f{-q,k,t) - f{q,k,t) 

V dt y^^^.j Te 2rclast 

Here scattering processes which change both momentum and energy are contained in the energy 
scattering time Tg and elastic scattering events, changing only the momentum, are taken into 
account by Tdast- While the Boltzmann equation was solved explicitly in Sec. 3.1.1, dynamical 
equations for the physical quantities of interest are derived here by taking the appropriate 
averages with the distribution function. At first consider the electron density (38). Performing 
the integral ,^,3 J dqd^k of Eq. (35) one obtains 

dn(t) ricq — n{t) 

dt Te 

where the periodicity of f{q) has been used to eliminate the term oc F. Therefore the elec- 
tron density is again given by ncq{p,,T), see Eq. (42). In the same manner one obtains (using 
integration by parts in the term oc F) 

dJ(t) e'^Fd ,^„, , , J(t) , ,, 

^-— 2|Tf|e«) = -J-i (54) 
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(53) 



with the momentum relaxation time l/r^ = l/Tg + l/rdast and the average 

2d l"^!'^ f 
c{t) = j—^ J dqj d\ cosiqd)fiq, k, t) . (55) 

Finally, the dynamical evolution of c(t) is given by 

dc(t) ^ Fd J ^ Ccciip,,T) - c{t) 

dt 2|Tf| Te 

This equation can be considered as a balance equation for the kinetic energy in the superlat- 
tice direction as j^p^ f dqd'^k E{q)f{q,'k,t) = E"'n{t) — 2|Tf|c(t). The stationary solution of 
Eqs. (54,56) gives 

•^^^^ = '^ h (n/reff)^ + (eFd)^ "-'^^^'^^ ^^^) 



with the effective scattering time t^s = y/TeTrn and (5 = y^W^- This is just the result (43) 
with the additional factor 5 reducing the magnitude of the current, as Tm < ^e- The prefactors 
ncq{fi,T) and Ccq(/U, T) are identical with those introduced in the last subsection. 

The relaxation time approximation has proven to be useful for the analysis of experimental 
data by fitting the phenomenological scattering times Te, Tdast- In [101] the times l/rg = 9x lO^'^^/s 
and l/Teiast = 2 X 10^^ /s has been obtained for a variety of highly-doped and strongly-coupled 
superlattices at T = 300 K. 

It should be noted that there is an instructive interpretation [102] of Eq. (57) which may be 
rewritten as 



with the low-field velocity 



where 



J eCeq(/i,T) Uif(F)+t;M(F) ' ^^^^ 



vidF) = '"""^ F (59) 

'^miniband 



1 d'^E'^iq) 2lTfld2 

vy^ _ I II (go) 



"Iminiband ?l^ dq"^ \g^Q k"^ 

This is just the standard expression for the linear conductivity which is dominated by momentum 
relaxation. The high-field velocity 

vm{F) = ^ (61) 

is determined by the maximal energy loss per particle given by Pioss = 2lT{*]/Te for the particular 
scattering term (52) where k is conserved. This relation results from the energy balance providing 
negative differential conductivity as already pointed out in [103]. While both expressions for the 
low-field velocity and the high-field velocity are quite general, it is not clear if the interpolation 
(58) in the form of a generalized Matthiessen's rule [102] holds beyond the relaxation time model. 
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3.1.3. Results for real scattering processes 

The relaxation time approximations discussed above contain several problems: 

- The scattering processes conserve k, which is artificial. An adequate improvement to this point 
has been suggested in [104]. 

- The magnitude of the scattering times is not directly related to physical scattering processes. 

- Energy relaxation is treated in a very crude way by assuming that in-scattering occurs from 
a thermal distribution. 

In [105] balance equations have been derived for the condition of stationary drift velocity and 
stationary mean energy. Here the distribution function was parameterized by a drifted Fermi- 
function similar to the concepts of the hydrodynamic model for semiconductor transport (see 
[106] and references cited therein). This approach allows for taking into account the microscopic 
scattering matrix elements for impurity and electron-phonon scattering and good results were 
obtained for the peak position and peak velocity observed in [83]. 

Self-consistent solutions of the Boltzmann equation have been performed by various groups. 
In [91] results for optical phonon and interface roughness have been presented where Boltz- 
mann's equation was solved using a conjugate gradient algorithm. Using Monte-Carlo methods 
[94] the Boltzmann equation can be solved to a desired degree of numerical accuracy in a rather 
straightforward way (at least in the non-degenerate case and without electron-electron scatter- 
ing). Results have been given in [107] for acoustic phonon scattering and in [108] for optical 
phonon and impurity scattering (using constant matrix elements). Modified scattering rates due 
to collisional broadening have been applied in [109] without significant changes in the result. 
Recently, extensive Monte-Carlo simulations [110,111,92] have been performed where both op- 
tical and acoustic phonon scattering as well as impurity scattering has been considered using 
the microscopic matrix elements. Results of these calculations are presented in Fig. 8a. The 
general shape of the velocity-field relations resembles the Esaki-Tsu result shown in Fig. 5b both 
here and in all other calculations mentioned above. This is demonstrated by a comparison with 
the two-time model (57), where the scattering times have been chosen to give good agreement 
with the Monte-Carlo simulations, see Fig. 8b. The increase of the scattering rate with lattice 
temperature can be attributed to the enhanced phonon occupation. In contrast, the high-field 
behavior does not strongly depend on lattice temperature. Here the drift velocity is limited by 
energy relaxation (61) which is dominated by spontaneous emission of phonons and thus does 
not depend on the thermal occupation of phonon modes. 

3.2. Wannier-Stark hopping 

If a finite electric field is applied to a semiconductor superlattice, the Bloch states are no longer 
eigenstates. Within the restriction to a given miniband u, Wannier-Stark states ^'^^{r,z) = 

<I>^(z) e''^''" /\/]4 with energy E'^j^ = -E^s ~ jeFd + E^ diagonalize the Hamiltonian as discussed 
in Sec. 2.4. (We apply a normalization area A in the (x, y) direction here yielding discrete values of 
k and normalizable states. For practical calculations the continuum limit ^j^ -^ A/{2it)'^ f d'^k 
is applied.) These states are approximately centered around well n = j. In the following we 
restrict ourselves to the lowest band v = a and omit the index z^. In a semiclassical approach the 
occupation of the states is given by the distribution function fj (k) . Scattering causes hopping 
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Fig. 8. (a) Drift- velocity versus electric field for a superlattice with miniband width 20.3 meV and period d — 5.1 nm 
with homogeneous doping Nsd = 10^^ /cm^. The calculation has been performed within the miniband transport 
model employing the Monte-Carlo method (Fig. 3.8 from [92]). (b) Fit by the two-time model (57) using scattering 
times Toiast = 4ps,1.6ps,l.lps,0.12ps and Te = 10ps,2ps,0.9ps,0.4ps for T = 4K,77K,150K,300K, respectively. 

between these states [33,112]. Thus, this approach is caUed Wannier- Stark hopping. Within 
Fermi's golden rule the hopping rate is given by 



R 



j,k^jr',k' 



27r 



(^,VI^ 



scatt I (Rfl \ 



' 6 {Ek' - jeFd -Ek + ieFd [±hupi,anon]) [l - /,(k')] (62) 



where the term [ib/iajphonon] has to be included if emission or absorption of phonons is considered. 
For details regarding the evaluation of scattering matrix elements see [33,92,113]. The current 
through the barrier between the wells m = and m = 1 is then obtained by the sum of all 
transitions between states centered around wells m < and those centered around m > 1, i.e. 

2 (for Spin)e 



JwSH 



A 



■ X] X] X] [^i,k^i,k'/j(k) - i?j,k'^j,k/i(k')] 

i<0 j>l k,k' 



(63) 



If the occupation /j(k) = /(k) is independent of the index i, i.e., the electron distribution is 
homogeneous in superlattice direction, one finds 



2e 
~A 



JwsR = =r 2^ L ^ [«o,k^/»,k'/(k) - i?/.,k'^o,k/(k')] (64) 

?i>lk,k' 

where Ri^]i-tjy = Ro,k^{j-i),k' has been used. Typically, in the evaluation of Eq. (64) thermal 
distribution functions /i(k) = npiE^^g + E^ — /Uj) are employed. The underlying idea is the 
assumption that the scattering rates i?i,k— >i,k' inside each Wannier-Stark state are sufficiently 
fast to restore thermal equilibrium. In this case Eq. (64) can be further simplified to 



■AvSH = ^ 2.^ 2.^ ^-Ro,k^?i,k''raF(-E'ws + -Efc 

/i>lk,k' 



/^) 



1 — exp 



heFd \ 



(65) 
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Evaluating this expression for various types of scattering processes one obtains a drop of the cur- 
rent density with electrical field as shown in Fig. 9. This is caused by the increasing localization 
of the Wannier-Stark functions (see Fig. 4) which reduces the matrix elements (•5j,k'|-f^^'^^*''l^i,k) 
with increasing field. This behavior can be analyzed by expanding the scattering matrix elements 
in terms of Wannier states |^^k) fro™ Eq. (31): 

^) ^-^ [^d] ^*vi^^^'"mk) • (66) 

As the Wannier states are essentially localized to single quantum wells the diagonal parts m = n 
dominate. Neglecting correlations between the matrix elements (^^k'l-^'^'^'^'^l^rnk) ^°^ different 
wells m one obtains: 



RoM^hy = -tYI 



m 



2T7\ (2Tf 



\ m,k' I I m,k 



eFd) \eFd^ 

X 6 {Ek' - heFd - Ef, [±nWp\,onon]) [l - /(k') 



(67) 



For eFd ^ 2Tf the Bessel functions behave as Jn{2x) ~ x^/n! giving a field dependence 
oc (Ti/eFd)"^''. Therefore the transitions 0^1 dominate and 






-1 



/iTfC I ijscatt|,T,a \ I /iTfi I uscatt I ,T,a 
\^0,k'l-W I^O.k) + \^l,k'l^ 1^1, 



(68) 



h \eFd, 

X 6 {Ey -Ek- eFd [±?ia;phonon]) [l - /(k')] • 

For high electric fields the wave vector k' must be large in order to satisfy energy conservation 
and thus the scattering process transfers a large momentum. If the scattering matrix element 
does not strongly depend on momentum (such as deformation potential scattering at acoustic 
phonons) JwsH oc 1/F'^ is found, while different power laws occur for momentum dependent 
matrix elements (such as JwSH oc 1/F^'^ for impurity scattering [113]). For optical phonon 
scattering, resonances can be found at heFd = /iWopt when hopping to states in distance h 
becomes possible under the emission of one optical phonon [110], see also Fig. 9. 

If the sum in Eq. (64) is restricted to h < /imax one obtains a linear increase of the current-field 
relation for low fields [33,112] and a maximum at intermediate fields before the current drops 
with higher fields as discussed above. In [113] it has been shown that this is an artifact and the 
correct J{F) relation is proportional to 1/F for low fields. Thus the linear response region for 
low fields cannot be recovered by the Wannier-Stark hopping approach. 

While most calculations have been performed assuming thermal distribution functions /(k) 
recently self-consistent calculations of the distribution functions have been obtained by solving 
the semiclassical Boltzmann equation for the Wannier-Stark states: 

^^^^ = E ^J,k'^.,k/i(k', t) - i?.,k^,,k'/i(k, t) (69) 

j,k' 

The self-consistent stationary solution of this equation can be used for the evaluation of Eq. (64). 
As can be seen in Fig. 9 significant deviations between both approaches occur for low lattice 
temperatures, when electron heating effects become important. 
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Fig. 9. Drift- velocity versus electric field for a superlattice with miniband width 20.3 meV and period d = 5.f nm 
with homogeneous doping N^d ~ 10^® /cm"'. The transport model of Wannier-Stark- Hopping has been applied. 
Both self-consistent (thick lines) and thermal (thin lines) distributions of electrons in the k direction have been 
used (Fig. 4.20 from [92]). 

3.3. Sequential tunneling 



If the barrier width of a superlattice is large, the structure essentially consists of several decou- 
pled quantum wells. In each quantum well n we have a basis set of wave functions ^5^(2) e /\A, 
where ^^^(-z) is the v eigenfunction of the quantum well potential. The states have the energy 
E'^ + Ek + ecpn, where the potential energy due to an electrical potential (pn has been considered 
separately. The notation is clarified in Fig. 10. If the wells are not coupled to each other (infinite 
barrier width or height) no current flows in the superlattice direction. For finite barrier width 
the states from different wells become coupled to each other which can be described by a tunnel 
matrix element Hn,rn in the spirit of [68] inducing transitions between the wells. In lowest order 
perturbation theory the transition rate is given by Fermi's golden rule 

2tt. 

T 

As the superlattice is assumed to be translational invariant in the (x, y)-plane, the matrix element 
is diagonal in k. Thus, transitions are only possible if E^ — E^ = e(j)m — &<Pn-, suggesting sharp 
resonances when the potential drop between different wells (j)m — <Pn equals the energy spacing 
of the bound states. This could be nicely demonstrated experimentally in Ref. [114] for simple 
superlattices and in Ref. [62] for superlattices with a basis where even tunneling over 5 barriers 
was observed. In the presence of a strong magnetic field along the superlattice direction further 



R{m, ;U, k ^ n, u, k') 



I "m'nl*^)*^ J 



'diE" + Ek' + e(/>„ -Ef'-Ek-e<f 



(70) 
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Fig. 10. Sketch of the levels in neighboring wells n = m + 1, where /i = a is the lowest and v = b is the first 
excited level in the respective quantum wells. The chemical potentials /im as well as the electrical potentials e<j>m 
are marked as well. For a constant electric field F one finds (;/>m+i ~ (j^m ~ Fd. 

peaks due to transitions between different Landau levels are observed [115]. 

The resonance condition from Eq. (70) implies vanishing electric field for transitions between 
equivalent levels (in particular the lowest level). As for zero field the current vanishes (provided 
the electron density is equal in both wells), it was concluded that only phonon-assisted tunneling 
processes are possible. Thus neither a linear increase of the current for low fields nor a peak at 
low fields was expected for weakly coupled superlattices [8]. This conclusion is in contradiction 
with experimental findings [22], where a drift-velocity in qualitative agreement with the Esaki- 
Tsu result (Fig. 5) has been obtained for weakly coupled superlattices as well. This discrepancy 
is due to the neglect of broadening in the argument given aboveJ_] 



3.3.1. General theory 

In a real quantum well the states ^^i.^) ^''' "" / v^ with energy E'^ + E^ + ecpn are not exact 
eigenstates of the full Hamiltonian due to the presence of phonons and nonperiodic impurity 
potentials. The respective scattering processes lead to an energy shift AE and a finite lifetime 
T of the states. These features can be treated within the theory of Green functions (see, e.g., 
Ref. [65]). While a general treatment is postponed to Chapter 4, a motivation of the concept 
and a heuristic derivation of the current formula (79) is given in the following. 

For a stationary fiuctuating potential V{x,y) due to impurities or interface fluctuations one 
finds 



j^/ ft, ft. 



' ^iEk - Ek') 



(71) 
(72) 



® Broadening had been included in earlier theories [73] but there a term was missing which is essential for the 
transition between equivalent levels. This point is discussed in Appendix A. 
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where the second order of stationary perturbation theory as well as Fermi's golden rule was 
applied for a stationary fluctuating potential V{x,y) due to impurities or interface fluctuations. 
For a particle, which is injected at t = 0, the time dependence of the wave function is then given 
by 

^5; ,^(t) = ^^'■^^^(k; t, 0) = -ie(t) e-'(^"'+-^'=+^<^"+^''°')*/'i (73) 

with S''^* = AE{k) - i/i/2r(k), so that |^5^^k(t)P = e"*/^. This motivates the meaning of the 
(retarded) Green function G^^^ and the (retarded) self energy S''''*, which are key quantities in 
the theory of Green functions. (A nice introduction can be found in Ref. [116].) The Fourier 
transformation is given by 

and the spectral function is defined by 

<(k,i?) = -2Im{Gr*(k,i?)} . (75) 

For infinite hfetime r ^ cx) one finds A'^{k, E) -^ 2tt6{E - [E" + E^ + e0„ + AE{k)]) which is 
(except for the factor 27r) just the contribution of the state (n, z^, k) to the total density of states. 
This relation is more general and ^^'^(k, E)/{2'k) can be viewed as the contribution of the state 
(n,z/, k) if the system is probed with an energy E. Here the expressions (71,72) approximately 
correspond to the Born approximation for Y7^^ . Results for higher order approximations as well as 
for phonon or electron-electron scattering can be obtained in a systematic manner by the theory 
of Green functions [65]. Typically 'S^'^^(k,E) becomes a function of energy and momentum, but 
the structure of Eqs. (74,75) persists. 

As the wells are considered to be almost uncoupled, the electronic distribution is described 
by np{E — fj-n — e^n) in each well with the local chemical potential /U„ measured with respect 
to the energy of the lowest level (see Fig. 10). This provides an electron density npiE — Hn — 
e4>n)A'^(M; E)/{27r) per energy in the state (n, z/, k). Summing all states and integrating over 
energy, the electron density in well n (in units 1/cm^) is given by 



nn = ^'^^^^^^ f dEAUKE)nFiE-fi^-e<l>n) (76) 

= Y, j ^E p^^{E)nF{E - ^xn) (77) 

with the density of states 

PI{E) = ^-^^ E <(k' E) = ^r d^^ <(k, E) (78) 

in well n belonging to the level v. Regarding the transitions from level // in well m to level v in 
well n, Eq. (70) is modified as follows: 

- The energy E is conserved instead of the free particle energy E'^ + Ek + ecpn ■ 

- The energy conserving 5- function is replaced by A'j!^(k' ,E)/2tt. 
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- At each energy the net particle flow from weh m to weU n is proportional to (np^E — 
ecprn — IJ'm) [1 " nF{E — ecpn — Mn)]) while the net particle flow from n to m is proportional to 
{np{E — e(j)n — f^n) [1 — nplE — ecpm — f^m)])- Therefore the total current is proportional the 
difference [np{E — ecpm — fJ-m) — npiE — ecpn — fJ-n)] in occupation between both wells. Defin- 
ing the electrochemical potential Epin) = e(j)n + Hn it becomes clear that the current is driven 
by the difference in the electrochemical potential. 

Then the current density from level /x in well m to level i' in well n is given by 

2 (for spin)e 
~2ttKA 



p-*^ ^ ^lior spmje ^p ^^^ 



dEA'^iKE)A';^ik',E) 

(79) 



k,k' 

X [nF{E - e(t)m - ^m) - npiE - e(/)„ - ^„)] . 

While the derivation given above is heuristic, a microscopic derivation is given in Sec. 9.3 of [65]. 
In appendix B.l it will be shown that Eq. (79) is the limiting case of the full quantum transport 
theory based on nonequilibrium Green functions in the limit 2Ti <C T. The same approach has 
been used for tunneling between neighboring two-dimensional electron gases [117,118]. It should 
be noted, that Eq. (79) only holds if the correlations between scattering events in different wells 
are not significant, otherwise disorder vertex corrections must be taken into account [119]. 

An important task is the determination of the matrix elements ^^'knk'' ^^^ possibility is 
to start with eigenfunctions of single quantum wells and consider the overlap between those 
functions obtained for different wells m, n. This is the procedure suggested by Bardeen [68] 
which essentially has been applied in [66,67] for superlattice transport. Another possibility is 
to start from the Wannier states (see section 2.3). If a constant electric field is applied to 
the superlattice structure, the electric potential reads (/>„ = —eFd and the matrix elements 
can be obtained directly from Hi in Eq. (24) or from Hl*^^ in Eq. (28). In subsection 3.3.3 it 
will be shown that the latter Hamiltonian is more appropriate (as already suggested in [73]) 
by a comparison with experimental data. These matrix elements are calculated for a perfect 
superlattice structure. Thus they conserve the parallel momentum k. In addition scattering 
processes at impurities, phonons, or interface fiuctuations may cause transitions between different 
wells. The respective matrix elements can be obtained from the respective scattering potential 
and the Wannier functions as well. Examples will be given in section 3.3.3, where it is shown 
that these processes give a background current while the peaks in the current-field relation is 
typically dominated the k-conserving terms of HY'^. 

An important feature of Eq. (79) is the fact that the current is driven by the difference of 
the electrochemical potential in both wells. This is in contrast to the findings of [73] (based 
on density matrix theory) where the current is driven by the difference /m(k) — /^(k') where 
/^(k) = npiE^ + E^ — fin) is the occupation of the state (n, u, k). In the latter case the current 
vanishes for equivalent levels (/U = i^) if the matrix element H^^.,^ ^, is diagonal in k. In Appendix 
A it will be shown how the factor [np{E — e(j)m — f^m) — np{E — e(j)n — f^n)] is recovered from 
density matrix theory. 

3.3.2. Evaluation for constant broadening 

Here we want to derive simple expressions for the current under the assumption of a constant 
broadening for the electronic states. We assume that only the lowest level ^ = a is occupied (i.e., 
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fim,kBT<^E^). 

First we consider nearest neighbor tunneling with k-conserving matrix elements. Then Eq. (79) 
can be rewritten as follows: 

[npiE - e4>m - fJ-m) - ripiE - e4>m+i - ^m+l)] 
with 

Jo d^fe^m(k,^) 

Here the effective field F = ((/>„ — (pm+ij/d has been introduced and the density of states p'^{E) 
from Eq. (78) has been applied. Let us assume a constant self-energy S''^'^*(k, i?) = — ir'^/2 in 
Eq. (81) for the sake of simplicity. Then the spectral functions become Lorentzians ^5:^(k, £■) = 
V'/liE - Ek- E" - e(j)mf + (r''/2)2]. As p^{E) is essentially zero for energies below E" we 
may restrict ourselves to E > E"' in the evaluation of Eq. (81). For these energies the function 
A'^^-^ (k, E) takes its maximum at E^ > eFd + E"" — E'^ which is larger than zero provided 
we restrict ourselves to eFd > E^ — E'^ (this is always the case for the a ^ a resonance and 
eFd > 0). In this case the integrand of Eq. (81) does not take large values for Ek < and it is 
justified to extend the lower limit of integration to — oo. A straightforward evaluation using the 
calculus of residues yields: 

(^™+i ) = 7 1^ 7 WTT^ with r'^''''^ = F'' + r'^ (82) 

which only depends on F. Numerical evaluations indicate that these simplifications are quite 
good for fields above the resonance, i.e. eFd > E'^ — E"". In the case oi f ^ a and eFd < 
E^ — E"" the choice F'^''^''^ = F*^ may be better. Note that this simple model with a constant 
self-energy cannot be used in the calculation of the density of states (78) as the integral for 
the electron density (77) diverges. Therefore we use the free-electron density of states pra{E) = 
PqQ{E — E"" — ecpm) in the following. With these simplifications Eq. (80) becomes: 

ja^u I m+1, ml J- / p r \ fQ'J\ 

Jrn^rn+1 " « ^ ,^^ - ^a - eFdY + {T^,^^ 12)'^'^''^^ n„, n^+ij \p6) 

with the effective electron density (here E = E — e(f>m) 

POD 

ncs{F, rim, rim+i) = dE po npiE - prn) - npiE - prn+i + eFd) (84) 

which describes the difference of occupation in both wells. The electron density is related to the 
chemical potential via: 



n™, = / dE pon F {E - prn) = poksT log 1 + exp f "^ j 



(85) 



Here we assumed that only the lowest level is occupied. The inclusion of higher levels is straight- 
forward, but leads to more complicated expressions. Inserting into Eq. (84) yields: 
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ncs{F,nm,nm+i) = rim - PoksTlog 






eFd 
'I^bT +1 



(86) 



In the nondegenerate limit limit we have n^+i ^ PoksT and Eq. (86) is further simphfied: 

eFd 



nesiF, rim, rim+i) =nm- rim+i exp 



ksT 



(87) 



It is interesting to note that the total current can be written as a discrete version of the drift- 
diffusion model in this case: 



Jm-^m+1 = e—v{F) - eD{F) — 



with the velocity 



„,(p)-^^\fj'^,a |2 

V 

and the diffusion coefficient 



-•v,e?i 



{E^ -E"- - eFdY + (r^'^ff/2)2 



^v,e?i 



eFd\ 



{E'' -E^- eFdf + (r^'^ff/2)2 



exp 



eFd 



(89) 



(90) 



which satisfy the Einstein relation D(F) = v'{F)kBT/e for F = 0. Remember that these ex- 
pressions only hold for eFd > 0. For eFd < the expressions for the opposite direction can be 
applied. 

In the last part of this subsection some explicit results for the current a ^ a are given. Here 
we assume equal densities Um = rim+i = n in the wells and set -ff^+i m ~ Ff according to 
Eq. (24). 

In the degenerate limit {n » poksF) we find n^s = n — {n — poeFd)@{n — poeFd). For 
eFd < n/po we obtain 



^ _, iTfP 2T-eFd 



h {eFdf + (r«)2 
In the nondegenerate limit we obtain for eFd ^ k^T 
, ,^, en iTfp 2r«eFd 



for 



for 



n » poksT 
n > PoeFd 

n ^ poksT 
ksT > eFd 



(91) 



kBT h {eFdY + (r«)2 
In both cases the field dependence is given by the Esaki-Tsu result (34) with F" = H/t. 



(92) 



3.3.3. Results 

A variety of calculations for different superlattice structures have been performed within this 
model. These calculations consist of the following steps: 
(i) Calculation of the miniband structure E^{q) and the associated wave functions (pq{z) ac- 
cording to Sec. 2.1. 
(ii) Evaluation of the Wannier-functions and respective couplings Tj^,R'^''^, see Sec. 2.3. 
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eFd [meV] 

Fig. 11. Current-field relation for the superlattice studied in [123] for constant electron density n-m = Wm+i = No- 
The dashed line shows the current from resonant transitions with momentum independent matrix elements, while 
both contributions from resonant and non-resonant currents contribute to the full line. An electron temperature 
Te = 20 K was used in the calculation. 



(iii) Evaluation of (intrawell) scattering matrix elements (^J^k'l-^'^'^ l^nk) ^"^^ ^^^ dominant 

scattering processes. For doped superlattices ionized impurity scattering dominates and the 

respective calculations including screening are presented in [120,35]. Scattering processes 

at interface roughness [35] and phonons have been considered as well. 

(iv) Calculation of the self-energies Y7^^^{\i.,E) within different approximation schemes such as 

the self-consistent single-site approximation [120-122] for impurity scattering, 
(v) Determination of the chemical potential for given electron density provided by the donors 

from Eq. (77). 
(vi) Renormalizing of the matrix elements for each electric field according to Eq. (29). 
(vii) Evaluation of the current density according to Eq. (79). 

Results for the superlattice studied experimentally in [123,124] (9 nm wide GaAs wells, 4 
nm AlAs barriers, doping density N£) = 1.5 x 10^"*^ /cm^, cross section A = 1.13 x 10~'^cm^) 
are shown in Fig. 11. Two pronounced peaks can be identified at low and moderate fields. The 
low-field peak is due to tunneling between the lowest levels (a -^ a), while the peak around 
eFd ^ E — E^ = 130 meV is due to a ^ 6 tunneling. While for the dashed line only k- 
conserving transitions are taken into account, the full line also includes the contributions of 
interwell scattering matrix elements {"^^+1 k' I -f^^'^^** I ^n k) evaluated for interface roughness, see 
[35]. (The respective matrix elements for impurity scattering are negligible.) These scattering 
events represent an additional current channel in Eq. (79) yielding a background current which 
dominates between the resonances, but is negligible compared to the resonant currents. The 
height of the a ^ b peak (/max = 1-75 mA) is in good agreement with experimental data 
exhibiting /max = 1-45 mA (Fig. 6 of [124]). The low-field peak is not resolved experimentally 
due to domain formation yielding a current of about 0.076 mA. 

Results for the superlattice studied in [31] (15 nm wide GaAs wells, 5 nm Alo.3Gao.7As barriers, 
doping density Nd = 6 x 10^/cm^, cross section A = 8/im^) are shown in Fig. 12. The general 
shape is like the results from the sample mentioned before. While the latter (highly-doped) 
sample did not exhibit a strong temperature dependence, the situation is different for the low- 
doped sample considered here. For low electron temperatures Tg < 15 K the electrons are located 
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Fig. 12. Current-field relation for the superlattice studied in [31] for constant electron density rim ~ nm+i = No- 
The bias was taken to be NFd for the theoretical results, where Ai" = 10 is the number of quantum wells (from 

[120]). 

in the impurity band (about 10 meV below the free electron states). Therefore a current exhibits 
a peak at eF^ighd ~ 10 meV when these electrons can tunnel into the free electron states of the 
next well. For higher temperatures the electrons occupy the free electron states and the maximum 
occurs at eFio^d ss F" « 2 meV as suggested by Eq. (92). Due to the same effect the peak at 
eFd ~ 50 meV due to a ^ 6 tunneling shifts with temperature. The experimental data (taken 
at a lattice temperature of 4K) are shown for comparison. While the low-field conductance is 
in good agreement with the Tg = 4 K calculations, close to the first maximum the agreement 
becomes better for the Te = 35 K curve, which can be caused by electron heating. The heights 
of both maxima are in excellent agreement between theory and experiment. The difference in 
position of the second maxima may be caused by an additional voltage drop in the contact, 
which is not taken into account in the theory, where the field was just multiplied by the length 
of the sample. Finally, the saw-tooth shape of the experimental current-voltage characteristic is 
due to the formation of electric field domains as discussed in section 5. 

In Fig. 13 results are shown for the superlattice structure from Ref. [125] (25 nm wide GaAs 
wells, 10 nm Alo.3Gao.7As barriers, doping density iV/j = 1.75 x 10^'^/cm^). Due to the large 
well width the level seperation is small and several resonances -E" -^ E^ can be observed with 
increasing field. The calculations (for T = 4 K and within the approximation (83) applying 
phenomenological broadenings V^ = 4 meV for all levels) have been performed with and without 
taking into account the renormalization of the matrix elements according to Eq. (29) including 
the lowest 6 levels. Fig. 13 shows that the result with renormalized matrix elements (full line), 
see Eq. (28), is in good agreement with the experimental result but exhibits higher peak currents. 
This may be due to an overestimation of the couplings in the calculation. E.g. assuming barriers 
of 10.6 nm, the current drops by a factor of 2. An increase of the Al-content in the barriers 
would give a similar trend. The results with the bare matrix elements (dashed line), see Eq. (24), 
deviates strongly from the experimental result. In particular, the peak currents do not increase 
for resonances at higher fields. This shows that the renormalization procedure is essential if 
higher resonances are considered. 
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U, NFd [V] 

Fig. 13. Current-field relation for the superlattice studied in [125]. The bias was taken to be NFd for the theoretical 
results, where A'^ = 30 is the number of quantum wells. Full line: Result from sequential tunneling with renormalized 
matrix elements. Dashed line: Result from sequential tunneling with bare matrix elements. Dots: Experimental 
data (courtesy of Yu. A. Mityagin). 

3.3.4. Tunneling over several barriers 

Up to now the discussion was restricted to next-neighbor couphng, which is described by 
matrix elements H^'^ in Eq. (28). The extension to tunnehng over h barriers can be treated 
analogously taking into account the matrix element Hj^'^ . The discussion of Sec. 3.3.2 can be 
performed in the same way applying the bias drop heFd. Thus we expect resonances at field 
strengths eFd = (E'^ — E"')/h. The total current is then given by 



•-max ^max 



JsT = E E ^J^ 



(93) 



h=l u=l 



where the individual current densities Jg— ^^ are evaluated according to Eq. (79). Here it is 
assumed that only the lowest level a is occupied. For the samples discussed in the last section, 
as well as for most other samples considered, the respective currents for h > 1 are negligible. 
This is different for the sample discussed in [126] (5 nm wide GaAs wells, 8 nm Alo.29Gao.71 As 
barriers, doping density Nd = 2.25 x 10^^/cm^), where the second miniband is located around 
the conduction band of the barrier and the subsequent minibands zv > 3 resemble free particle 
states. Results of the calculation with /imax = 2, i.e. taking into account tunneling to next- 
nearest neighbor wells, are shown in Fig. 14. If the calculation is restricted to the two lowest 
levels (t'max = 2) the current-field relation resembles the findings of Fig. 11. There is a peak at 
low fields due to a ^ a tunneling and a peak at eFd ^ E^ — E"- = 0.177 meV due to tunneling 
a -^ &(1) into nearest neighbor well. The matrix element Hif is small, thus no transitions 
a -^ 6(2) to the next nearest neighbor well can be seen. This changes completely if the third 
level (t'max = 3) is taken into account for the renormalization of the energy levels and couplings. 
First the strong coupling to the third level diminishes the value of E'' by 15 meV close to the 
a -^ 6(1) resonance. Thus the position of this resonance is shifted to eFd = 162 meV where the 
new resonance condition is fulfilled. Secondly a new peak arises at eFd = 85 meV~ {E'' — E°')/2 
due to next nearest well tunneling a -^ 6(2) (remember that the renormalized level energies E'' 
are field dependent and thus the local field must be taken into account at each comparison). 
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Fig. 14. Current-field relation for the superlattice studied in [126] for constant electron density rim ~ No- Transi- 
tions between nearest and next nearest wells have been taken into account and the number of levels was i^max ~ 2 
(full line), i^max ~ 3 (dashed line), and i^max = 4 (dotted line). T — 20 K was used in the calculation. 

The reason is the strong admixture of i?2 (which is quite large for the superlattice structure 
considered) in the renormalization of the matrix element i^|^ due to Eq. (29). For the same 
reason a third peak appears at eFd = 136 meV~ {E^ — E"')/2 due to a ^ c(2) tunneling. 
If the fourth level (fmax = 4) is taken into account as well, the result is hardly changed, thus 
providing confidence into the results. These findings are in agreement with the experiments [126] 
where a strong increase of the current was observed at field strengths of eFd ~ 80 meV and 
the current density becomes larger than 0.15 kA/cm^. Nevertheless, no current peak has been 
observed so far in this field region. Tunneling over more than one barrier has also been observed 
experimentally in [62,127]. Current peaks at eFd = {E — E^)/h corresponding to resonances 
between h next neighbors have also been found in the calculation by Zhao and Hone [128]. The 
hight of these peaks was quite small, probably due to the neglect of interband couplings i?^ in 
their calculation. 

These findings show that next-nearest neighbor tunneling is possible in superlattice structures. 
Nevertheless the quantitative description is still an open issue. The inclusion of results from Zener 
tunneling [129] may be helpful in future research here. 



3.4. Comparison of the approaches 

Let us now compare the results from the different approaches miniband transport (MET), 
sequential tunneling (ST), and Wannier-Stark hopping (WSH). 

MBT-ST: Comparing Figs. 8 and 11 one notices that the global behavior with linear increase 
of the current for low fields and a maximum at moderate fields is in qualitative agreement 
for the MBT and ST approach. While the current scales with the square of the coupling for 
ST, the Esaki-Tsu drift velocity is proportional to Ti. This discrepancy is resolved if either 
the temperature or the electron density is high and a jTip dependence of the current density 
is recovered for MBT as well, see Eqs. (46,51). Comparing these results with Eqs. (91,92) we 
find that the simplified expression become identical for MBT and ST if either ksT or n/pQ 
are large with respect to both 2Tf and eFd. This explains the 1/T dependence of the current 
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density observed experimentally in [98,99] for superlattices exhibiting a rather small coupling 
Tf « 1 meV. As the experiments are performed at T > 77 K, the estimations (51,92) hold 
simultaneously and the findings cannot be taken as a manifestation of miniband transport. 

ST-WSH: Both approaches exhibit negative differential conductivity for high electric fields. 
Let us restrict ourselves to the a ^ a resonance and consider a superlattice with nearest neighbor 
coupling Tf. For large electric fields eFd > 2Im{E"''*'*} = T the term A^{k,E)Af{k,E) from 
Eq. (79) exhibits a two-peak structure 

21m {T,1'''\k, Ek + E"")} 



A-,{k,E)AUKE) 



2Tr6{E 
2Im{Sg 



Eu - ^"1 



+ 



\k,Ek + E'' 



{eFdy 
-eFd)] 



(94) 



{eFdy 
Within the Born approximation for the scattering 



27r5(^ -Ek-E" + eFd) . 



Im{srHk,i?)}=Ef k^VI^ 



scattig^a \ 



5{E -E" - Ek' + meFd) 



(95) 



Applying the approximation (68) for eFd 3> 2T" the expression (79) for the current in the ST 
model becomes after several lines of algebra 



la- 
^0- 



2e 
hA 



yZ -Ro,k^l,k'«F(-Efc + Ea- fl) 



k,k' 



1 — exp 



eFd 



(96) 



This is the dominating term of the current for Wannier-Stark hopping (65). Thus, the expressions 
of ST and WSH become identical in the limit of eFd ^ T and eFd ^ 2|Ti|. This is just the 
overlapping region between the ranges of validity of both approaches as depicted in Fig. 7. 

The transition between WSH and MBT is even more difficult. MBT typically exhibits a \/F 
behavior for eFd ^ F as predicted for the Esaki-Tsu relation, while l/F"^ with various exponents 
r > 2 is found for eFd ^ 2|Tf | from the WSH model, see Sec. 3.2. As mentioned there, a 1/F 
behavior can be recovered from the WSH approach by summing all contribution h in Eq. (64) 
for eFd <C 2|rf |. In [110,113] it is shown that in the field region F <C eFd < 2|rf | the results 
of both approaches agree fairly well. Again this agrees with the joint range of validity depicted 
in Fig. 7. 
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4. Quantum transport 

In semiconductor superlattices the miniband width A, the potential drop per period eFd, 
and the scattering induced broadening T are often of comparable magnitudes. This requires the 
application of a consistent quantum transport theory combining scattering and the quantum 
mechanical temporal evolution. Different formulations applying nonequilibrium Green functions 
[130], density matrix theory [131], the master-equation approach [132], or Wigner functions [133] 
have recently been used to tackle this general problem for a variety of different model structures. 

Here the formalism of nonequilibrium Green functions is applied to study electrical transport in 
superlattices. This approach allows for a systematic study of both quantum effects and scattering 
to arbitrary order of perturbation theory. Although the calculations involved are quite tedious 
(as well as the acquaintance with this method) such calculations are of importance for two 
purposes: On one hand it is possible to derive simpler expressions like those studied in the 
preceding section from a general theory, thus shedding light into the question of applicability. 
On the other hand there are situations where no simple theory exists and thus one has to pay 
the price to work with a more elaborate formalism. 

A variety of different quantum transport calculations for semiconductor superlattices have 
been reported in the literature: In [134] an analysis within the density matrix theory has been 
presented, which was simplified to different approaches for low, medium, and high field. The 
same method was applied to study transport in a perpendicular magnetic field [135]. A similar 
approach was performed in [136], where the quantum kinetic approach was solved in the limit 
of Wannier Stark hopping and the nature of phonon resonances were analyzed. The formation 
of Landau levels in a longitudinal magnetic field causes additional resonances [137]. A transport 
model based on nonequilibrium Greens function [138] has been proposed as well, although explicit 
calculations could only be performed in the high temperature limit and within hopping between 
next neighbor Wannier-Stark states there. 

This section is organized as follows: At first the general formalism of nonequilibrium Green 
functions for stationary transport is briefly reviewed in a form which can be applied to a variety of 
devices. The special notation to consider transport in homogeneous semiconductor superlattices 
as well as the approximations used are described in the second subsection. In the third subsection 
the standard approaches (miniband transport, Wannier-Stark hopping, and sequential tunneling 
as discussed in section 3) will be explicitly derived as limiting cases of the quantum transport 
model. This proves the regions of validity given in Fig. 7. Finally, in the fourth subsection results 
are presented for different samples. The results from the self-consistent quantum transport model 
will be compared with simpler calculations within the standard approaches applying identical 
sample parameters. This will demonstrate that the standard approaches work quantitatively well 
in their respective range of applicability. The reader who is less interested in the theoretical 
concept and underlying equations may skip subsections 4.1-4.3 and continue with the results in 
subsection 4.4. 
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4.1. Nonequilibrium Green functions applied to stationary transport 

In this subsection the underlying theory of nonequilibrium Green functions is briefly reviewed. 
The notation of [39,65] is fohowed here and the reader is referred to these textbooks for a detailed 
study as weh as for proofs of several properties addressed here. 

We consider a set of one-particle basis states \a) and aa(t),aa(t) are the corresponding anni- 
hilation and creation operators. The time dependence stems from the Heisenberg picture. Most 
physical one-particle observables can be expressed by the one-particle density matrix 

PaAt) = (4(i)o/3(t)) = Tr [pai{t)af,{t)} (97) 

which is the corresponding quantum mechanical expectation value with the density operator p. 
In particular, the occupation of the state |a) is given by Pa,a{t)- The task of any many-particle 
quantum theory is the evaluation of pap (t) in the presence of a Hamiltonian 

H = Ho + U + iJscatt (98) 

where 

Ho = Y,Eaai{t)a^{t) (99) 



is diagonal in the basis \a), 



U = Y,Uc.AtWa{t)ap{t) (100) 

a,/3 



describes an additional potential term, and -ffgcatt refers to interactions with phonons, random 
impurity potentials (which are treated within impurity averaging), or interactions between the 
particles. 

Within density matrix theory the temporal evolution of Pa,/3{t) is studied directly by applying 
Heisenberg's equation of motion for the product aa{t)af^{t). E.g., the time dependence of the 
occupation Pa,a{^) is given by: 

-^^{ai{t)a^{t)) = '-{[H,al{t)a^{t)]) 

^r^ i r + + -^ . 1 (101) 

= 2Zn [f^/3,a(4(i)aa(i)) - f^a,/3(4(i)a/3(i)) + ([^scatt,aL(i)aa(t)]) 

where [a, h] = ab — ha denotes the commutator. In order to satisfy the equation of continuity, 
the particle currents j{t) between the basic states have to be identified by 

/-"(t) = ^Re [iUp^^{a\{t)a^{t))}+jiZ'^{t) (102) 

which satisfy j^^°'{t) = —j°''^^{t) as each part of -ff is a Hermitian operator. The scattering 
induced current iscatt'(0 can be determined once ^scatt is specified. This term typically con- 
tains higher order density matrices like (aL(t)al(t)aQ,(t)aQ,'(t)) in the case of electron-electron 
scattering. Thus there is no closed set of dynamical equations for Pa,i3{t) and the dynamical 
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evolution generates a hierarchy of many-particle density matrices, which has to be closed by 
approximations, see, e.g., [139,140] as well as references cited therein. 

A conceptually different approach to many-particle physics constitutes the theory of nonequi- 
librium Green functions which has been developed by Kadanoff and Baym [141] and indepen- 
dently by Keldysh [142]. In this theory the time dependence of aa{ti) and a^(t2) is considered 
separately, thus two different times appear in the calculation. The corresponding generalization 
of the density matrix is the correlation function (or 'lesser' Green function) 

G^,,a,ih,t2)=i{ai^{t2)aaAti)) (103) 

which describes the occupation of the states together with the respective correlations both in 
time and state index. Note the unusual order of indices, which will be helpful in later stages of 
the theory. The notation follows [39,65] here. Sometimes (e.g., [141,143]) the factor i is dropped 
so that C" agrees directly with the density matrix for equal times. Next to this correlation 
function the retarded and advanced Green functions are defined by 

GZc.,ih,t2) = -ieih-t2){{aaAh),al^it2)}) (104) 

G^%,{h,h)=ie{t2-h){[a^Ah),aUt2)}) = [GZaAt2,ti)]* (105) 

respectively, where {a, b} = ab+ba denotes the anticommutator which is appropriate for fermion 
operators a^ considered here. These functions describe the response of the system at time ti in 
state «! which is excited at time t2 in state a2- 

4.1.1. Temporal evolution 

The time dependence of these Green functions is given by the following set equations 

''i^h— -Eqi j 001,02(^1 '*2) - X.^«l./3(^l)'^fta:j(^l'^2) 

i^Q^ -^"2 JGo J, 02(^1,^2) - / ,G'^^g{ti,t2)Ui3,a2{t2) 

= E /f [^^t/3(*l'*)SA.2(*'*2) + G<^,^(tl,t)S^^^2(*.i2) 

i/i J- - i?.,) G-*/„f ^(ti, t2) - E Uo.,,p{h)G';';if''{h,t2) 

= h5{t, - 12)5„„„2 + E / y s^T^(*i' *)G?.f ^(*' *2) 
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(106) 



(107) 



(108) 






(109) 



which are derived in section 5 of [39]. The same result is obtained from the matrix equations 
(3.7.5) and (3.7.6) of [65] if the definitions of the retarded and advanced Green functions are 
inserted. 

While the equations for G^^^ and G exhibit a 5{ti — ^2) inhomogeneity typical for Green 
functions, this is not the case for G^ which is, strictly speaking, not a Green function (although 
this term is often used). The self-energies S describe the influence of scattering (compare the 
simplified description in Sec. 3.3.1). They can be expressed by functionals of the Green functions 
which depend on the approximation used, such as the self-consistent Born approximation. Here 
it is crucial to pay attention to the fact that S'''^*, S*^*^^, and S^ belong to the same quantity 
S (a matrix [65] or a self-energy defined on the complex temporal contour [39]) and thus the 
same approximation has to be performed for each quantity. In this way one obtains a closed set 
of integro-differential equations, which governs the temporal evolution of the Green functions. 
On the right-hand side of Eqs. (106,107) the lesser Green functions and lesser self-energies only 
depend on the time arguments (t'l, ^2) with t'^ < ti and t2 < ^2 due to the properties of retarded 
and advanced Green functions. Therefore these differential equations for G^{ti,t2) can (at least 
in principle) be solved explicitly by forward integration. 

Equations (106,107) are solved by the integral equation (sometimes called Keldysh relation) 



G<„.,{h,h) = E / Y [^GZp{i^,t)J:<^,{t,t')Gpl^{t',t2) (110) 



which has a nice interpretation: S5o,(t,t') can be considered as an in-scattering term, which 
creates a correlated one-particle excitation at times (t, t') as a result of a scattering event. The 
retarded and advanced Green functions provide the action of this excitation at the later times ti 
and t2) at which the correlation function G^ is observed. The relation Eq. (110) is a particular 
solution of the differential equations (106,107). The general solution contains a further term 
(proportional to the free evolution of G^ without scattering) to satisfy initial conditions, see, 
e.g., Eq. (5.11) of [39]. Typically, the contribution of these terms decays in time if scattering is 
present, so that Eq. (110) holds in the long time limit. 

If we consider a stationary state without any time dependence of the external potential [/, 
all functions depend only on the time difference ti — t2 and it is convenient to work in Fourier 
space defined by 

Fa„a,{E) = ^j Ate'^'I^F^,^^,{t + t2,t2) (111) 

Fe.„aAtiM) = ^ / d^ e-^(*^-*2)/^F«,,«,(E) (112) 
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both for self-energies and Green functions]^. Then the foUowing relations hold: 

{G^:^p{E)y = GZ{E) and G<^(i^) = - {G<Ji^)}* (113) 

Eqs. (106,107,108,109) yield 

{E - Ea^) G^^^^^{E)- 22 UauisGp^^^iE) 

P 
= E [^tAE)Gto.,{E) + J:<^^^{E)G^%iE) 

{E - E^,)G<^^^^{E)- Y, Gl^^{E)Ufs,a,{E) 

= E [GT.dE)^,c.,{E) + G<^^p{E)T.Z^{E) 
P 
{E - E^,)GtJ,^^^E)- Y, U^,,pGf^if''{E) 

= 6^,,^, + Y Kl[t{E)G-^lf'{E) 



(114) 



(115) 



(116) 



P 

■,ret/adv / 



(117) 



{E - E^,) GtJ,^^^{E)- J2 gTJ {E)Up,., 

P 

P 

If Uf^^a and ^^aj^^ ^ {E) are symmetric matrices, then G'^ai,a2^{E) = G^a2,ai^ (E) holds as well. 
This can be shown by subtracting Eq. (117) from Eq. (116), where ai and 02 are exchanged. 
In the same way the Keldysh relation becomes 

G^r,aAE) = E GZpiE)^^AE)Gp:a2iE) . (118) 

P,P' 

A quite elementary derivation of Eqs. (116,118) is given in section 8 of [38]. From G'^{E) the 
density matrix can be evaluated directly via 

/J J? 
-^Gl^[E) (119) 

which provides us with the one-particle expectation values for most quantities of interest. 

4.1.2. Self-energies 

The self-energies can be obtained from the usual diagrammatic rules which are derived in most 
textbooks on many-particle theory such as [65]. In the following the results are given within the 
self-consistent Born-approximation. 



^ Different definitions have been suggested which produce gauge invariant equations when U is due to a com- 
bination of an electric and magnetic field [144,39]. This is important for various approximations to treat slowly 
varying fields. 
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For impurity scattering one considers an impurity potential V{f; {fi}) which depends on the 
locations fj of the impurities. The respective Hamiltonian is given by 



/3i,/32 



(120) 



with V/3-^^^/3^{{fi}) = {l3i\V{f; {fi})\P2)- For large systems the average {V . . . y)imp over all possible 
impurity configurations {f^} has to be taken (interface roughness can be treated in a similar way). 
Within the self-consistent Born-approximation only correlations between two scattering matrix 
elements are taken into account. Thus one finds 

/3i,/32 

For phonon scattering the respective Hamiltonian reads 



(121) 



/3l,/?2 



fl22) 



-^ phonon / ^ 

where h^^i, b^^ are the (bosonic) annihilation and creation operators of the phonon mode I (such 
as acoustic/optical or longitudinal/transverse) with wave vector p. Using Langreth rules (section 
4.3 of [39]), one obtains the retarded self-energy within the self-consistent Born approximation: 



K,a,iE) =i^ Y. M^,,pMl)Mp,,aMl) j 



P,l /3i,/32 
~fret 



cLEi 

"2^ 



Gflp^{E-E,)Dr{p,E_ 



+ Gf^^p^{E - E,)D<{p, E,) + G< ,^,(i? - E,)Drip, Ei) 



(123) 



where Di(p, Ei) refers to the phonon Green function. Now we replace the phonon Green function 
by its unperturbed equilibrium values (section 4.3 of [39] and Eq. (2.9.9)F] of [65]). 



Dr\p,E^ 



El -huji{p)+ iO+ El + fyjji{p) + iO+ 
1 1 



El - hji{p) - \{)+ El + huJi{p) - iO^ 



(125) 



D<\p,Ei) = - 2T:\{nB{nwimKEi - ^lip)) + [nB{huJi{p)) + 1] 6{Ei + huji{p})} , (126) 



i.e., D<^{p,Ei) = nB{Ei)[Df^^{p,Ei) - Df^'^{p,Ei)] where UBiE) = [exp{E/kBT) - 1]-^ is 
the Bose distribution. As G^^*^ {E — Ei) has only poles for Imjii^i} > 0, we find 

'•^ d^i 



27r 



-Gflp,{E - Ei)Dr\p, El) = -i [GZp,{E - f^m) - GZp,{E + h^m)] 



from the residua of the contour over the complex plane with Im{£'i} < 0. Putting things together 
we obtain: 



Note the sign error for Df^ in Eq. (2.9.9) of [65] 
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P,l PiA 



[nBinu;iip)) + 1] G^t/32(^ - ^liP)) 



+ nB(M(p1)G^t/32(^ + f^iiP)) + lGl^f^,(E - fu^iip)) - \g<^^p^{E + h^i{p)) 



.i/4|G<,,(.-.,(p{^ 



hwi{p) 



V 



E\ +hoJi{p) 



(127) 



The lesser self-energy reads 

/J rp 
P,(. m,P2 
= E E M^,^f3,{p,l)Mp,^^,{p,l)[nB{nwi{p))G^^^p^{E - hwiip)) 

P,l ft,/32 

+ [nB{hu:m) + 1] G<^^p,^{E + M(p1) 



(128) 



which describes the in-scattering from correlated states [3i , [32 by phonon absorption as well as 
stimulated and spontaneous emission of phonons. 

The Coulomb interaction can be easily included within the Hartree-Fock approximation, which 
provides an additional potential (depending on J dEG'^{E)), see, e.g., chapter 8 of Ref. [38]. 
Higher order approximations (describing electron-electron scattering) are difficult to implement. 
All these effects have been neglected in this work. 

The combination of these functionals for the self-energies with Eqs. (116,118) for the Green 
functions provides a coupled set of equations which has to be solved self-consistently to obtain 
the functions G^ (E) for the stationary state. Afterwards the physical quantities of interest can 
be evaluated by Eq. (119). 

4.1.3. Thermal equilibrium 

In thermal equilibrium the electron distribution is governed by the Fermi distribution. In the 
language of Green functions this can be written as 

G<^^^^{E)=inF{E)A^,,^,{E) 

with the spectral function 

Aaua2iE) = i G'^^^a^iE) - G'^^^^^{E) 



(129) 



(130) 

which is derived in section 3.7 of [65], e.g.. As discussed before, the occupation of the state a is 
given by 



Pa,, 



^ j dEG<,{E) =1 dE l.A^^^{E)nF{E) . 



USD 



Thus Aa,a{E)/2Tr represents the energy resolved density of state a. Here it is crucial to note 
the difference to the classical value of the occupation np{Ea), which is only recovered in the 
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free-particle case U = T, = when Aa^aiE) = 2tt6{E — Ea) holds. In contrast pa,a and npiEa) 
will in general differ, if scattering induced broadening leads to a finite width of the spectral 
function. 

These effects can be estimated assuming [/ = and Y7^^^^{E) f» —\Q{E)5ax^a2^/'^ which 
mimics the fact that there are no scattering states below a band edge sA. E = Q. Then one finds 
from Eqs. (116,130): 

r 

and in the limit E^ ^ k^T one obtains 

CksT T 

p^.^ - 2vr {E^y + ryi ^^^^^ 

for a non-degenerate distribution npiE) ~ C e~ ' ^ . Thus, the occupation of the high energy 
states is larger than one would estimate from a semiclassical distribution np{Ea) ~ (j Q-^a/kBT^ 
Finally, it should be pointed out, that the different Green functions C", G^'^^, and G are 
related to each other in thermal equilibrium which allows for a description in terms of a single 
Green function. Thus, the theory of equilibrium Green functions is significantly simpler than its 
nonequilibrium counterpart discussed here. 

4.1.4. Spatial boundary conditions and contacts 

Although we are concerned with homogeneous infinite systems in this section, some remarks 
concerning boundary conditions in real structures are appropriate. They are needed in the dis- 
cussion of transmission through superlattices [145]. 

In order to solve the system of equations discussed above in a finite system, boundary con- 
ditions have to be specified. Here two types can be distinguished. On the one hand there are 
regions where the device is terminated by an insulating layer. Here it is appropriate to neglect 
states in these regions, as their energy is significantly larger than the relevant energies in the 
device. A far more interesting point is the treatment of contacts, which act as a source or drain 
for the electric current. 

We separate the system into a central region with index C and lead region with index L. The 
matrix Ga,/3 can than be divided in submatrices of the type Gql, where the index a belongs to 
the central region and f3 to one of the leads. Then the matrix equation (116) can be written in 
the form 

'G^SciE)G^Sl{Ey 

Gtc{E) GtdE)^ 

^ . Ucc Ucl\ _ Ig^Sc{e) g^SUe) \ 

[ulc Ull) ' {GTciE) GTdE)) 

^^ShiE) ^dliE)\ _ (g^SUe) G^SUE) 

J:Tc{E) ^tdE)) ' \Gtc{E) GfdE), 
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Now let us assume that 'Scl = "^LC = 0, i.e., there is no scattering between the lead regions 
and the central region. Then we can write: 

{E-El + iO^)Grc{E) = Vlc ■ G^SciE) + \^ll ■ GTc{E) + ^TdE) ■ Gtc{E) (135) 

This equation is solved by 

Gtc{E) = GToiE) • Vlc • G^SciE) (136) 

where G^^q (E) is the Green function of the lead satisfying the equation 

iE-BL + iO+) • GToiE) = 1 + Vll- G^iE) + S^KE) • GToiE) . (137) 

It is important to note, that G^iE) is not exactly the Green function of the pure lead. In fact 
Y^TiE) has to be evaluated from the full Green function G'"*^*(£') and not only from G^iE). 
This difference vanishes under the usual assumption that scattering is negligible in the leads. 
Now the part of equation (134) for Gqq{E) can be written as 

{E-Bc + iO+) • GTciE) =1 + Vcc ■ G^ciE) 

+ i^TciE) + Vol ■ GTiE) • U^c] • G^ciE) . ^ ^ 

This is a closed equation for the matrix G^^(-E') concerning the states inside the structure. The 
term TJcL • G^iE) • Vlc can be viewed as an additional self-energy, due to the transitions 
between the central region and the lead. 
In a similar way Eq. (134) yields: 

GTiE) = G%{E) . [1 + Vlc ■ G^LiE)] (139) 

From Eq. (117) a similar structure as Eq. (134) can be obtained which provides: 

GTl{E) = GTciE) ■ VcL ■ GToiE) (140) 

Furthermore all relations hold for the advanced Green functions in the same way. 
With these ingredients the Keldysh relation(118) can be rewritten as 

G<ciE) =GTciE) • ^^ciE) • G^^SiE) + GTLiE) • 5]<Ji?) • G^SiE) 

=GTciE) • i^^ciE) + VcL • G<o(^) • Vlc] • G^^^iE) 

where G^q(-E) defined by 

GUE) = GToiE) ■ ^<TE) ■ Gf^iE) (142) 

is not exactly the lesser Green function of the pure lead as 'S^TE) depends on the full Green 
function G'"*^*(£'). E.g., this reflects the fact, that the presence of a current through the central 
region will in principle effect the electron distribution in the lead. Nevertheless this reaction 
is typically negligible. Again the complications vanish under the assumption, that scattering is 
neglected in the leads. In the same way one obtains 

G<^iE) =GTciE) • ^^ciE) • G^^UE) + GTLiE) • ^UE) • Gf/(^) 
=GTciE) ■ ^^ciE) ■ G^^SiE) . VcL ■ Gf^iE) 



adv,.. (141) 



+ GTciE) ■ VcL ■ GToiE) ■ ^<TE) ■ Gf^iE) [l + Vlc ■ G^^E) 
--G<ciE) ■ VcL ■ Gf^iE) + GTciE) ■ Vcl ■ G^E) 
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(143) 



where Eqs (139,140,141) have been subsequently applied (partially in the form for advanced 
functions). 

Let us consider a typical structure where the central region is connected to several independent 
leads i, which are translational invariant in their current direction. It is assumed that each lead 
i is disorder-free so that the eigenstates can be separated into transverse and longitudinal parts, 
0A(j(^) = Xex{i^)^l{z), where z is the spatial coordinate in the direction towards the central 
structure and r is a two-dimensional vector perpendicular to z. (A different coordinate system 
is applied for each lead.) The index A numbers the transverse modes within a given lead and 
q denotes the behavior far away from the central region where ^i{z) ~ e'"?^ is assumed. The 
corresponding matrices G^o are diagonal with matrix elements Gi\q{E) in this basis. 

The electric current from lead i and mode A into the central region can be obtained from 
Eqs. (102,119). 

2 r dE 

h. ^ 



lex 



2 (for spin)e 



■ ^ ^ Re j UiXq,aG'^i>^q 



{E)} 



(144) 



where then sum ^^ runs over all states belonging to the central region. Eq. (143) provides: 



4e 



Ux 






■ zZ X/ ^l3/XqUaq,oL \G^p{E)G'^q{E) - G'^p{E)Gf^q{E) 



h / 27r 



Q «/3 



+ G%{E)G<^^{E) - Gf^{E)G<^^{E) 



=T I ^^' {^cciiX, E) . [iG<c{E) + hx{E)i (g^SUE) - G^^^iE))]} 

(145) 

In the second line, the real part was taken by adding the complex conjugated expression (with 
a and /? exchanged). In the last line the definition 



r^„(^A, E) = i2j t^/3/Agf^ttQ,Q Gfxg{E) - G'ifgiE) 



(146) 



is introduced and it is assumed that the occupation of the modes A in the lead i can be treated 
by a distribution function fix{E) with 



GfxqiE) = hx{E) G^liq, E) - G\%{q, E) 



(147) 



Equation (145) has been derived in [146,38,39]. Together with Eqs. (138,141) one obtains a closed 
set of equations for the Green functions in the central regions which is a convenient starting point 
for the simulation of quantum devices, see. Equation (145) can also be used as a starting point 
to derive the Landauer-Biittiker formalism [147] which can be easily applied if S = inside the 
central region [148]. Otherwise some complications arise as discussed in [145]. A generalization 
to take into account time-dependent phenomena is straightforward [149]. 
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4.2. Application to the superlattice structure 

While the discussion of nonequihbrium Green functions was quite general in the preceding 
subsection, the general formalism will now be applied to a superlattice structure, which is as- 
sumed to be infinitely long. Therefore we restrict ourselves to the central region here. We use the 
basis given by the products of Wannier states multiplied by plane waves in the direction parallel 
to the layers, ^^ni^) ^ "" / v^- Then the general states are given by |a) = \n, u, k). The Hamilto- 
nian Hq is given by Eq. (23) and U is given by Eqs. (24,25). For simplicity we restrict ourselves 
to the lowest level ^ = a (and omit the respective indices) and nearest neighbor coupling Ti in 
the following. Furthermore we set E"" = 0. (The inclusion of higher levels is straightforward but 
tedious.) We assume that the superlattice is (after impurity averaging) spatially homogeneous 
in the x,y plane. Then the expectation values ([am,(ki)]^a„(k2)) must vanish for ki ^ k2. Thus 
the Green functions are diagonal in the wave vector k and can be written as G„^m(k; ti, ^2) in 
the following. 

4.2.1. Basic equations 

The Wannier functions are essentially localized within a single quantum well. Therefore the 
scattering matrix elements connecting states of different quantum wells are small compared to 
those describing intrawell scattering. Thus, we restrict ourselves to scattering matrix elements 
Vmk,nk' {{^i}) 1 ^mk.rak' (P) which are diagonal in the well indices m,n in the following (no 
interwell scattering). 

In this case the scattering-induced currents between different wells vanish in Eq. (102) and 
the total electric current density well m to well m + 1 is given by 

J,„_„,, , 2(fo,-Sp.n)e ^2 f ^^ ^T^a<,,„JK E,} (148) 

where Eq. (119) has been applied. Similarly the electron density in well m is given by 

«,» = ^^S^E/fl-{Gl.(><,iJ)} (149) 

k -^ 

We consider impurity and phonon scattering within the self-consistent Born approximation. 
Furthermore we neglect correlations between the scattering matrix elements in different wells. 
This means that (T^mk,mk'({^i})^nk',nk({^j}))imp vanishes in Eq. (121) for n / m. This assump- 
tion is realistic for short range potentials of random impurities but becomes problematic if 
interface roughness scattering is considered where significant correlations between neighboring 
wells may occur [150,151]. For phonon scattering this approximation makes sense if localized 
phonons such as interface phonons are considered which are different in each well. Under this 
assumption the self-energies become diagonal in the well index and can be written as 

S</-t''-P(k, E) = j;(Kk,nk'({^l})Kk',nk({r-;}))impG'</:^*(k', E) (150) 

k' 

and in the same way for phonon scattering with phonon modes l,pwe find from Eqs. (127,128): 
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5.r.t,phonon(k,^) = ^ | M„k,nk' (P, P 



p,l,k' 

fret /I / 



[nB(M(p1) + 1] G'^'i{k',E- nLOem 



1 1 (151) 

+-G< (k', E - fuviip)) - -G<{k', E + fuviip}) 



2 "•"' ' --VX-// 2 






and 



S<'Pi>°°°°(k,i?) = Y, \Mr,^,n^.'{p,l)\'[nB{fu^i{p))G<„{-k',E - fuviip)) 

p,i,k' (152) 

+ [fiBihwiip)) + 1] G<„(k', E + M(p)] 

Then Eq. (116) becomes 

[E - Ek + eFdmi - S^*(k, S)) G'^^^^,^{\s.,E) - TiG'^^_^„^^{\<i, E) - TiG''^^^i^^^{\<l, E) 

and the Keldysh relation becomes 

G^^^^^{k,E) =Y,GZ,n{^^E)T.<{k,E)Gt;^^%E) . (154) 

n 

Together with the functionals for the self-energies, Eqs. (153,154) form a closed set of equations 
which will be solved in subsection 4.4, where explicit results are presented. 
In the same way the difference between Eq. (114) and Eq. (115) gives 

(m - n)eFdG<,„(k, E) =n [G<_i,„(k, E) + G<+i,„(k, E) - G<,„„i(k, E) - G<,„+i(k, E)] 

+ E-*(k, E)G<^^{k, E) + S< (k, E)G^^{k, E) 

- G-t„(k, i?)S<(k, E) - G<,„(k, E)i:f^{k, E) 

(155) 

which will be the starting point for the derivation of the miniband conduction and Wannier-Stark 
hopping model in appendix B. 

4.2.2. Constant scattering matrix elements 

For numerical calculations the k-dependence of the self-energy is difficult to handle, as a two 
dimensional array of Tiik^E) has to be evaluated and stored in each calculation cycle. Thus 
the problem becomes much simpler if constant (i.e. momentum independent) scattering matrix 
elements are assumed. Then the self-energy for impurity scattering can be written as 

s</-*(k,£;)=^(K„({f,})K„({f,}))i^pG</r(k',i?) 

'^' ^ /.oo (156) 

=A{Vnn{{r,})Vnnm})U/-^ j^ di?fc,G</r (k', E) 
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which does not depend on k. (Note, that Gn,n only depends on |k| = yj2Ekm/fi due to the 
rotational symmetry in the (x,y)-plane.) As G5^*'^(k', i?) ~ —\/Eki for large E^i, the respective 
integral exhibits a logarithmic divergence. This can be either cured by applying a finite cut-off 
^max Qj, i^y adding \/{Ej.i + T) in the integrand with an arbitrary but constant value F > 0. 
Throughout this work the total scattering rate of free-particle states is used here. This procedure 
adds a constant real term to the retarded and advanced self-energy, which effectively renormalizes 
the energy scale to E + \A{Vnn{{rj})Vnn{{rj}))impPo log(£'™'^^/r + 1) but does not change the 
physics, which only depends on energy differences. 

For free-particle Green functions G^^\^i\i.\ E) = 1/{E — E^' + neFd + iO"'") the integral yields 



Po 



J:T{k,E) = ^(Kn({r-,})Kn({rS}))imp^ 

implying a scattering rate (for E > —neFd) 



log 



E + neFd 



m@{E + neFd) 



1 ^ --Im|E''=H = ^/^o^(^"»(i^i})^""(i^j}))imp 



(157) 



(158) 



which will be used as the parameter for the impurity scattering strength in the subsequent 
calculations. [^ 

Phonon scattering is treated in the same way with a constant matrix element Mn{p, I) = Mi. 
Here the scattering rate is defined via 



1 



TTpoAM^ 



h 



(159) 

''"phonon, i 

which is the free-particle spontaneous phonon emission rate if the final state is energetically 
available and Pauli blocking is negligible at low densities. 



4.3. Solution for constant self-energy 



The nature of the retarded Green function determined by Eq. (153) can be analyzed by 
an analytical evaluation for the constant self-energy S^*'*(k, £■) = — ir/2. Then one finds with 
Eq. (9.1.27) of [72]: 



G-;„(k,i?)=^ 



'Jm-j \eFdJ "'"-J KeFd) 

E + jeFd -Ek + i^ 



(160) 



which is a superposition of broadened Wannier-Stark states, compare Eq. (31). The time- 
dependent Green function can be evaluated with help of Eq. (9.1.79) of [72]: 



G^*„(k;ti,i2) =-ie(ti -t2)i"-'"e-'(^^-^-™)(*i-*2)Ae-J^(ti-t2) 

^ 'Jm—n 



4Ti . feFd, 

sm I -z^[ti - t2) 



eFd 



2h 



(161) 



Note that the product A{V„n{{rj})Vnn{{rj})}iinp typically does not depend on the sample area A. 
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This formulation is quite helpful in order to study the different limits. At first one notes that 
the off-diagonal elements of G^^^ start to contribute significantly if 



-^T-.sm.\ -TT^ih -t2) 



> 1 (162) 



eFd V 2/i 

Now only the times < (ti — ^2) ^ 2h/T contribute to G^'^^ due to the exponential decay. If 
eFd > r the sine may oscillate over many periods with the average absolute value of 2/tt ~ 1/2. 
Then Eq. (162) gives 2|ri| > eFd. If otherwise eFd < T we may replace sin(x) ~ x and Eq. (162) 
gives 2|Ti| > r at ti — ^2 = ^/r. Thus we conclude that the states are essentially delocalized if 

2|ri|>r and 2|Ti| > eFd (163) 

On the other hand for 

2|Ti|<r or 2\Ti\<.eFd (164) 

^ntn becomes small for m ^ n and the states are essentially localized. Furthermore for 

eFd > r (165) 

the poles at E = —iT/2 + jeFd are clearly resolved in the energy dependence of Eq. (160) which 
indicates the persistence of the Wannier-Stark ladder under scattering. The ranges (163,164,165) 
correspond to the regimes of validity for miniband transport, sequential tunneling and Wannier- 
Stark hopping, respectively, as given in Fig. 7. In Appendix B it will be shown explicitly that 
the respective transport equations can be recovered from the model of nonequilibrium Green 
functions in these ranges. 

4.4. Results 

Now the results from explicit calculations of the full quantum transport model are presented. 
Let us first summarize the underlying assumptions as applied successively in subsection 4.2. 

- Stationary transport in a homogeneous electric field is considered where only the lowest sub- 
band of the superlattice is taken into account. 

- The scattering matrix elements are assumed to be diagonal in the well index (no interwell 
scattering) and no correlations between scattering matrix elements belonging to different wells 
exist. 

- The scattering is treated within the self-consistent Born approximation (thus neglecting effects 
due to weak localization [152,153]) using momentum-independent matrix elements. (The latter 
assumption is equivalent to a localized scattering potential, which may be problematic for 
phonon scattering) 

Within these assumptions it is possible to solve the system of equations Eqs. (150,151,152, 
153,154) self-consistently. Note, that in all these calculations temperature only enters the occu- 
pation of the phonon modes, which are assumed to be in thermal equilibrium. The electronic 
distribution is calculated self-consistently without any assumption of (heated) equilibrium. 

Throughout this subsection we use the scattering time Topt = 0.125 ps and the phonon energy 
^opt = 36 meV for optical phonon scattering, which is the dominating energy relaxation process 
in III-V materials. In order to guaranty energy relaxation for particle energies below 36 meV, 
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Fig. 15. Current- voltage characteristic for the superlattice studied in [154]. Full line: Result from quantum trans- 
port for A'' = 100 wells. Dots: Experimental data (courtesy of E. Schomburg). Dashed line: Theoretical result with 
a serial resistance of 7.2 x lO^^ficmV^ (T = 300 K). 



acoustic phonon scattering has to be taken into account. We mimic these phonons by a second 
phonon with constant energy HuJac = /i'^opt(v5 — 1)/10 ~ 4.4498 meV. Here hujs_c should be 
less than ksT (~ 6.4 meV at 77 K), so that this mechanism can be efficient close to thermal 
equilibrium. Furthermore the ratio Wac/'^opt was chosen irrational in order to avoid spurious 
resonances. The respective scattering time is chosen t^c = 5 ps. 

At first consider the strongly coupled superlattice studied experimentally in [154]. It consists 
of 100 periods with 3.45 nm GaAs wells and 0.96 nm AlAs barriers, yielding a coupling Ti = 
—20.5 meV. The doping density provides Nd = 3.6 x lO^'^/cm^ and an impurity scattering 
time Timp = 0.12 ps is applied, which can be estimated from the impurity scattering rate for 
this doping range [120]. Results for the lattice temperature T = 300 K (which only enters the 
phonon occupation number ub) are shown in Fig. 15 (full line). One encounters the typical 
shape for superlattice transport. While the current peak is in good quantitative agreement with 
the experimental observation the peak position is shifted significantly. The agreement becomes 
excellent if a serial resistance of the order of lOil (for the experimental sample area) is included 
which may result from contacts, leads, or the substrate. 

Several simulations have been performed for a model superlattice with Ti = —5.075 meV, d = 
5.1 nm, and Nd = 5.1 x 10^/cm^, which has been extensively studied by S. Rott [92] (see also 
Figs. 8 and 9 in section 3). The rather low doping is taken into account by an impurity scattering 



time T\ 



imp 



0.333 ps in agreement with semiclassical Monte-Carlo simulations. The results for 
the drift velocity f drift = J /{eNo/d) for T = 77 K and T = 300 K are shown in Fig. 16 (full 
line). For comparison, Monte-Carlo simulations within the miniband transport model (Sec. 3.1) 
have been performed by S. Rott and A. Markus (Institut fiir Technische Physik der Universitat 
Erlangen) applying the same scattering matrix elements. For low and moderate field strengths up 
to eFd ~ 20 meV the relation obtained from miniband transport agrees extremely well with the 
full quantum transport result. At larger fields, eFd exceeds the miniband width and miniband 
transport is no longer valid as discussed in subsection 4.3. In Fig. 17 the electron distribution 
functions (for T = 77 K) are depicted for various electric field strengths. Comparison of the 
dashed and full lines for miniband conduction and quantum transport shows that the distribution 
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Fig. 16. Drift velocity versus field for a superlattice with A = 20.3 meV, d = 5.1 nm, Nd ~ 5.1 x 10^/cm^, and 
Timp = 0.333 ps. Full line: Quantum transport. Dashed line: Miniband transport. (From [155]). 



functions agree very well for eFd < 20 meV in the energy range E^ < 30 meV. For higher values 
of E). the quantum mechanical distribution function is larger than its semiclassical counterpart. 
Here the occupation is quite small, so that the tail from the broadened spectral function ~ '^ /Ef. 
becomes visible [156,157], in good agreement with the estimate of Eq. (133). In compensation 
the quantum result is slightly smaller than its semiclassical counterpart for low values of E]^ , as 
the total density has to be the same. Nevertheless, these effects are small compared to the typical 
occupation numbers and one can conclude that the Boltzmann equation for miniband transport 
gives reliable results both for the current density and the electron distribution. Actually, the 
quality of agreement is stunning regarding the crude assumptions necessary for the derivation 
of the Boltzmann equation in Appendix B.2. 

The velocity-field relations from Fig. 16 can be directly compared with those in Fig. 8, where 
the correct matrix elements have been applied within the semiclassical miniband transport 
model. There are no qualitative differences, so that one may conclude that the assumption 
of constant matrix elements as well as the artificial acoustic phonons do not cause unphysical 
results. 

Let us now study the details in the velocity-field relation and through a glance at the respective 
distribution functions. At very low fields one encounters a linear increase of the drift velocity 
with the electric field. Here the distribution function resembles a thermal distribution with the 
lattice temperature (dotted lines for eFd = 0.03 meV). The slight shift can be easily treated 
within linear response, yielding a field-independent mobility. 

For higher fields {eFd = 0.3 meV) one encounters a sublinear increase of the drift velocity. Here 
significant heating occurs, which can be seen from the respective distribution function in Fig. 17. 
The distribution functions resembles a shifted Fermi distribution with an electron temperature 
Te ~ 140 K for E]^ < /ii^opt • For higher energies optical phonon scattering is still efficient and one 
encounters a steeper decrease of f{Ek) with energy. This shows that the concept of an electron 
temperature makes sense in this range both for quantum and semiclassical transport even without 
electron-electron scattering. For electric fields close to the current maximum [eFd = 2 meV) the 
distribution functions strongly deviates form any thermal or heated distribution function. In the 
semiclassical picture this can be viewed as a result from frequent Bragg-scattering where particle 
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Fig. 17. Distribution functions /(g, k) and /(-Efc) = d/{2n) f dq/(g, k) of the carriers for different field strengths 
for the superlattice studied in Fig. 16 at T = 77 K. Full line: Diagonal elements of the density matrix p{q, k) as 
calculated from the quantum transport model. Dashed line: semiclassical distribution function as calculated from 
the Boltzmann equation for miniband transport. Dotted line: Thermal distribution npiEq + Ek — jJ,) applying an 
electron temperature Te = 77 K for eFd = 0.03 meV and Te = 140 K for eFd = 0.3 meV. 
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Fig. 18. Drift velocity versus field for a superlattice with Ti — 1 meV, d — 5.1 nm, No ~ 5.1 x 10^/cm^, and 
Timp ~ 0.0666 ps. Full line: Quantum transport. Dashed line: Miniband transport. 

reach the states with q = n/d. If the electric field becomes even stronger, the electrons may 
traverse the g-Brillouin zone several times without scattering for eFd ^ h/r in the semiclassical 
miniband picture. This leads to a flat electron distribution in q as can be observed for eFd = 50 
meV in Fig. 17. The same holds for the quantum distribution function although the reason is 
different. In this fleld range (2|Ti| <^ eFd) the Wannier-Stark states are essentially localized to 
a single well. Thus a semiclassical occupation of the Wannier-State creates a flat distribution in 
g'-space as well. 

While both approaches explain the flat distribution in g-space, significant differences in scat- 
tering arise. The electron running through the q states exhibits different energies E{q) during 
passage. This provides different selection rules for scattering then the presence of a Wannier- 
Stark states with fixed energy. Therefore the /(-Efc)-distribution calculated by the quantum 
transport model exhibits pronounced features on the energy scales eFd and ^opt as well as 
the difference eFd — hiOopt because at these energies new scattering channels appear. If these 
scales match, the phonon resonance at eFd = hujopt appears in the velocity-field characteristics, 
see Fig. 16. In contrast, the semiclassical result exhibits a rather fiat f{Ek) distribution and no 
phonon resonance can be observed. 

The situation changes for small coupling and strong scattering, when F > 2|Ti| holds. Fig. 18 
shows the velocity field relations for Ti = — 1 meV and Timp = 0.0666 ps. Here the semiclassical 
miniband transport calculation strongly deviates from the full quantum result both in the low- 
field and in the high-field region. Nevertheless, the agreement gets better for higher temperatures. 



In Fig. 19 the quantum transport calculations are compared with results from the models of 
Wannier-Stark hopping and sequential tunneling. In all calculations identical scattering matrix 
elements are used. The calculations for Wannier-Stark hopping have been performed by S. Rott 
and A. Markus (Institut fiir Technische Physik der Universitat Erlangen) using self-consistent 
distribution functions, see section 3.2 for details. One finds that the Wannier-Stark hopping 
model provides good results (including the phonon resonance) in the high-field region where 
eFd ^ F (F ~ 7 meV and 15 meV for the left and right superlattice, respectively, for energies 
at which phonon emission is possible). The agreement deteriorates significantly if the simple 
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Fig. 19. Drift velocity versus field for a superlattice with d — 5.1 nm, Nd ~ 5.1 x 10^/cm^, Ti = 5.075 meV, 
Timp ~ 0.33 ps (left) as well as Ti = 1 meV, ri^p = 0.0666 ps (right). Full line; Quantum transport. Dashed line: 
sequential tunneling. Dotted line: self-consistent Wannier-Stark hopping. 

version of Wannier-Stark hopping without self-consistency is applied. 

The sequential tunneling model clearly fails at low and moderate field strengths for the strongly 
coupled superlattice (left part of Fig. 19) as 2|ri| > F. In the high field region eFd ^> 2|Ti| the 
results becomes valid and the phonon resonance can be observed within the sequential tunneling 
model. For weakly coupled superlattice (right part of Fig. 19), the sequential tunneling model 
should be applicable and the low-field conductance is in good agreement with the full quantum 
transport result. Significant deviations occur at intermediate fields because electron heating 
is not included in the sequential tunneling model. As the current density drops with electron 
temperature (see section 3.3.2), the drift velocities are too high, if heating is neglected. These 
heating effects have been recently taken into account within the sequential tunneling model 
[158,159]. 

In Fig. 20, a comparison between quantum transport and sequential tunneling is shown for 
the weakly coupled superlattice studied in [123,124], see also section 3.3.3. Excellent agreement 
is found between both approaches and the results agree well with those presented in Fig. 11 
applying realistic impurity scattering matrix elements and neglecting phonon scattering. 

The observations presented here correspond to the boundaries of validity given in Fig. 7, 
which can be considered as a reliable guide. Furthermore the results indicate that the models of 
miniband transport and sequential tunneling give qualitatively reasonable results even outside 
their range of applicability. Thus they can be used as a rough guide to obtain a correct order of 
magnitude of the current. 
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Fig. 20. Current-voltage characteristic for the superlattice studied in [123,124] with parameters T\ = —0.02 meV, 
No = 1.5 X lO'^'^/cm'^, d = 13 nm, rimp = 0.0666 ps, and T = 4 K. Full line: Quantum transport. Dashed line: 
sequential tunneling. 
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Fig. 21. Experimental results for a semiconductor superlattice with 40 periods consisting ol 4 nm AlAs barriers 
and 9 nm GaAs wells. The doping density is Nn = 1.5 x 10^^/cm^ and the sample area is A = 1.13 x 10~*m^. (a) 
Current-voltage characteristic for negative bias applied to the top contact, where a stationary current is observed 
(from [179]). (b) Current-time signal for a constant positive bias exhibiting self-sustained current oscillations (from 

[181]). 

5. Formation of field domains 

As shown in sections 3 and 4, semiconductor superlattices typically exhibit ranges of negative 
differential conductivity. This occurs after the first current peak for miniband transport (or 
sequential tunneling) as well as after the subsequent resonances at higher fields when different 
levels align. The shape of the current-field relation is thus typically of N-type and domain 
formation effects are likely to occur (see [45] for a general overview). The prototype of an 
extended device with N-shaped current-field relation is the Gunn diode which exhibits self- 
sustained current oscillation due to traveling field domains [160-166]. A similar behavior has 
been suggested for semiconductor superlattices [18,167], and oscillatory behavior has indeed 
been found experimentally in the last years [19,168] with frequencies over 100 GHz [20,169]. 
In contrast to the Gunn diode, semiconductor superlattices frequently exhibit the formation 
of stable stationary domains which lead to a characteristic saw-tooth pattern in the current- 
voltage characteristic (see Fig. 21a) as observed by many different groups [16,17,123,124,170- 
174]. (See also [175,176] for domain formation in a parallel magnetic field.) The measurements 
are typically performed by applying a continuous sweep of the bias. As the branches overlap, one 
observes different parts of the branches for sweep up and sweep down of the bia^^^l Sometimes 
a quite complicated behavior is observed as well. In particular for the superlattice structure 
of Fig. 21 stationary field domains [179], self-sustained oscillations [19], as well as bistability 
between stationary and oscillatory behavior [180] has been observed under fixed bias conditions 
within the first plateau of the current-voltage characteristic. 



^° Ofcourse the variation of the bias must be slow, so that the field domains can follow adiabatically. Otherwise, 
the sawtooth structure disappears [177,178]. 
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Fig. 22. Sketch of a superlattice structure with an inhomogeneous field distribution. 

In this section such complex behavior wih be analyzed by a comparison between numerical 
results based on the sequential tunneling model from section 3.3 with analytical studies. It will 
be shown that most of the observed effects can be understood as the result of a competition 
between two mechanisms: (i) The motion of fronts connecting low- and high-field domains and 
(ii) the dynamical evolution of the field close to the injecting contact, i.e., the cathode. 

This section is organized as follows: First a model is described which is able to reproduce most 
of the experimental findings for weakly coupled superlattices. In subsection 5.2 numerical results 
are presented reproducing the behavior shown in Fig. 21. In the subsequent subsections the key 
elements, namely the traveling fronts (subsection 5.3) and the cathode behavior (subsection 5.4) 
are discussed separately. In subsection 5.5 it will be demonstrated how their combination explains 
the behavior observed both in experiment and simulation. Finally, the findings are summarized 
and an instruction for the analysis of nonlinear superlattice transport behavior is presented. 

5.1. The model 

In weakly-coupled multiple quantum wells the electronic states are essentially localized in 
single wells forming energy levels E'^. Transport then occurs by sequential tunneling between 
neighboring wells. The current from well m to well m + 1 is modeled by a function J„^(„4.x) = 
J{Fm,nm, nm+i), where F^ is the average field drop between the respective wells and Um denotes 
the electron density (per unit area) in well m as depicted in Fig. 22. Considering a superlattice 
with N wells embedded between N + 1 barriers, the dynamics is determined by the continuity 
equation 



drirr, 
"dT 



Ji 



(m-1)- 



Ori 



*{m+l) 



for m = 1, . . . N . 



The electric field satisfies Poisson's equation 



ereoiFm - Fm~i) = e{nm - Nd) ior m = 1, . . . N 



(166) 



(167) 



where Nd is the doping density per period (per unit area), and e^ and eg are the relative and 
absolute permittivities. Finally, the total voltage U is determined by 
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N 



U{t) = Y, Fmd (168) 



m=0 



where we have neglected a possible voltage drop at the contacts for simplicity. Equations (166- 
168) can be transformed into an equivalent set of equations: 



ereo-T^ = (^(i) - ^m^(m+i)) for m = 0, . . . iV (169) 

n^ = ND + ^-^{Fm-Fra-i) for m = 1, . . . iV (170) 

e 

{N + i)j{t) = x; J™_(„+i) + ^^ (171) 

i=o 

This shows that the dynamical evolution of the local fields is driven by the total current density 
J(t), which itself is determined by the global behavior of the sample. Thus, J{t) represents a 
global coupling. In order to obtain a closed set of equations, the currents across the first (cathode) 
and the last (anode) barrier, Jq^i and Jjv^(Ar+i), respectively, have to be specified. Here the 
following approaches have been taken previously: 

In [66,182-184] these contact currents were calculated within the assumption of two fictitious 
additional wells, one before the first and one after the last barrier. Then the current Jo— >i is 
given by J(Fo,no,ni) for tunneling between two wells, where the fictitious density uq has to be 
specified, usually assuming no = (1 + c)N£), with c > —1. The current across the last barrier is 
treated analogously by introducing a fictitious density un+i- This approach will be referred to 
as constant density boundary condition in the following. 

Alternatively, one may assume that the current is proportional to the local field, i.e., AJq^i = 
aF^d with an Ohmic conductance a. For the anode condition J^v^rAr+i) one has to take into 
account, that the current must vanish if un tends to zero, as otherwise the density njv can become 
negative. This can be ensured by an additional factor uj^/Nj:, and ^Jjv— >(Ar+i) = (jFf^dnN/N£) 
is a reasonable choice. This will be referred to as Ohmic boundary condition. 

The actual potential distribution at the boundary could be taken into account within a 
transmission-type formalism [67]. Here I will restrict myself to the two approaches sketched 
above, which are simpler and provide an understanding of most of the underlying physics. They 
contain the main ingredients to understand more complicated and physically better motivated 
versions. 

5.2. Numerical results 

Let us consider the GaAs/AlAs superlattice structure from Fig. 21 (which exhibits a rather 
small miniband width of 0.08 meV) as a model system for the subsequent calculations. In the 
following the function J{Fm, nm, n„i+i) is calculated from the nominal sample parameters by the 
model described in section 3.3. (See also Fig. 11.) The result for a homogeneous electron density 
Um = nm+1 = Nd is shown in Fig. 23. Here, the current / = sgn(e)^J is shown applying the 
experimental sample area A = 1.13 x 10~®m^ in order to facilitate comparison with experimental 
data. In this section the factor sgn(e) is included, so that the sign of the current equals the 
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Fig. 23. Current-field relation I{eFd,nm = ND,nm+i ~ Nd) calculated from the microscopic model for the 
superlattice of Fig. 21. The relation for the cathode current /o^i = sgn(e)aFd is also shown for two different 



values of a (broken lines) . The hatched area indicates the current range F, 
are found. 
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particle current in transport direction which simplifies the following discussion. Fig. 23 shows 
that the current-field characteristics exhibits a linear conductivity for low fields, a maximum at 
eFd = eFmax^, a range with negative differential conductivity for eFd > eF^naxd and a second 
sharp rise of the current for higher fields eFd > eF^-amd due to resonant tunneling from the 
lowest level to the second level of adjacent wells. 

Different approaches to obtain the local current density J{Fm,nm,nm+i) can be applied as 
well and the overall results do not depend on this choice, provided the general shape resembles 
Fig. 23. E.g., in Refs. [19,66,67,182,174,185] one can find similar results like those presented in 
this section, where varios types of local current density functions are applied. 

By simulating equations (166-168) for a fixed bias U until a stationary state is reached the 
stationary current- voltage characteristic shown in Fig. 24 are obtained. Here the initial condition 
for the calculation is taken from the result of the previous voltage point, simulating a sweep-up 
of the bias U. In Fig. 24a we use the constant density boundary condition with c = 0.5, while the 
Ohmic boundary condition with a = 5 mA/V was applied in Fig. 24b. Both characteristics are 
almost identical and in reasonable quantitative agreement with the experimental data presented 
in Fig. 21a. As shown in the insets, the branches are due to the formation of electric field domains 
inside the sample, where a low-field domain is located at the cathode and a high-field domain is 
located at the anode. Close to the contacts we find a small transition layer, which depends on the 
contact boundary conditions. The domain branches span a fixed current range /^°™ < I < I^"^ 
with /^<?™ = 24.7;uA and I^^ = 85.5^A, indicated by the hatched area in Fig. 23. 

For different boundary conditions oscillatory behavior can be obtained for the same sample 
parameters. This is shown in Fig. 25 for the Ohmic boundary condition with a = 0.5 niA/V. The 
current signal resembles the measured signal displayed in Fig. 21b. Nevertheless different current 
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Fig. 24. Current-voltage characteristics exhibiting stationary field domains as shown in the insets for differ- 
ent boundary conditions: (a) /o^i = /(Fo, 1.5A^D,ni) and /jv^jv+i = /(-Fo,n-iv, I-SA'^d), (b) Jo^i = aFod and 
In^n+1 ~ (TFNd{nN /Nd) with o = ^mA/V . (Parameters as in Fig. 21). 

signals have been observed as well, both experimentally and in numerical simulations (see, e.g., 
Fig. 10.12 in [35] for results obtained from the same model with different boundary conditions). 
The oscillations reported here are due to traveling high-field domains which is different from 
the calculations presented in [181,186,187] where traveling monopoles have been reported. This 
is due to the difference in the boundary conditions and in the transport model used in both 
calculations, see also Refs. [185,188], where this point is analyzed. 



5.3. Traveling fronts 

In this subsection traveling fronts will be examined which will form one of the building blocks 
to understand the global dynamics. Eq. (169) shows that the dynamical evolution of the electric 
field is determined by the total current I{t). The corresponding dynamics can be studied most 
easily for a constant total current corresponding to current controlled conditions. For Iram = 
7.1/uA < I < /max = 326/LtA we have three intersections of / with the homogeneous current-field 
relation shown in Fig. 23. They correspond to three stationary homogeneous field distributions 
eF < eF < eF . Linearization of Eq. (169) shows that the field F is unstable under 
current-controlled conditions as dJ/ dF < 0. Thus a homogeneous initial field distribution will 
either tend to F or to F in its temporal evolution. 

For an appropriate inhomogeneous initial field distribution a part develops to F while another 
part reaches F . Thus a front between these two spatial regimes appears. Calculations on the 
basis of Eqs. (169) - (170) show that the front develops on a typical time scale of less than 
1 //s and afterwards travels through the sample with unaltered shape, as shown in Fig. 26. Here 
two types of fronts appear. Accumulation fronts connect a low-field region on the left side to a 
high-field region on the right side as shown in Figs. 26a, 26c. Eq. (170) shows that the carriers 
accumulate {um > No) in the transition region. If, on the other hand, the high-field region is 
located on the left side, see Figs. 26b, 26d, a depletion front is present. 

As can be seen from Fig. 26 the front velocities depend on the external current and differ for 
accumulation fronts with velocity Cacc and depletion fronts with velocity Cdcp- These velocities 
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Fig. 25. Self-sustained current oscillations for fixed bias conditions sgn(e)L'^ — 1.2 V and an Ohmic boundary 
condition a = 0.5mA/V. (a) Current I versus time t. (b) Density plot of the electric field distribution. Tfie fiigh 
field region is black, (c) Electric field profile at different times. (Parameters as Fig. 21). 

have been determined from a series of simulations as a function of / and are given in Fig. 27. 
Note that Cacc(-f) becomes zero in a finite range of currents 24.7/iyl < I < 85.5 fiA corresponding 
to the range /^°Jf < I < /ma™ of the stationary domains discussed above. (Such a stationary 
front is shown in the inset of Fig. 24.) For higher currents Cacc(/) is negative, i.e., the front 
travels upstream against the direction of the average drift velocity of the electrons as shown in 
Fig. 26c. In contrast Cdep(/) is positive for all currents for the sample parameters used. These 
functions Cacc (/) ) Cdep (/) have been shown to be very helpful to understand and analyze Gunn 
oscillations [166,189]. Here they are applied in the context of semiconductor superlattices where 
some peculiarities can be found. In the following two subsections the special shape of Cacc(/) and 
Cdcp(/) will be explained. 



5.3.1. Continuum limit and fronts traveling backwards 

In order to understand the simulations, let us first discuss the continuum limit of vanishing 
superlattice period. We can approximate 



J, 



m— >(m+l) 



e — —V[rrr 



eD{F„ 



77.772+1 ^ri 



d d'^ 

with the average drift velocity v{F) and an effective diffusion constant D{F) given by 



(172) 
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Fig. 26. Temporal evolution of different initial field profiles under constant current conditions. (Parameters as in 
Fig. 21). 

d^dJ{F,ND,n2) 



v{F) 
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■J{F,Nd,Nd) and D{F) 



(173) 



This approximation becomes exact in the nondegenerate limit, see Eq. (88), which is not appli- 
cable for the sample discussed here. Replacing finite diff'erences by derivatives, Eqs. (169,170) 
are transformed to 

dF{z,t) J 



-^-„(F)-„(F)5£|ri^ + o(F)«'^<^-''> 



(174) 



dt e^eo derSo ^ ^ ^ ^ dz ' ^ ^ q^2 

This is the standard equation for the Gunn efi'ect, for which dcacc(^)/ d/ < and dcdcp(^)/ d/ > 
holds [163]. Neglecting diffusion, the fronts can only travel in the direction of the drift velocity, 
i.e. c(/) > 0, which can be shown by the method of characteristics [190]. We conclude 

without diffusion term: cjep > and Cacc > for /min < I < Imax 

In contrast, if the drift term v{F) dF/dz is neglected, one obtains a nonlinear reaction-diffusion 
system. Here solutions of traveling accumulation fronts F{z, t) = f{z — Cacct) exist with a mono- 
tonic function f(z) and /(— cxd) = F^ and /(oo) = F^^^ . Furthermore one finds Cacc > if 
Aa.rea.iJ) < where 
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Fig. 27. Front velocities for accumulation fronts (cacc, full line) and depletion fronts (cdep, dashed lines). The 
velocities arc approximately linear functions of the current outside the range shown here, /i and I2 denote the 
current where Cacc = Cdcp and 2cacc = Cdep, respectively. (Parameters as in Fig. 21). 

see section 2 of [191]. The depletion fronts can be obtained applying the symmetry operation 
z -^ —z, thus Cdcp(^) = — Cacc(<^)- Typically, the integral Aarea is positive for I^o < I < /max and 
negative for I^^^ < I < ho, where Ico is the current, where both field domains coexist. (This 
value satisfies the equal area rule [192] ^arca(-^co) = 0.) Then one finds the generic scenario: 



without drift term: 



Cdcp < and Cacc > for I„iin < I < Ico 

Cdcp = and Cacc = for / = Ico 

Cdcp > and Cacc < for Ico < I < /max 



Including the drift term, it is obvious that both velocities will increase and the intersection 
point Cdcp = Cacc is shifted to a positive velocity in agreement with Fig. 27. If the diffusion 
constant is not too small, a range of negative velocity is likely to remain. This can be estimated 
in the following way: Let Fq satisfy J = eN£,v{FQ). Linearization of Eq. (174) for F{z,t) = 
Fq + F{k) exp[i(/cz — ujt)] gives the dispersion 

u; = -i7 + vok - '\Dok^ (176) 

with vq = v{Fq)^ Dq = D{Fq), and the dielectric relaxation rate 7 = eN^v' (Fq) / {der€o) ■ Thus 
an initial condition F{z, 0) = Fq + e6{z — zq) has the solution 



Fiz,t) 



y/47rDot 



exp 



-It 



{z - Zq- Vptf 

ADot 



(177) 



which essentially travels in the direction of the drift velocity. Nevertheless F{z, t) also grows 
exponentially in the range z < zq provided —7 > v^/iiDQ. This can be interpreted as the 
occurrence of an absolute instability [193] for 
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when a local perturbation from a homogeneous field spreads in both directions. If D is larger 
than the bound given by (178), spatial variations not necessarily travel through the structure and 
either stationary domain structures or backward traveling field domains are likely to occur. This 
condition has been applied to superlattice transport in [194] and seems to be in good agreement 
with data from a variety of strongly coupled superlattices [195]. Nevertheless the condition for 
a stationary front (c = 0) is only valid for a specific current, as the functions c(/) are strictly 
monotonic for the continuous drift diffusion model (174). This is in contrast to the results shown 
in Fig. 27 where Cacc(-f) is zero in a finite interval. Thus we can draw the following conclusions 
from the continuum limit: 

- Backward traveling fronts can occur for large values of DN^i both for accumulation and 
depletion fronts. 

- The front velocities c(/) are decreasing with / for accumulation fronts and increasing for 
depletion fronts. 

Finally note that the simple model for the current used in [181,182,184] (where the Um+i- 
dependence of J{F,nm,nm+i) is neglected) implies D = in Eq. (172). These models are 
appropriate for the second plateau but differences occur for the first plateau discussed here. In 
particular these models cannot reveal fronts traveling backwards. 

5.3.2. Discreteness of the superlattice and stationary fronts 

While the considerations of the previous subsection apply to a continuous system, weakly- 
coupled superlattices form a system where the discretization due to the finite superlattice period 
is essential. As explained in [35,184,196] this leads to stationary domain states, provided the 
transition region between the two domains becomes of the order of the superlattice period. In 
this case the accumulation or depletion front gets trapped within one well. As this pinning can 
occur within a certain range of currents /£!?™ < I < I^^ one observes extended branches in the 

<~> niiii iiidx 

current-voltage characteristic (Figs. 21,24). As shown in [35,184] the sufficient condition for a 
stationary accumulation front reads 

Nd > A^r ^ ""T — (-^min - Fmax) • (179) 

(Similar results have been given in [182,183,196] as well.) For the superlattice structure con- 
sidered, one obtains N'^'^ = 1.2 x 10^'^/cm^, which is smaller than the actual doping density 
Nf) = 1.5 X 10^^/cm^. Therefore stationary accumulation fronts can exist in a certain range of 
currents for this sample. Let us remark that the condition (179) strongly resembles (178) if the 
diffusion constant D{F) = v{F)d/2 for shot noise [194] is used. Again the discreteness of the 
structure is responsible for the occurrence of stable domains. 

For more strongly doped superlattices the velocity of the depletion fronts may become zero 
as well. The corresponding condition is given by [35] 

^max ^inin ^ 

For the superlattice structure considered this relation give N^^ = 4.4 x 10^^ /cm^, which is larger 
than the actual doping. Therefore no stationary depletion fronts occur in this sample. Together 
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with the monotony arguments depicted in the previous subsection these estimations explain the 
shape of the c(/) functions shown in Fig. 27. Rigorous proofs of some of the features discussed 
here are given in [197]. 

In general one can distinguish three types of superlattice: For No < Nff^ there are no station- 
ary fronts. For Nff^ < Nd < Nj^"^ accumulation fronts are stationary within a certain current 

range. Finally for highly doped superlattices with Nd > ^^j^"^ both accumulation fronts and 
depletion fronts can be stationary, i.e., pinned at a certain well. 

5.4. The injecting contact 

Now the influence of the contact boundary condition at the cathode Jo^i(-Fo) shall be inves- 
tigated which is essential for the dynamical behavior. For a given current J, the evolution of the 
field at the cathode, Fq, is given by: 

d P 
ereo^ = J - Jo^iiFo) (181) 

Let Fc{J) be the solution of Jo^i(-F) = J, which forms an attracting point of Eq. (181), i.e., 
dJo^i/dF > 0. Then FQ{t) will tend towards Fc{J). Provided the relaxation time at the contact 
is much smaller than the corresponding time scale on which J changes, one may use the boundary 
condition Fq = Fc{J{t)) to describe superlattice dynamics. If Fc{J) is close to the value of 
F^{J) or F^^\j), a low-field domain or a high-field domain will be injected into the sample, 
respectively [190]. For a pure drift system the condition for the injection of a low-field domain 
is given by eFc{J) < eF^^{J) [45]. This gives a qualitative bound but this does not hold strictly 
in the superlattice system due to discreteness and diffusion. For the Ohmic boundary condition 
depicted in Fig. 23 we thus find that a low- field domain forms close to the injecting contact for 
/ < /c and a high- field domain forms there for I > Ic, where Ic is the current at the intersection 
of /o^i(-Fo) with the homogeneous current-field characteristic. 

In contrast to the cathode where electrons are injected into the sample, the anode contact 
conditions do not play a major role. A boundary layer exists there, which is typical stable (see 
e.g. [190] for a discussion within a drift model). This can be understood from the fact, that 
the perturbations mainly travel through the sample in the direction of the current flow, see 
Eq. (177), even if there might be some response in the opposite direction as well. 

5.5. Global behavior 

The behavior found numerically in Sec. 5.2 will now be explained within the interplay between 
the dynamics of fronts and contacts. Let us restrict ourselves to rather long superlattices where 
the N^£,L criterion 

NdN > iVB^ = 2.09^^^^^ (182) 

is satisfied so that the homogeneous field distribution is unstable [161,162]. In this case either 
stationary domain states, self-sustained periodic current oscillations, or aperiodic behavior occur. 
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Calculations yield N"^ ~ 10^^ /cm^ for the superlattice structure under consideration, which is 
much smaller than 40 A'^^. This criterion has been successfully applied to a variety of superlattice 
structures [195]. 

5.5.1. Formation of stable stationary field domains 

For the sample considered, accumulation fronts become stationary in the current range /^°™ < 
/ < /max • Obviously, these domain states are only possible if the low-field domain is maintained 
at the cathode, i.e. eFc{I) < eF^\l), as discussed in the previous subsection. For o" = 5 mA/V, 
we find Ic = 110/iA > 1'^°^- Thus these domain states are allowed in agreement with the numer- 
ical findings from Fig. 24. The same behavior can be found for the constant density boundary 
condition with c > 0. In this case the field Fq remains in the low field region as long as / < 
(1 + c)/max, which is hardly exceeded unless the superlattice is operated in the second tunneling 
resonance. Therefore similar behavior is observed for both boundary conditions. 

For the superlattice structure considered here, stationary depletion fronts do not occur. There- 
fore one expects that the high-field domain is always located at the anode for stationary behavior. 
This agrees with experimental findings for this superlattice [198]. For superlattices with a higher 
doping the situation is different and field distributions have been observed where the high-field 
domain is located at the cathode [172]. These experimental findings are in good agreement 
with the criterion (180). Simulations [122] give such a field distribution for appropriate contact 
conditions. A further example is shown in Figs. 28c, d. 

These stationary fronts appear within a finite current range I^^ < I < /m'ax- If the maximum 
accumulation occurs in well w-dom) the total bias for accumulation fronts is given by 

U{I, mdom) = mAornF'{I)d + (iV + 1 - mdom)i^'''(/)d + Uc{I) + ^7front(/) ■ (183) 

Here Uc and f/front are correction terms due to the inhomogeneous field profiles at the contacts 
and in the front region, which do not depend on the front position rTT-dorm provided this position 
is not to close to either of the contacts. Due to the periodicity of the superlattice, a shift 
of the front by one period does not change the current. Thus, one obtains up to A^ branches 
U{I^ m-dom) of Eq. (183) with mdom = 1, . . . A^. This results in the characteristic sawtooth pattern 
in the current-voltage characteristics for Figs. 21a, 24. In addition, two field distributions exist, 
which exhibit a homogeneous (except for boundary effects) field distribution F^ = F {I) and 
-Fm = F {I) for all m. The actual number of branches can be smaller due to the inhomogeneous 
field distributions at the boundaries as shown in the insets in Fig. 24. The voltages corresponding 
to neighboring branches differ by [F {I) — F {I)]d. As this difference depends on the current, 
these branches are not exactly reproduced for different mdom) but change their shape. Neglecting 
boundary effects, their slopes are given by 

dU ,dF^(/) ,,, ,dF^^^(I) , ^, 

— = mdomc? ^^ + (iV - mdom + l)t? ^j^ (184) 

which varies between L (1F^{I)/ d/ for small biases (accumulation front close to the right contact) 
and L dF^^^{I)/ dl for large bias (accumulation front close to the left contact). The behavior 
for depletion fronts is identical except that they are located close to the left/right contact for 
small/large bias and the fields F^^^ and F^ have to be exchanged in Eq. (183). 
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5.5.2. Self-sustained current oscillations 

For a = 0.5 mA/V we find Ic = 22fiA, which is sinaher than /^°™. Therefore no stationary 
domain states are possible. In contrast, self-sustained oscillations occur, as shown in Fig. 25. The 
mechanism of these oscillations can be understood by using ideas reminiscent of the Gunn effect 
asymptotic [166]: for t = 0.3)US the total bias is distributed between a high-field domain and a 
low-field domain. As the total bias is constant, there is a certain part of each oscillation period 
during which both boundaries must travel at the same velocity [166,189]. Then Cacc(-^i) = Cdcp(A) 
giving Ii = 0.0114 mA (see Fig. 26). As /i < Ic the cathode remains in the the low-field domain. 
We observe a current signal which is constant in average and exhibits fast oscillations with the 
period c/d due to well-to-well hopping of the accumulation front. This feature has been discussed 
in the analysis of switching behavior [199-201] and is explained in detail in [185]. 

After the leading edge of the high-field domain has reached the anode, the size of the high-field 
domain shrinks as the trailing edge travels further with unaltered velocity. In order to maintain 
the total bias, the fields increase in both domains, leading to an increase of I{t) via the global 
condition (171), as can be seen at t = 2.8/US. When I{t) becomes larger than Ic the field at the 
cathode injects a high-field domain into the superlattice (t = 4.1^s). As Cacc is quite large in 
this range of currents the newly formed domain expands relatively fast and the fields drop in 
order to maintain the bias. Therefore the current shrinks again below Ic and the cathode injects 
a low-field domain into the superlattice (t = 5.6;^s). 

In this situation three boundaries are present in the superlattice. The old accumulation front 
(around well 35) and the depletion and accumulation front limiting the newly formed high-field 
domain (from well 5 to 10). In this situation the sum of the extensions of both high- field regions 
is kept constant if the depletion front travels with twice the velocity of the accumulation fronts, 
i.e. 2cacc{h) = Cdcp(-^2) yielding I2 = 0.015 mA (see Fig. 26). This is just the average current 
observed in the simulation observed in the plateau around t = 5.6/Us. The oscillatory part of 
the current exhibits a lower fundamental frequency compared to the frequency observed around 
t = 0.3fis as Cacc(^2) < Cacc(-^i)- In addition one observes two peaks per period resulting from 
the presence of two accumulation fronts. This behavior is maintained until t = T.l^s when the 
old accumulation front reaches the anode and only the newly formed high-field domain remains. 
Afterwards the cycle is repeated. 

Different scenarios for oscillatory behavior are also possible. For the constant density boundary 
condition with c < oscillations due to traveling depletion fronts have been found for the same 
model, see Refs. [181,35] and Fig. 29d for a further example. In this case the field at the cathode 
remains above eFmin, once this range of fields has been reached. For elJ < NeF^-am the high-field 
domain cannot extend over the whole superlattice and thus a range of fields with eFd < eFmind 
must be present. This requires the existence of a depletion front, which cannot be stationary 
for the superlattice considered here as Nu < Nj^^. Thus the front travels towards the anode 
increasing the width of the high-field domain. In order to keep the bias constant, the field inside 
the domain shrinks during this process. If F becomes lower than Fmin the domain becomes 
unstable and the field tends to the low field value F . As the field Fq remains high due to the 
boundary condition a new depletion front is formed and the cycle repeats itself. An equivalent 
oscillation mechanism is discussed in detail in [181,190] for accumulation fronts. (This mechanism 
is also active in the oscillation shown in Fig. 30b.) It becomes relevant for lower doped samples, 
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Fig. 28. Results for a superlattice with 5 nm Alo.3Gao.7As barriers and 8 nm GaAs quantum wells with a high 
doping density A^d = 5 x 10"/cm^. (a) Current-field relation (full hne) evaluated by Eq. (83) with T" = T* = 8 
meV. The dashed lines depict the cathode current densities used in parts (b), (c), and (d), respectively, (b) Cur- 
rent-bias relation for an Ohmic boundary condition with large a. (c) Current-bias relation for an Ohmic boundary 
condition with small a. (d) Current-bias relation for a constant-density boundary condition with c = —0.2. The 
insets display examples of the respective field distributions. 

i.e., for Nd < Nff^, when these fronts cannot become stationary. 
5.6. Summary 



Let us summarize the findings from the previous subsections in order to obtain a general 

outline for the analysis of semiconductor superlattices. To visualize the general trends, various 

types of behavior are shown for a superlattice test structure in Figs. 28-30 for different doping 

densities. 

(i) Determine the local current-field relation J{Fm,nm,nm+i) from some transport model, 

see, e.g., section 3. The corresponding results for rim = "n-m+i = ^D are displayed in 

Figs. 28a,29a,30a. 

dcp 



rdcp 



ii) Determine Nff^ and A'^^ ^ from Eqs. (179,180). According to their definition Nff^ < N'^ 
holds. From Fig. 29a we find N^^'' « 1.7 x lO^O/cm^ and A^''^ Prf 3 x lO^/cm^. 
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Fig. 29. Same as Fig. 28 for a medium doping density Nd = 10^^ /cm^. (c) and (d) show current oscillations at a 
fixed bias. 



(iii) Compare theses quantities with the actual doping density N^, 

Nj^^ < Nf): Stable domains form both for good and bad contacts at the cathode. The 
domain boundaries can be formed by accumulation layers, see Fig. 28b, or depletions 



layers, see Figs. 28c, d. Both may even coexist if the respective ranges [/j 



dom rdom 
min ' max 



overlap. 



In all cases the current-voltage characteristic exhibits the typical saw-tooth behavior, 
where the number of jumps roughly equals the number of periods. 

jyacc ^ jY^ ^ -^^T^' '^^^ behavior depends crucially on the boundary condition at the in- 
jecting contact: If the current can be injected at a fairly low electric field in the cathode 
(good contact), stable domains are found which exhibit accumulation layers, i.e., the 
high-field region is located at the anode, see Fig. 29b. Otherwise, for bad contacts, one 
observes self-sustained current oscillations, see Figs. 29c,d. The shape of these oscillations 
as well as their frequency depends on the type of contacts. 

Nj:) < Nff^: Stable domains associated with a saw-tooth current- voltage characteristic do 
not occur. The behavior of the sample is either dominated by current oscillations, see 
Figs. 30b, c or a rather smooth stationary current-voltage relation, see Fig. 30d. According 
to the A''3/)L-criterion (182), the latter case dominates for short samples and lower doping. 
Here it seems to be crucial that the inhomogeneous field distribution essentially takes 
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Fig. 30. Same as Fig. 28 for a low doping density Nd = lO^^/cm^. (b) and (c) show current oscillations at a fixed 
bias. 

values from a range where dv/dF is small, suggesting large values of N^ . (Similar 
behavior is found for a = lO^A/Vcm^ if the number of periods is reduced.) 

Experimentally, this scenario can be verified by varying the electron density by optical irra- 
diation [202,203]. Transtions between oscillating and stationary domains can also be provoked 
by changing the lattice temperature [204]. This effectively alters the shape of the j{F) relation 
and thereby via Eqs. (179,180) the critical doping densities. There are several further aspects 
complicating the picture sketched above, which are not addressed in the discussion given here: 

Real superlattices are not perfect structures, where the properties of each well are repeated 
exactly. In contrast there will be fluctuations in the doping density, the barrier and well width 
as well as the material composition from period to period. Some information can be obtained by 
X-ray analysis and indeed it has been possible to relate some global properties to the extend of 
disorder [195,205] in the respective samples. By extensive simulations it has been shown that the 
presence of such disorder inside the superlattice affects the behavior significantly [187,206]. The 
nature of the current oscillations, the actual shape of the domain branches, as well as critical 
doping densities can be affected by the amount of disorder. Frequently, one observes a direct 
correlation between some global current signals and the actual realization in a particular sample 
[159,207]. 
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Recently, an additional S-type current- volt age characteristics has been found in strongly cou- 
pled superlattices due to electron heating [208]. The combination of N-type and S-type negative 
differential conductivity may provide additional interesting effcets. 

In all calculations performed here the transport model for sequential tunneling has been ap- 
plied, which provides rate equations between the quantum wells. For strongly coupled superlat- 
tices it is questionable if the electrons can be confined to accumulation layers extending over a 
few wells. Therefore it is not clear, in how far these stationary domain structures persist. The al- 
ternative is to start from the miniband model. Such calculations have been performed in [209,210] 
using the drift velocity from the relaxation time model. A more microscopic approach can be 
performed within the hydrodynamic model [167,211]. Nevertheless one has to be aware that 
the miniband transport model becomes questionable for large field strengths (compare Fig. 7), 
which are typically reached within a high-field domain. Thus a quantum transport calculation 
would be desirable to clarify the situation. Except for the stability analysis of the homogeneous 
state with respect to spatial fluctuations in [138] I am not aware of any quantum transport 
simulations concerning inhomogeneous field distributions in superlatticesP^. Thus it remains an 
open question in how far quantum effects modify the behavior discussed in this section. 



^^ Very recently stationary inhomogeneous field distributions could be modeled for short superlattices within the 
quantum transport model discussed in Section 4. First results are in qualitative agreement with those from the 
sequential tunneling model [212]. 
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6. Transport under irradiation 

In this chapter we consider superlattice transport under irradiation by an external microwave 
field with frequency i7/27r. In this case a further energy scale, the photon energy HQ of the 
radiation field, comes into play. For frequencies in the THz range (1 THz= 4.14 meV) this energy 
is of the same order of magnitude as typical miniband widths, scattering induced broadening, 
and the potential drop per period. This provides an interesting field to study various types of 
quantum effects. 

Transport under irradiation has first been studied theoretically within the simple Esaki-Tsu 
model. In this context it was shown that negative differential dynamical conductance [5,213] 
occurs. In the rectified response replica of the current peaks appear at field strengths obeying 
eFd = eFpca.kd + ihfl (with i £ Z) indicating the quantum nature of the radiation field [100,214]. 
Furthermore absolute negative conductance (i.e. a negative current for positive bias) is possible 
under certain conditions [215,216]. 

With the development of the free-electron laser as a high power THz source it became possible 
to study these effects experimentally and indeed the photon-assisted replica of current peaks 
[217,218] as well as absolute negative conductance [31] were observed. The superlattices used in 
these experiments exhibited rather small miniband widths, so that the application of miniband 
transport (as done in the theories mentioned above) is questionable as discussed in chapters 3 and 
4. Nevertheless, these findings could be explained qualitatively [31,219-221] within the standard 
theory of photon-assisted tunneling [222-224] as well, which is applicable for sequential tunneling. 
A quantitative description of these experiments [120] was possible within the sequential tunneling 
approach described in Sec. 3.3. 

Photon-assisted peaks in the current-voltage characteristic could also be observed in strongly 
coupled superlattices [225] although the results are less clear in this case. Furthermore the 
reduction of current due to the irradiation could be nicely demonstrated in experiments for 
superlattices with large miniband widths [226], which gave an excellent agreement with the 
simple miniband models mentioned above. 

In this section the basic ingredients of the transport theory under irradiation are reviewed. 
The first subsection deals with the simple quasi-static response, which holds for low frequencies. 
This will be the basis for the discussion of quantum effects in the subsequent subsections. The 
main results within the miniband transport will be reviewed in the second subsection. Here a 
form will be chosen which simplifies the comparison with sequential tunneling which is discussed 
in the third subsection. 

Throughout this chapter we consider a homogeneous electric field F{t) along the superlattice 
structure which can be separated into a dc-part F^c and a cosine-shaped time dependence with 
amplitude Fac ^ 0, i.e., 

F{t) = Fdc + F^,cos{nt). (185) 

The transport problem has been considered for Fg^c = in the preceding chapters where the 
relation /dc(^dc) was obtained. Now we are looking for periodic solutions with period 2it/Q, 
neglecting transient effects as well as the possibility of aperiodic behavior. Then the general 
current response can be written as 
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I{t) =Iq + Y^ [If"' cos{hnt) + Ifi'' sm{hnt)] . (186) 

h=l 

In general, four different aspects of transport under irradiation can be identified: 
(i) The rectified response Iq is considered in many different experiments. It is easily accessible 

and can be used for the detection of high frequency signals, 
(ii) The active current If"' provides the direct interaction with the irradiation field. If dlf°'/ dF^^c < 

we observe gain at the given frequency fi. 
(iii) The reactive current /f ° describes the response out of phase, which can be described by 
an inductance (/f° = Uac/LQ) or a capacitance (/f ^ = — CTit/ac) in standard circuit the- 
ory. There are two possibilities to define inductive and capacitive effects: (i) One assumes 
that C and L are always positive. Then positive/negative /f ^ is referred to as an induc- 
tive/capacitive effect, respectively [227]. (ii) One regards the low frequency limit. It will be 
shown later that typically If ^ oc 0, for low frequencies. This resembles the behavior of a 
capacitor which can either be positive or negative [228]. 
(iv) Harmonic generation Ih for h > 2: These terms describe the occurrence of higher harmonics 
and can be used to generate higher frequencies [229-231]. If i^dc = only the odd multiplies 
h are present for symmetric structures with Idc(-F) = — /dc(— -^)- 
In this chapter we restrict ourselves to homogeneous field distributions with the time depen- 
dence (185). In this case the current is homogeneous over the superlattice direction and no charge 
accumulation inside the structure occurs. Such complications can be treated within the general 
formalism discussed in [232,233]. Domain formation effects in superlattices under irradiation 
have been studied in [234]. If higher harmonics are present in the time dependence of the field, 
the superlattice may act as a rectifier [235]. Furthermore the response of an external circuit is 
neglected here. A detailed discussion of the latter issue can be found in [236]. 

6.1. Low frequency limit 

In the range of radio frequencies (say fi ^ 1 THz) the frequency 0, is slow with respect to 
the internal degrees of freedom, such as carrier heating or AE/h (where AE describes typical 
energy scales of the transport problem). Than one can assume, that the current follows the field 
instantaneously: 
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This shows that all terms /^"^ vanish in the radio-frequency limit. Furthermore, we obtain the 
following expressions in lowest order of the irradiation field: 
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They will be compared with the results discussed in the next subsections. 



6.2. Results for miniband transport 

In the miniband transport model the time dependence of the electric field enters the Boltzmann 
equation (35) which complicates the problem tremendously. The easiest way to deal with this 
situation is the relaxation-time approximation. Assuming Tm = Tg = t, Eqs. (54,56) can be 
solved for an arbitrary time dependence of the electric field. One obtains [215,237] 
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Let us now consider the field dependence (185). The crucial parameter in the following will be 
the ratio between the ac-field strength and the photon energy 



a 






(193) 



which will appear in most of the following results as the argument of the integer Bessel functions 
Jji- In order to avoid confusion with the symbol J for the current density, the results are given 
in terms of the current / = AJ in the following. Then one finds: 
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From Eq. (57) one can identify 

2eA|Ti|ccq(Ai,T) 
IdcieFd) = ^-fi Im ■ 

with r = h/r. Furthermore we define 
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where the Kramers-Kronig relation, connecting the imaginary and real part of h/{T — i<S), has 
been applied. Then we can evaluate the components of Eq. (186) in the form 



Y^ Je{a) {Je+h{a) - Je-h{a)) K {eF^^d + Ih^) . 



Tsm 



(197) 
(198) 
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Let us first consider the rectified response Iq, which is a sum of several dc-curves (as shown 
in Fig. 5b shifted by integer multiples of the photon energy. This explains the occurrence of 
photon-assisted peaks at biases eF^cd ~ T + ihrt as shown in Fig. 31a. If the coefficient for ^ = 
becomes small, i.e., close to a zero of Jo (a) (the first zero occurs at a = 2.4048 . . . ) the terms 
with ^ = ±1 dominate, which can provide absolute negative conductance for hVi > P as shown 
in Fig. 31b [216]. 

For small ac-fields eFg^d <C hO- we may set Jo(a) ~ 1, J±i(a) ss ^eF^^d / {2MI) , while the 
higher order Bessel functions are approximately zero. Thus J/j ~ for h > 2 and 
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Fig. 31. Rectified current-field relation from Eq. (197) within the Esaki-Tsu model for a — 1.2 (a) and a — 2.4 
(b). The thick line gives the total current /o(Jdc)- The thin lines show the contributions (Jo(a))^/dc(eFdcd) (full 
line), {J±i{a))^Idc{eFdcd ± hQ) (dashed line), and {J±2{a))^ Idc{eFdcd. ± 2hQ) (dotted line). Parameters F = 4 
meV, hQ = 6 meV. 
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which can be viewed as a resistor l/gdyn and a capacitor q in parallel yielding a complex 
admittance z~^ = QAvn + i^Ci (small letters indicate quantities per period and the engineering 
convention I{t) oc e is used here). Similar expressions were derived in Ref. [5] for the model 
with two scattering times. Both ^(dyn and q may take positive or negative values. Furthermore 
note that q is not the sample capacitance which has to be added in parallel as well, but originates 
from a quantum effect. (For small frequencies q ~ ^ — ,X vanishes in the limit /i ^ 0, if the 
functions Idc{£) and K{£) are kept constant.) 

Eq. (200) shows that the derivative in the low- frequency response (189) is replaced by a finite 
difference on the quantum scale. A straightforward calculation for the Esaki-Tsu model gives 
5dyn(^) < for |eFdc(i| > i/F^ -\- [hO,y. In this range the superlattice structure can provide gain. 
Nevertheless one has to note, that this occurs in the range of negative differential conductivity, 
where the homogeneous field distribution is typically unstable as discussed in section 5. 

In the limit of small scattering F ^ 0, the functions I^d^Fd) and K{eFd) vanish unless 
F ~ 0. From Eqs. (197,198,199) one finds for Fdc ~ all components lQ,If^^,If^^ of the current 
vanish for Jo(a) = 0. This can be interpreted as a dynamical localization [6] or the collapse of 
the miniband [7]: for a certain strength of the irradiation field the periodic structure does not 
conduct any current. In addition, a finite conductivity appears at eF^^-d ~ iKl opening up new 
transport channels, which are not present for a = in the limit of F ^ 0. 

All these results have been obtained within the simple Esaki-Tsu model. Calculations within 
the energy balance model [238] provide similar results. If the dielectric relaxation is included, 
chaotic behavior [239] as well as a spontaneous generation of dc current [240] can be found. 
Recently, Monte-Carlo simulations of the Boltzmann equation under THz-irradiation have been 
performed, too [241]. 
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6.3. Sequential tunneling 

For sequential tunneling between two neighboring wells {m and m+1) Eqs. (197,198,199) hold 
again with the dc-expression from Eq. (79): 

— i^(k,i?)i:;+i(k,ii; + eFd) 

K,.,^ -2^^ (202) 

X [npiE - firn) - npiE + eFd - Hm+i)] 
and the quantity 

KieFd) = - 4e J] jF^'^l^ / — - upiE - ^„)<(k,S)Re{G^^I(k,ii; + eFd)} 



k,i.,/i ''-'^ 



+Re{G^'^*(k,F)}nF(F + eFd - Aii+i)i^+i(k,F + eFd) 



(203) 



TT LS - eFd 



where G'm(k, F) equals Gm,(k, F — meFd) in the limit of decoupled wells. (The definition used 
here differs by an ?n,-dependent shift of the energy scale from the one used in Sections 3.3 and 
B.l.) These expressions have been derived in [223] for a constant matrix elements H'^'^ . A similar 
derivation is provided in Appendix C. Furthermore it is shown there that Eq. (197) also applies 
to the case of a field-dependent matrix element H'^'^ = eFdR'^'^, which is relevant for tunneling 
between nonequivalent levels. Note that the THz field also couples the different subbands via 
the term Rq'^ , which may cause further effects [242]. 

For sequential tunneling Eq. (197) has a simple interpretation. In the evaluation of the dc- 
current, eFd gives the energy mismatch between the levels in well m and m+1. Under irradiation 
photons of energy hQ can be absorbed or emitted during the tunneling process, which provides 
the energy mismatch eFjcd it ihO, for the absorption/emission of I photons during the tunneling 
process. The Bessel functions represent the probability that an ^-photon process occurs assum- 
ing a classical radiation field (i.e. neglecting spontaneous emission). Eq. (197) then consists of 
the weighted sum of all possible photon-assisted tunneling processes. Similar results have been 
obtained in [243], where a strictly one-dimensional tight-binding lattice coupled to a heat bath 
has been considered (the approximation of incoherent tunneling dynamics applied there corre- 
sponds to sequential tunneling). An extended discussion of the methods used there can be found 
in [244]. 

As discussed in section 3 the structure of the first peak at low electric fields eFd ~ F is 
similar for miniband transport and sequential tunneling. Therefore the discussion for miniband 
transport given in the preceding subsection holds for sequential tunneling as well. In addition, 
Eqs. (202,203) also describe the current peaks at resonances between different levels (a,b) in 
neighboring quantum wells. There one typically finds, see Eq. (83) 

\zj'ab\2 -pbyCS 

/dc(eFd) = eNoA ^ ^ ' -— -r-73 ,^. ^ , ,„ . (204) 

dcv ^ ^ f^ (eFd -F F« - F*)2 + (F'''^f^/2)2 ^ ^ 

With the Kramers-Kronig relation one obtains 
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K{eFd) = -cNdA ^K^^ , ^ , ^y"' ^^1 f.? ^ , ,, (205) 

where the field dependence of H has been neglected. 

In the following some results are presented for the superlattice structure studied experimentally 
in [31,120,219] (15 nm wide GaAs wells, 5 nm Alo.3Gao.7As barriers, doping density Nd = 
6 X 10^/cm^, cross section A = 8/um^). Results for this structure have been presented already in 
Fig. 12, where the complicated temperature dependence was discussed. Here we will focus to the 
behavior under irradiation. For these calculations a constant electron temperature Tg = 35 K is 
assumed, which provides best agreement around the first maximum. Due to the presence of the 
irradiation field, it seems realistic that the electron gas does not even reach thermal equilibrium 
for vanishing dc-bias. 

6.3.1. Rectified THz response 

In Fig. 32 results are shown for the rectified current response loleFdcd) under difi'erent 
strengths and frequencies of the irradiation field. Both from experiment and theory one observes 
a range of absolute negative conductance (i.e., /q < for F^^ > 0) for low biases. Furthermore 
photon-assisted peaks are visible at field strengths eF^^d ~ e-F^^^^d + (.Ml. These findings are 
in qualitative agreement with the discussion of Fig. 31. Quantitative agreement between the- 
ory (Fig. 32a) and experiment (Fig. 32b) is found for h^ = 6.3 meV (1.5 THz) for different 
strengths of the laser field. The low-field peak occurs at f/jir = A^-^max^ ~ 20 mV correspond- 
ing to direct tunneling. Photon replicas can be observed at f^ ~ f/jir + NMl/e = 83 mV and 
U ~ L^dij. -|- 2NhQ/e = 146 mV. For low bias and high intensities (a = 2.0) there is a region of 
absolute negative conductance [31], which will be discussed in the following. 

In Fig. 32d the laser intensity has been tuned such that maximal absolute negative conductance 
occurred for each of the different laser frequencies. Then one observes a minimum in the current 
at [/ ~ — C^dir + NhQ/e which is just the first photon replica of the direct tunneling peak on 
the negative bias side. This replica dominates the current if the direct tunneling channel is 
suppressed close to the zero of Jo (a) in Eq. (197), i.e., a « 2.4, as used in the calculation of 
Fig. 32c. Both the theoretical and experimental results show that absolute negative conductance 
persists in a wide range of frequencies but becomes less pronounced with decreasing photon 
energy. In the calculation absolute negative conductance vanishes for hQ < 1.8 meV which is 
approximately equal to hQ < eF^^^xd- (The latter relation has been verified by calculations for 
different samples as well.) For hQ = 5.3 meV a smaller value of a = 2.1 (thin line) agrees better 
with the experimental data (in the same sense the value a = 2.0 agrees better for hQ = 6.3 
meV, compare Fig. 32a. This may be explained as follows: If strong NDC is present in doped 
superlattices, the homogeneous field distribution becomes unstable and either self-sustained 
oscillations or stable field domains form as discussed in section 5. Then the current-voltage 
characteristic deviates from the relation for homogeneous field distribution, where U = NFd, 
and typically shows less pronounced NDC. Therefore maximal negative conductance is observed 
at a laser field corresponding to a value of q < 2.4, where the NDC is weaker and the field 
distribution is still homogeneous. The presence of an inhomogeneous field distribution could also 
explain the deviations between theory and experiment for U > 150 mV. Quantitative agreement 
between theory and experiment regarding the rectified response was also obtained for a different 
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Fig. 32. Rectified current response for tlie superlattice of [31] displaying absolute negative conductance, (a) 
Theoretical results for hfl = 6.3 meV and different field strength eF^cd = ahQ of the irradiation, (b) Experimental 
results for hfl = 6.3 meV and different laser intensities increasing from the top to the bottom. The actual values 
_Fac inside the sample are not accessible, (c) Theoretical results for a = 2.4 and different photon energies. The 
thin line depicts hi} = 5.3 meV and a — 2.1. (d) Experimental results for different photon energies. The laser 
intensity was tuned to give maximum negative conductance. (From [120]). 
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Fig. 33. Dynamical response for the superlattice of [31]. (a) Current-field relation without irradiation, (b) Current 
response, full line /i°'', dashed line /f" for fixed eFdcd — 8 meV and hfl — 5 meV. (c) dito for fixed eFacd — 1 
meV and ftfi = 5 meV. (d) dito for fixed e_Fkcd = 1 meV and eFdcd = 8 nieV. 

superlattice structure, where up to seven photon replica of the first current peak could be 
observed [245]. 

6.3.2. Negative dynamical conductance 

A semiconductor element is able to give gain at the given frequency if d7™^/dFac < 0. In 
the low frequency limit (subsection 6.1) one obtains Ii^rf = d-Idc{F) / dF x Fac if F^c is not too 
large. This means that gain is related to negative differential conductance in the static current- 
field relation I^dF). As such a situation is typically unstable with respect to the formation of 
inhomogeneous field distributions, it is difficult to apply this concept for a real device. Thus it 
would be of interest to have a system exhibiting dli°^/ dFg^c < in the THz frequency range 
considered but a positive differential conductance at fi — > 0. Inspection of Eq. (200) for the 
standard expressions (195,204) shows that 1^°^ is always positive as long as d/dc/dF > holds. 
Thus, it is a nontrivial task to find the opposite case in a semiconductor superlattice. 

Fig. 33a shows the calculated current-field relation Idc{eFd) for the superlattice structure 
discussed here. It exhibits a kind of plateau in the range 2 meV< eFd < 9 meV, due to the 
presence of impurity bands which is observed experimentally as well (see Fig. 12). Experiments!^ 
show that this plateau is almost unchanged for lattice temperatures between 4 K and 35 K. At 



' private communication from Stefan Zeuner 



84 



eFd = 8 meV we have positive differential conductance, but the finite distance derivative for 
hi} = 5 nieV is clearly negative as indicated in the figure (see also [246]). Thus /J°^ will be 
negative as long as eFg^cd is not too large and the terms with i > 1 become important in 
Eq. (198). From Fig. 33b one obtains a negative dynamical conductance ^fdyn = — 33/uA/V 
at hQ = 5 meV. Due to the higher order terms the dynamical conductance becomes positive 
for ac-field strengths larger than 8.5 mV per period. Fig. 33c, d show that negative dynamical 
conductance persists over a wide range of dc-bias and frequency. Nevertheless one must note 
that for eF^cd > 9 meV the static conductance d/dc/ d(Fd) becomes negative yielding domain 
formation as observed experimentally (see Fig. 12). 

Unfortunately, the negative dynamical conductance of the device considered here is compen- 
sated by the contact resistance. Measurements of the temperature dependent conductance of 
the sample yield values up to 20;uS (around T = 30 K) at zero bias [120]. As a part of the 
resistance is from the superlattice itself one may conclude that the contact resistance Re is def- 
initely smaller than 50 ki7 in the sample. In comparison the dynamical resistance for 10 wells is 
-^dyn = lO/ffdyn = "300 kil, which sccms to dominate. But capacitive effects have to be taken 
into account: The sample capacitance per period is given by Cg = Aeeo/d ^ 46 fF. This gives a 
total impedance 



1 ^^ ,^ 1 - injcs + Ct)/gdyn 

ffdyn + i^(Cs + Q) "" ^'' I + {Q{Cs + Ct) / gdyn)"^ 



Z = Rc + N -— ■ - = i?, + i?dyn , , ,^; ', ; ^2 (206) 



As Cs/gdyn = —1.4 ns, the negative dynamical resistance will be compensated even by a small 
contact resistance at THz frequencies. Thus, significant larger values of (/dyn (i-e. higher current 
densities) are necessary for the observation of gain. Nevertheless, the effect discussed here is quite 
general and therefore gain should be observable in superlattice structures with special shapes 
for the dc-characteristic. 

6.3.3. Tunneling capacitance 

Now we want to investigate the reactive current If^ from Eq. (199). As can be seen in Fig. 33d 
as well as in Fig. 34b) If™ oc fi for low frequencies in the linear response region (i.e., small 
irradiation fields). This can be interpreted as a tunneling capacitance 

/r ~ -ctVLF^d (207) 

as shown in Eq. (201). In Fig. 34a the reactive current /f° is displayed as a function of -Fdc 
using a small irradiation field and a small frequency, so that Eq. (207) holds. The quantity If™ 
shows a very characteristic behavior around the a ^ b resonance with a minimum at eF^cd = 
50meV, where the differential conductance is positive. This behavior can be understood within 
the approximation (205), yielding 

eh d^KjeFd) ^ ^, ^^^|^a6,2 2(ei^rf + ^" - E") [{eFd + E- - E^f - 3{T'^-^ /2f] 



2 d(eFd)2 [{eFd + E<'-Eby + {Tb'^f^/2yf 

(208) 

which provides just the structure observed around ei^dcC ~ 50 meV. From Fig. 34b one obtains 
ct = 0.67 fF at eFdcd = 50meV, which is about 1.7% of the sample capacitance c^. It would 
be interesting if this tunneling capacitance can be measured. Here it may be useful to study 

85 




T \ 1 \ 1 \ 1 \ 1 \ 1 1 

10 20 30 40 50 60 



T \ 1 1 1 r 

2 4 6 



~\ \ 1 

8 10 



eFdcd [meV] 



h^ [meV] 



Fig. 34. Dynamical response for the superlattice of [31]. (a) Current response (full line Ii°^, dashed line /f") for 
fixed hfl — 1 meV, (b) for fixed eF^cd = 50 meV. In both cases eFi^^d = 0.1 meV. 




eFacd [meV] 



eFacd [meV] 



Fig. 35. Generation of harmonics |Ji| (full line), |73| (dashed line), and lis] (dotted line) for two different frequencies 
for the superlattice of [31]. As Fdc = 0, the even harmonics vanish. 

superlattice structures with higher doping and larger coupling, which enhances the ratio between 
ct/cs-, as can be seen from Eq. (208). 



6.3.4. Harmonic generation 

In order to investigate harmonic generation the quantities |//i| = \ {I'l^^Y + i-^h^^Y have been 
plotted for i^^c = in Fig. 35. The basic frequency |/i| dominates for low irradiation fields, 
while the other quantities vanish like I/3I ~ (eFac^)^ and \I^\ ~ {eFi^d)^ as expected from the 
low- frequency limit (190). The current of the third harmonic is in the range of the dc-current 
at the first peak indicating strong harmonic generation. In the range eFg^cd > 50 meV the 
a ^ b resonance becomes of importance which can be understood from the classical behavior 
Idci&Fi^cdcos{Qt)). This provides a possibility to probe the second resonance even if it is not 
directly accessible due to sample heating. A similar idea to probe a peak by its response to an 
applied frequency has been performed in [247]. 
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6.4. Discussion 

In this chapter it has been shown that transport under irradiation is essentiaUy governed by 
Eqs. (197,198,199). This scheme holds both for miniband transport within the relaxation time 
approximation and sequential tunneling for constant coupling matrix elements, albeit with dif- 
ferent functions IdcisFd) and K{eFd). Therefore the question arises, if this structure might be 
general. Although this issue has not been settled finally, I believe that this is not the case. If one 
considers next-nearest neighbor tunneling processes within the model of sequential tunneling, it 
becomes obvious, that photon-assisted peaks occur at field strengths 2eFdcd = 2eFpeakd + ihi^ 
in contrast to the structure of Eq. (197). (This might provide an interesting tool to investigate 
tunneling processes where tunneling occurs between wells separated by more than one barrier 
[62,127,126].) The same behavior is found in a miniband model containing higher Fourier compo- 
nents in the band structure [248]. Therefore it seems that the structure of Eqs. (197,198,199) is 
limited to next-neighbor tunneling, or a cosine-shaped band structure E(q), respectively. Further 
deviations may occur if the THz-field causes additional heating of the electron gas| ^'^ |, which had 
been neglected in the models applied here. 

The generic structure of the functions I{eFd) and K{eFd) is given by Eqs. (195,196,204,205). 
These imply the following typical effects under irradiation: 

- Photon-assisted peaks in the rectified current response at characteristic field strengths eFd « 
eFpeak*^ + ^^^ with integer values £ for next neighbor coupling. 

- Absolute negative conductance if the normalized ac- field strength a is close the zeros of Jo (a) 
and Ml > T. 

- Gain in most of the region of negative differential conductivity d/dc/ d-F < 0. 

- A quantum capacitance with the characteristic dependence (208) on the dc-field close to the 
resonances. 

- Generation of higher harmonics. 

These effects rely on the structure of Eqs. (197,198,199) and should therefore (at least approxi- 
mately) hold both for wide and narrow minibands. In addition, it has been shown that gain in 
the THz range is also possible in the region of positive differential conductance for appropriate 
shapes of the Idc{eFd)-icelation. 



^^ This was pointed out by A. A. Ignatov (private communication, 1998) 
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7. Summary 

In this review the transport properties of semiconductor superlattices have been analyzed. 
Strong emphasis has been given to the microscopic modehng of the stationary transport for a 
homogeneous electrical field as well as the formation of inhomogeneous field distributions leading 
to stationary field domains and self-sustained current oscillations. 

The three different standard approaches, miniband transport, Wannier-Stark hopping, and 
sequential tunneling, have been reviewed in detail. Although the concepts applied are quite 
different, each approach provides negative differential conductivity for sufficiently high electric 
fields. While miniband transport and sequential tunneling provide an Ohmic behavior for low 
electric fields together with a maximum of the current around eFd ~ T, Wannier-Stark hopping 
fails in the low-field region. In particular, good quantitative agreement with various transport 
measurements in weakly coupled superlattice structures has been obtained within the sequential 
tunneling model, both with and without irradiation. 

The relation between the standard transport approaches could be identified by considering 
quantum transport based on nonequilibrium Green functions. It has been explicitly shown that 
the equations used for these simplified models can be obtained from the quantum transport 
model by applying various approximations. This justifies each of these approaches and sheds 
light on the respective ranges of applicability sketched in Fig. 7. Good quantitative agreement 
was found between self-consistent solutions of the quantum transport model with each of the 
simplified approaches in their respective ranges of applicability. 

We have shown how different aspects of nonlinear pattern formation in semiconductor super- 
lattices, such as the formation of stationary field domains as well as self-sustained oscillations, 
can be understood by the properties of traveling fronts. Fig. 27 shows that depletion and accu- 
mulation fronts travel with different velocities, which depend on the total current acting as a 
global coupling. This situation is similar to the Gunn diode albeit the discreteness of the super- 
lattice structure allows for stationary fronts in a finite interval of current, when the front becomes 
trapped. As this trapping may occur in each of the wells, one obtains several branches (of the 
order of superlattice periods) in the global current voltage characteristics, which exhibits a typ- 
ical sawtooth shape. The calculated current-voltage characteristics as well as the self-sustained 
oscillations are in reasonable quantitative agreement with experimental findings. Furthermore, 
the conditions for the occurrence of self-sustained oscillations have been discussed. 

Under strong irradiation by a THz field, photon-assisted peaks appear, and absolute negative 
conductivity is possible. Both effects are predicted in a similar way by the sequential tunneling 
model and miniband transport. It has been demonstrated that there is a possibility for gain even 
in the region, where the low-frequency conductivity is positive (i.e. no field-domain formation 
effects are expected). In addition, the tunneling processes are connected with characteristic 
variations in the capacitance of the structure, which have not been observed so far. 



8. Outlook 

Even 30 years after the proposal by Esaki and Tsu, semiconductor superlattices continue to 
be a hot topic of ongoing research. Presently the following directions stand out: 
Superlattices with lower dimension: In this work superlattice structures have been consid- 
ered, where a free electron behavior is present in the two directions parallel to the layers. 
Presently, the first experiments have been performed where these lateral directions are con- 
fined. In Ref. [249] the measurement of the conduction through a stack of 50 InAs quantum 
dots was reported (see also [250]). Such structures can be regarded as a superlattice structure 
consisting of zero-dimensional boxes. Negative differential resistance was observed recently in 
a superlattice formed by quantum wires, fabricated by the method of cleaved edge overgrowth 
[251]. Due to the restricted phase space perpendicular to the transport direction, scattering 
should be strongly reduced in these structures, and strong effects related to the miniband 
structure are likely to be observed. Theoretical approaches to transport in such structures can 
be found in Refs. [252,253] for quantum box superlattices and in Ref. [254] for a superlattice 
formed by quantum wires. 
Self-sustained current oscillations: As discussed in section 5, traveling field domains cause 
self-generated current oscillations in superlattices. While the first experiments reported fre- 
quencies in the MHz range [19] at 5 K for a weakly coupled superlattice, frequencies up to 150 
GHz at could be observed even at room temperature in specially designed superlattices with 
a large miniband width [20]. While it may be difficult to increase the fundamental frequency, 
the use of higher harmonics can be helpful for possible devices in the THz-range, although the 
prospects of this approach are under debate [255]. 
Field domains in strongly coupled superlattices: The formation of field domains is usu- 
ally considered within the model of rate equation presented in Section 5.1. This model is 
justified for weakly coupled superlattices, where the electrons can be considered to be lo- 
calized in single quantum wells. This localization becomes questionable for strongly coupled 
superlattices and it would be desireable to develop a quantum transport model, which takes 
into account both the inhomogeneity and the time dependence of the field distributions. (See 
also the discussion at the end of Section 5.) 
Chaos and nonlinear dynamics: If the superlattice is driven by an ac-bias, chaotic behavior 
is likely to occur due to the presence of incommensurable frequencies, which is well known 
for the Gunn diode [256]. Simulations for semiconductor superlattices [257,258] yield similar 
results. These chaotic oscillations could be observed experimentally [259,260] and provide a 
nice example, where many aspects of chaotic systems [261] can be observed. There are also 
reports indicating undriven chaos under dc bias conditions in superlattices [259], which is not 
understood yet. In addition, frequency locking associated with the the occurrence of a Devil's 
staircase has been observed [262]. These examples provide an interesting field to study general 
phenomena of nonlinear physics. 
Generation of THz signals: Since the original proposal by Esaki and Tsu, it has been tempt- 
ing to use the Bloch oscillator as a source for THz radiation. Unfortunately, only transient 
signals have been observed so far. Another possibility lies in the occurrence of a negative dy- 
namical conductance. Here the problem arises that domain formation is a competing process 
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in this region of parameter space. Two different scenarios have been proposed as a solution 
to this problem: In Ref. [246] it was proposed that for special shapes of the local current- 
field characteristic gain is possible even in the region where the low-frequency conductivity is 
positive, i.e., no field-domain formation effects are expected (see also section 6.3.2). A differ- 
ent possibility could be related to the fact that the low-frequency conductivity can become 
positive for large amplitudes of the ac fields [255], thus stabilizing the oscillation. It is not 
clear by now, if one of these effects may be useful to establish THz devices on the basis of 
semiconductor superlattices. 
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Appendix A. Sequential tunneling with density matrices 

As discussed in Section 3.3 the electrical transport in weakly coupled superlattices can be 
described by sequential tunneling. The idea of this concept is that the sates are essentially 
localized in single wells and the residual coupling causes transitions between neighboring wells. 
As this coupling (Ti for equivalent levels) is small, one can restrict the theory to the lowest order 
Tf, which provides essentially Fermi's golden rule. Already in Sec. 3.3 it was mentioned that 
scattering induced broadening is essential to recover the correct behavior for finite fields. This 
complication was treated within the theory of nonequilibrium Green functions which provided 
Eq. (79) as derived in Appendix B.l. 

An alternative way to treat quantum transport is density matrix theory which also gives an 
exact treatment of both quantum effects and scattering. A recent overview can be found in 
[140]. Density matrix theory has been applied to superlattice transport in [73] for the evaluation 
of currents between nonequivalent levels. To my knowledge no such calculations exist regarding 
transport between equivalent levels with identical particle densities. In this case the calculation of 
[73] provides zero current independent of the electric field as shown below. A possible resolution of 
this problem will be presented in this appendix and a form similar to Eq. (83) will be derived. This 
demonstrates the equivalence of both approaches and highlights the differences in performance 
of both approaches. 

A.l. The model 

We use the Hamiltonian (23,24) in the basis of Wannier states and restrict ourselves to the 
lowest band for simplicity. Furthermore impurity and phonon scattering is taken into account 
within the restriction of intrawell scattering. Then the Hamiltonin reads H = Hq + U + Mmp + 
Vphon with: 

Ho = y^{E^ + Ek-eFdn)ai{k)an{k) (A.l) 



n,k 



U-- 



J2Ti [al+i(k)a„(k) +aUk)G„+i(k)] (A.2) 



n.k 



Vimp= J2 n\4(k')«n(k) (A.3) 

n,k,k' 

V^phon = Yl ^pbiip)bnip) + Yl ^p'«n(k + p) [&n(p) + ^^-p)] ani^^) ■ (A.4) 

n,p n,k,p 

In the following we assume that phonon scattering is strong enough to establish thermal equi- 
librium. Thus a thermal distribution function with chemical distribution Hn can be assumed for 
each well. Furthermore, let us assume that impurity scattering is stronger than phonon scat- 
tering and dominates the broadening of the states, which simplifies the following calculations 
essentially. Correlations between the scattering matrix elements V^ in different wells are ne- 
glected. (These drastic approximations make sense as we are mainly interested in the structure 
of the theory, not in quantitative results. The inclusion of different scattering mechanism should 
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be possible in an analogous way.) In this case the self-energy within the Born approximation for 
impurity averaging is given by 



|2 



1 



Sn {KE)-Y,\ykk'\ E_E,, + neFd + iO+ 

k' (A.5) 

« - iTT ^ |VS,|^ 5{E - Ek' + neFd) « -iT/2 
k' 

where we neglected the real part and assumed that the scattering rate is energy independent for 
simplicity. Using Eq. (102) one obtains the current density from well n to well n + 1 

Jn^n+i = '(^°^ J")^ ^ ^Im{T,(aU,(k)a„(k))} . (A.6) 

k 

The task is the evaluation of the current in lowest order T^ with respect to the coupling (sequen- 
tial tunneling). Under these conditions the theory of nonequilibrium Green functions provides 
Eq. (79) which can be simplified to Eq. (83) for the self energy (A.5). 

A. 2. Density matrix theory 

The key point is the temporal evaluation of (a|^_|_-|^(k, t)o„(k, t)). Similar to Eq. (101) the 
dynamics is given by 

Y^(4+i(k,t)an(k,t))=([i^,a|;+i(k)a„(k)) 



= - eFd(ai+i(k)a„(k)) + ^ ((at,(k)a„(k)) - (4+i(k)a„+i(k)) 
+ E K'r(4+i(k')«n(k)) - ni,(al+i(k)a„(kO)" . 

k' 

(A.7) 

Here (and in the following) terms containing density matrices (ajj^_2(k')a„(k)) will be neglected 
as they provide terms of order Tf in the current. As stationary states will be considered, the 
time dependence is dropped in most of the density matrices. On the right side new expressions 
of the type (a„_|_]^(k')a„(k)) appear. Their temporal evolution is given by 

Y^K?r(4+i(k',t)«n(k,t)) = {Ek> -E,- eFd)y,\+i(al+i(k')a„(k)) 

+ T.^k%'^kTi^Uii^'>ni^)) - E^.l"n"r(4+i(k')«n(k")) (A.8) 

k" k" 

+ Ti [F,'?r(at(k')a„(k)) - V,T{aUiik>n+ii^))_ 
and 
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Y^n"fc'(4+i(k,i)an(k',t)) = {Eu - Ey - eFd)y,l,(al+i(k)a„(k')) 



+ E V^,'^4'^/rfc'(4+i(k")«n(k')) - E nV^fcl'(4+i(k)an(k")) (A.9) 

k" k" 

These equations have again to be solved in the stationary state. In the fohowing the key quantities 
of interest are the occupation /n(k) of the mode k in weh n and the corresponding polarizations 
P„(k) which provide the current. They are given by 

/„(k) = (a];(k)a„(k)) and P„(k) = (al+i(k)a„(k)) . (A.IO) 

In the following, the more complicated density matrices appearing in Eqs. (A. 8, A. 9) have to be 
related to these quantities. 

A. 2.1. Lowest order calculation 

In [73] only the lowest order terms within the impurity averaging process and the cou- 
pling have been considered for the evaluation of Eq. (A. 7). As the current is already of order 



Ti(ajj_|_-|^(k)a„(k)) only terms up to T\ or V will be taken into account. Due to impurity averag- 

/n+l-r/n+l _ |T/n4 



ing only the term VTIt Vr,t = \yu/t P remains then in Eq. (A. 8) and we find in the stationary 



state 

where the term iO^ ensures that the correlations vanish for t ^ —oo. (An alternative way to 
obtain the factor iO^ is the application of the Markov limit, see [140].) Using the same argument 
we obtain from Eq. (A.9): 

ni'W„»(k)a„(k')) = „^_Ej,ir,^,„, \VS,\'P.ik} (A.12) 

Inserting into Eq. (A. 7) gives 

-Ap{k,t) = -eFdPnik) + Ti [/„(k) - /„+i(k)] + irP„(k) (A.13) 

1 at 

where the approximation (A. 5) has been used. This equation has also been applied in [263] to 
study time dependent phenomena. The stationary solution yields 

^"('^) = -eFd + iT ^' ^^""^^^ " ^"+^^^^1 ^^-^^^ 

and we obtain the current density via Eq. (A. 6) 

T'^ 2r 

k 

Taking into account that p{E) = 2(for spin)/ A'^^S^E — Ek) is the density of states, this 
expression is almost identical with Eq. (83). Nevertheless there is a significant difference: While 
in Eq. (83) the transport is driven by the difference of the occupation at the same energy^ now 

93 



the transport is driven by the difference of the occupation of the state k in both wells. The 
latter difference becomes zero, if both wells have the same electron density. This has lead to the 
conclusion that no resonant tunneling peak occurs in weakly coupled superlattices for tunneling 
between equivalent levels [8], in contrast to the findings of section 3.3. Regarding tunneling 
between the ground state and excited states (which are typically empty) the respective formula 
provides a finite current which essentially agrees with the corresponding result of Eq. (83). The 
corresponding calculations have been presented in [73], where interwell correlations between 
scattering matrix elements were also taken into account. 

A. 2. 2. Improved treatment 

This problem can be circumvented by taking into account the last term in Eqs. (A. 8, A. 9) as 
well. Neglecting terms containing the coupling Ti (which provide terms ~ Tf in the current) the 
dynamics of V'^"^(ajj(k')a„(k)) is given by: 

-4^fe"l(4(k',t)a„(k,t)) = {Ek> - Ek)V,i,{aUk')an{k)) 

1 f A 16) 

+ E^fc"fc'^fc'^(«"(k")«n(k))-j;F,V<(at(k')a„(k")) 

k" k" 

By impurity averaging all terms V'^ vanish unless n' = n and k" = k oi k" = k'. Thus, we find 
in the stationary state: 

V^'f,{ai{k')an{k)) = -__^i__^|y«,|2^„„, [/„(k) - /„(k')] (A.17) 

Then the evaluation of Eq. (A. 8) yields 

^.\^'(4+l(k')On(k)) 

{K%'\'Pnik) + E^,_ll^^o+ KT\' [/™+l(k) - /n+l(k')] } 



Ek' - Ek- eFd + iO 

-1 
Ek'-Ek-eFd + iO+^^'''^ 



Ti ,...^i,o / 1 1 



+ eFd I ^fc"^' ' ' ( Ek' -Ek + iO+ Ek' - E^ - eFd + iO+ ) [^"+' ^^^ fn+ii^')] 



and Eq. (A. 9) gives 



(A.18) 



n".(4..(k)a„(kO) = ^,_^,_,^,^,o^ l^."..|^^»W 



Tl „.n ,2^1 1 



^fc^fcl^ P. _ P., . ^n+ ~ p. _ p.,_.p^^;n+ [f^^(^') - /«(k)] 



eFd ' '"" \Ek - Ek' + iO+ Ek - Ek' - eFd + iO+ 
This provides further terms for Eq. (A. 7), which becomes 
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(A.19) 



y^P„(k,t) = -eFdPnik) + Ti [/„(k) - /„+i(k)] + irP„(k) 



' E [E,,-E, + iO+ " E,,-Ek-eFd + m) '^^'^^''' [^"+^^^'^ " ^"+^^''^] 



eFd 



(A.20) 

Now we apply the relation l/(x + iO'^) = V{l/x} — iir6{x) and restrict ourselves to the imaginary 
parts. The isotropy in k-space gives /n(k) = fn{Ef;) and we find with Eq. (A. 5) the stationary 
solution 

P„(k) = ^/ ^^ {Ti [UEk) - fn+l{Ek) 

(A.21) 
eFd 2 L^^+i^^'^) - /«+i(^fc + e^^) - /n(^fc) + fniEk - eFd)] I 

and the current (A. 6) becomes 

Tf sr^ 2r 




1 \~~^ 



M ^^ (eF(i)2 + r2 

^ (A 22) 

fn{Ek) + fn{Ek - eFd) fn+l{Ek + eFd) + fn+liEkT 



In contrast to Eq. (A. 15) the current is now driven by the occupation difference taken at different 
values of Ek so that the total energy E^ — neFd is equal in both wells. This structure agrees with 
Eq. (83) and a finite current is found for identical occupation function fniEk) = fn+i{Ek)- This 
shows that the terms V^"^(a]i(k')a„(k)), which had been neglected in the derivation of Eq. (A. 15), 
are of crucial importance. They describe the internal correlations in the single quantum wells. 
These correlations correspond to the broadening of states in the Green function formalism, where 
they are taken into account by treating the probe energy E and the energy of the bare state E^ 
separately. As the broadening of the states is of crucial importance for the tunneling current, it 
becomes clear, that density matrix theory gives a wrong result if the corresponding matrices are 
neglected. 
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Appendix B. Derivation of the standard approaches 

In this appendix the relations between the quantum transport equations of section 4 and the 
standard approaches (miniband transport, Wannier-Stark hopping, and sequential tunneling as 
discussed in section 3) are examined. It will be shown, that the transport equations for the 
different standard approaches can be derived explicitly from the quantum transport model using 
various types of approximations. In each case the respective approximations can be justified 
within the range of validity of the given standard approach sketched in Fig. 7 and motivated in 
subsection 4.3. 



B.l. Sequential tunneling 

In the parameter ranges 2Ti ^ F or 2Ti <^ eFd the diagonal elements of Gn,m dominate (see 
subsection 4.3) and an expansion of Eqs. (153,154) in Ti is appropriate. In this way the formula 
for sequential tunneling (79) will be recovered as the leading order in Ti of the general Eq. (148). 

For Ti = we obtain Gmi^^O^^ ^) — '^mi,m2^mi C^' ^) which are determined by 



E-Ek + eFdm - S^*(k, E) ) GZ\K E) = 1 (B.l) 
and 

G^%E) = &^\k,E)t<{k,E)G^^^{k,E) (B.2) 

where the self-energies S^ are evaluated applying the Green- functions Gm- These equations 
decouple in the well index. As no current flows in this case, one obtains the equilibrium solution 

(?< (k, E) = ii„(k, E)nF{E - ^„ + meFd) (B.3) 

S< (k, E) = [t^^K E) - S^t(k, E)\ npiE - ^i^ + meFd) (B.4) 

with the spectral function 

i„(k, E) = i [gZ\K E) - G^\k, E)] = -2Im [gZ'% E)] = 2Im {^^^^(k, E)] (B.5) 

To first order in Ti, Eq. (153) gives together with Eq. (B.l): 

G-*±i,„(k, E) = (5^*±i(k, E)T^&:^'% E) + 0{Tl) (B.6) 

G-*^(k, E) = (5-*(k, E) + 0{Tl) (B.7) 

Gl%iKE)=OiT^) for \m-n\>2 (B.8) 

and Eq. (154) gives 

96 



+ GZXiiK E)t<_,,{k, E)G:^l,{k, E)T,Gt\k, E) + 0{tI) 



--Ti 



G':,%iiKE)G<{k,E)+G<^,{h,E)GT{KE)\ +0{Tt) (g_g^ 

GZ\i% E)Ara% E)nF{E -ijira + mcFd) 



+ A„+i(k, E)nF{E - fim+i + (m + l)eFd)G''^''{k,E) + OiT{ 



^2\ 



Then the current is evaluated by Eq. (148) 



2e^-^|ri|2 I- dE ~ ~ (B.IO) 

=l"Z^~/r / ^^™+i(k,^)Ar«(k,^) 

X [npiE - fXra + meFd) - npiE - Hm+i + {m + l)eFd)] + 0{Tf) 
which is just the expression (79) used for sequential tunneling. 

B.2. Miniband conduction 

If 2|Ti| ^ T,eFd, the states in a semiconductor superlattice are essentially delocalized, as 
shown in subsection 4.3. In this case it makes sense to work in an extended basis like the Bloch 
states q. The key point is the idea, that the occupation of these states can be treated as a 
semiclassical distribution function f(q, k) neglecting quantum mechanical correlations. We will 
show that the Boltzmann equation for f(q, k) as well as the formula for the current density (see 
Sec. 3.1) can be derived from the full quantum transport model under this assumption. 

We consider a superlattice with a homogeneous electric field F and a homogeneous carrier 
distribution. In this case it makes sense to use a local energy scale £ = E — e(j)n which refers to 
the bottom of the respective quantum well (here e(pn = —neFd holds). We define 

k,£:-eFd^— j (B.ll) 

and S„^m(k, iS) is the same way. Now, G'm^„(k, iS) only depends on the difference m — n and 
S^(k,iS) = X](k, iS) due to the homogeneity of the system under stationary transport. Therefore 
it is helpful to define the spatial Fourier transform 

G{q,k,£) = Y,^-''^''Gn+hA^,£) (B.12) 

h 
This corresponds to the Fourier representation 

G{q,k,£) = Y.l[ ^t eK^+-^4^)ie-'^'^'^G„+;,,„(k;t + t2,t2) (B.13) 

h •' 
which is the special case for a homogeneous electric field of the general gauge invariant version 
used in section 7 of [39] to obtain gauge invariant quantities [144]. The respective diagonal 
element of the density matrix is defined by 
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^^'^'^^ = ^/ d^^^('?'k,^) = ^/ d^ E«~''"'^n+M(k'^) (B.14) 



Applying Eq. (155) we obtain with these definitions 

eFd 



d - 
ieF—G<(q,k,£) = 2[sm(qd)Ti 
oq 



G^ q,k,£ + 



G<{q,k,£ 



eFd 



+E' 



-iqhd 



s< ( k,f + ^] Gi^^{k,£) - Gfo(k,Oi:< (k£-^ 



(B.15) 



=iadv 



k f _ !!^££] _ f;rct /^ -^ /ieFd\ 



J-S'-^*^k,f + ^-JjG<o(k,f) 



which is still exact. In the same way Eq. (153) yields: 



£-Ek + '-^-^ ]G''\q,k,£) - e'l'^TiG'"' 



2 dq 



hK£ + 



eFd 



-'^'^T^G'^^{q,k,£ 



eFd 



1 + ^ ,-^'iM ^ret j^k, f + ^^ G^4(k, f ) 



(B.16) 

Assuming that the self-energy does not depend strongly on £\^ within the energy scale eFd, 
i.e. S'^* {K£ + heFd/2) ^ t.""^ {k,£), this equation is solved by 

^^d^G'^\q,k,£) 



G'^'''{q,K£) 



+ O { {eFdy 



(B.17) 



£-Ek-2Ticos{qd)-T.'''^{k,£) ' " |^^"^ ' d£^ 

On the energy scale 2|Ti| the Green function G^'^^{q,k,£) essentially resembles the free particle 
Green function 1/{E -E^- 2Ti cos(gd) +iO+) if S''<='^ ~ T < 2|ri| and eFd < 2|Ti| (so that the 
last term is negligible), which is just the condition (163). 

Now we integrate both sides of Eq. (B.15) over energy £, so that the Ti terms on the right-hand 
side cancel each other. In the terms containing the self-energies the following approximations 
are performed: 

- The heFd/2 terms in energy dependence of the self-energies are neglected. 

- The expression [G'^'^''{q,k,£) — G'^'^^{q,k,£)] is approximated by 2TTi5{£ — 2Ti cos{qd) — Ek) 
which holds exactly for the free particle Green functions. 

- We use G^ {q, k, £) ^ f{q, k)2TTi6{£ — 2Ti cos{qd) — Ej.). This means that the energetical width 
of the respective states is neglected. As information about quantum mechanical correlations 
is stored in the energy dependence (the Fourier transform of the time difference ti — ^2)1 this 
approximation provides quasiclassical particles with specific momenta q, k. 

Then one finds 



d 



ieF—f{q, k) = S< (k, E, + Ek) - ( ^^"^ (k, E, + Ek) - S^^ (k, Eg + Eu) ] f{q, k 



(B.18) 



^"^ This implies that the density of final states for scattering processes does not vary strongly with energy. This 
variation occurs typically on the energy scale of the miniband width. Therefore this assumption can be justified 
for 2Ti » eFd. 



98 



Now the same approximations are used in the evaluations of the self-energy for impurity scat- 
tering (150): 

S<(k,i?g + Ek) = V(Kk,nk'({^l})V;k',nk({r-;}))imp7^ /" dq' G< {q' ,k' , Eg + Ek) 



v^ d rl"^ (B.19) 

= y ,7r dg 27r(T4k„k/({fj})T4k/ „k({fj}))imp 

^ 27r J_^/a 

xf{q'X)id{Eg+Ek-Eg,-Ek') 

which is just (up to the factor i) the in-scattering term from Fermi's golden rule. In the same way 
the term 5]adv_^rct pj-Qvides the out-scattering rate. It is instructive to note, that the insertion of 
Eqs. (151,152) gives the respective phonon scattering terms including the Pauli blocking factors 
(1 — f{q, k)). Thus, the Boltzmann equation can be derived as a limiting case of the full Kadanoff 
Baym equation if the scattering induced broadening as well as the field dependence is neglected 
in the scattering term. The same results essentially holds if the full space and time dependence 
is maintained, as shown in many textbooks, such as [39]. Some of these approximations can be 
relaxed. In particular, a significant improvement can be made by the generalized Kadanoff Baym 
ansatz [264-267]. 

Finally the current density is evaluated via Eq. (148) which takes the form: 

^^^^.-^El/^I^e{T4£^^d,e^^^G<(.,k,.)} 

r dqRe{T,e"^''if{q,k)} (B.20) 

J-w/d "^ J 

J-TT/d l^ 



2e^p2^ r''^ 



A ^^ h2TT 
k 



2e^p 1 /■'^/'^ ^ -2Tidsin(gd) 

k 



which is just Eq. (37). 

B.3. Wannier- Stark hopping 

In subsection 4.3 it was shown that the Wannier-Stark states become resolved for eFd ^ F. 
Here it will be shown that the equations (64,69) for the self-consistent Wannier-Stark hopping 
model can be derived from the of the quantum transport in this limit. Similar to the derivation 
of miniband transport, the key point is the idea, that the occupation of the Wannier-Stark states 
j can be treated as a semiclassical distribution function fj(k) neglecting quantum mechanical 
correlations. Similar derivations have been presented in [112,252], where essentially the same 
approximations have been applied like in the quantum transport model discussed in this work. 

According to Eq. (31) the Wannier-Stark states are \^j) = ^„ J„_j(/3)|^„) with f3 = 2Ti/eFd. 
The respective Green functions are given by 

G'wSii,J2(k, -E) = 2^ Jmi-ji{l3)Jm2-J2{P)Gmi,ni2{KE) . (B.21) 



mi,m2 
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Within this new basis Eq. (153) becomes 

{E-Eu+ heFd) G^*s,i,,.(k, E) = 6,,,,, + J^ Sws,,,(k, E)G\^%^^^{k, E) (B.22) 



with the definition 

^WSji ,J2 (k, -E) =2_^ Jm-ji (/3) Jm~J2 (/3) ^m {K E) 



(B.23) 



m,J3,J4 k' 

where the Born approximation was apphed under the assumption, that scattering is diagonal in 
the Wannier basis, independent on the well number m, and that no correlations exist between 
different wells m. Furthermore we restrict ourselves to impurity scattering with matrix element 
(^m,k'|-^'^'^'^**l^m,k) = ^'k for simplicity. In the same way Eq. (155) gives 



Ui-h)eFdG<,s^^^^^ik,E)=Y^ 



^wsji,i(k, E)Gyj^jj^{'k, E) + i;^gj-^_^.(k, S)Gwsj,j2(k, E) 



- G'wsji,i(k, -^)^wsi,i2(ki E) - G^g^-^^j(k, £')S^sjj2(k' E) 

(B.24) 

The typical energy scale is given by eFd. If the self-energies (of the order of F) are small in 
comparison to eFd, Eq. (B.22) gives 

1 



^WSil,J2(kl-^) ^ ^jl,J2 



(B.25) 



'^^E-Ek+3ieFd + \Q+ 

describing free-particle Wannier-Stark states. The occupation of these states is governed by a 
semiclassical distribution fj^ (k) yielding 

G<,s^.^^^.^(k,i?) « 2m5,,,,J),{k)5{E - E^ + jieFd) (B.26) 

Let us first consider the case ji = J2, when the left-hand side of Eq. (B.24) vanishes. Applying 
the approximations (B.25, B.26) in the scattering term on the right-hand side and performing 
the integration l/(27r) J dE we obtain 



= i 



^'^%,,j, (k, Ek - jieFd) - ^Tsn,n (^^ Ek - JieFd)\ /,, (k) + i^^sjun (^' ^^ " ^'i^^*^) 

(B.27) 

Inserting Eq. (B.23) one obtains together with the approximations (B.25, B.26): 
= ^ [Jm-jAP)f l^k'kl' [Jm-jM? ^T^HEk - heFd - {Ek' - J2eFd)) [/,,(k) - /,,(k')] 

j2,"^ik' 

=fi^ [^ii,k^i2,k'/ji(k) - /?J2,k'->ji,k/j2(k')] 



k'i2 



(B.28) 
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where Eq. (67) has been inserted. This is just the condition for self-consistency (69) in the 
stationary case for the Wannier-Stark hopping approach. 

The current is determined from Eq. (148) which can be rewritten as 



k '' 

^ E E mM(3)Jm-n if^) I ^Re {neFdG<,s,,,oiK E) } 



Ah 

il,k m 

(To verify this idendity insert all definitions into the second line and use J dEG^^(\i.^E) = 
J dEG^_^Q(k,E) for the homogeneous system) Now Eq. (B.24) can be inserted, where the 
approximations (B.25,B.26) are applied to the scattering terms. This yields: 



Ah 

jr'l ,k m 



/o(k)Im{S^ts,^,o(k,i?fe)} - hm[^<,s,,,,{k,Ek)} 
hm {s<,g^„o(k, Ek - jieFd)} + /,, (k)Im {s^g,„o(k, E^ - jieFd)} 



(B.30) 
where it has been used that S^ is purely imaginary. Now we define the auxiliary function 

/aux(/ - j) = E l^k'kl V<5(i5;fc - jeFd - Ek' + j'eFd) [/,(k) - /,v(k')] (B.31) 

k,k' 

which is an odd function of the difference j' —j for a homogeneous situation, when the occupation 
functions are independent from the index j. Inserting Eq. (B.23) the first two summands of 
Eq. (B.30) yield: 

^ Y, mJ„mJm-n W)Jn-h iP)MP) [Jn-n (/?)]' /aux(j2) 
ji,m,n,J2 

= ^ E '^ [«^n(/3)]^ [Jn~J2W)f /aux(j2) (B.32) 

= %, E h[Jn{P)f[Jn-Km''Un.{h) 
n,h>l 

where in the last line the relation J2n "■ [Jn{P)] \Jn-j2 (/?)] = ^ Y^n [Jn{P)] [Jn-j2 (/5)] and the 
symmetry for j2 = =b/i was applied. Comparing with the relation (67) one obtains: 

^ E [^o,k^?.,k'/o(k) - i?/.,k'^o,kA(k')] (B.33) 

/i>l,k,k' 

which is just half the current density for Wannier-Stark hopping (64). Similarly the last two 
summands of Eq. (B.30) yield: 
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— ^ mJrn{fi)Jm~h iP)Jn-n (.P)Jn(.P) [Jn~J2 iP)? /aux(j2 " j'l) 
ji,m,n,J2 

^ Ah ^ {m'+ji)Jm'+ji{(3)Jm'{P)Jn'(.P)Jn'+jAP)[Jn'-h(.P)fUux{h) (B.34) 

jl,m',n' ,h 

n' ,h 

where in the second hne n' = n — ji, m' = m — ji, and h = J2 — ji have been introduced. In the 
derivation of the last hne, the term with prefactor j'l vanishes by performing the m' sum. The 
final expression is identical to the second line of Eq. (B.32) and thus the full current density 
•AvsH for Wannier-Stark hopping (64) is recovered from Eq. (B.30). 



102 



Appendix C. Quantum transport under irradiation 

In this appendix the quantum transport in superlattices under irradiation with a THz field Fac 
is investigated. Together with a static field F^c one obtains a potential with diagonal elements 

Un{t) = —eFdcdn — ah^cos{Qt)n (CI) 

in the Wannier basis, where a = eFg^cd/h^ is the ratio between the radiation field strength 
and its energy quantum. Similar to section 4 the theory of nonequilibrium Green functions 
is applied. The respective equations are derived in subsection C.l for transport in the lowest 
miniband. In subsection C.2 it is shown that Eqs. (197-199) together with Eqs. (202,203) hold 
for sequential tunneling between equivalent levels a ^ a, where the tunneling matrix element 
Tf does not depend on the field. In contrast, the matrix element H = eF{t)dRi = [eF^cd + 
hQa cos{Qti)]Ri for the transition between different levels [see Eq. (24)] depends on time which 
provides further complications. This will be analyzed in subsection C.3, where it will be shown 
that Eq. (197) for the rectified response still holds in this case. 



C.l. General formulation 

In this subsection essentially the same approximation as in section 4.2 are applied. In particular 
the self-energies Sm (^3,^4) are assumed to be diagonal in the well index. For simplicity we 
set -E" = here. Neglecting scattering and coupling between the wells, the bare Green function 
reads 

g^^^ih ^2) = -iQ(h - ^2) e~'*^'^*=~^-^'l<^'^'^)*^*l~*2) gian[sin(Qti)-sm(nt2)] (C.2) 

for the potential (C.l). Similar to [268] we define the on-site evolution 

S^{t) = e-i"™^''^(™) (C.3) 

of the states and apply the following Fourier expansion for retarded functions 

C;;::„(ti, ^2) = 5t,(ti)i^ y di? 5^ e--^*^ e-'^(*-*^) C-:„^,(i?)5„(t2) 



(C.4) 



In particular this definition gives: 

Now investigate 
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T / dt3 5^*(ti,t3)riGJ!^±i„(t3,t2) 



X 



-t3)/?i. 



^ E + rh^-Eu + meF^^d + iO+ ^^ ^ •^±(^^-^) («)G'r±i,n;r, (i?) • 

1'2 



Similarly 

■p / dt3 / dt4 5'Ji^ (ti,t3)S^ (t3,t4)G5^„(t4,t2) 

With these relations Eq. (108) gives the recursion for the retarded Green function: 
{E + rhn -Ek + meFdcd + iO+)GJ^*„.,(S) = Sm,nSr,o 

+ ^1 2^ |.<^(r2-r)(o)G'm+l,n;r2(-^) + •^{r-r2) (")G'm-l,n;r2 (-^)J 
'"2 

+ 2^ ^m;r-r3 (^ + ^3^^)G'J^*„;r3 (-^) 
»"3 

For the advanced functions we define 

C;^d;(ti, f2) = 5t,(tO^ I di? j; e-™^ e-^(*i"*^)/^ C;^'^;^,(i?)5„(t2) 



(C.7) 



(C.8) 



(C.9) 



so that G^'^.^{E) = |Gjf^.^(£')} holds. The Fourier expansion of the lesser functions is defined 
by 
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C<Jti,t2) = 5t,(ti)i^ /di?j;e--^(*^+*^)/2e-^(*-*^)/^C<„^,(i^)5„(i2). (CIO) 

r 

Then we find the product 

^ / dts / dt4G^*„^(ti,t3)i;<^(t3,t4)G^jj^„(t4,t2) 

(C.ll) 



so that the Keldysh relation becomes 



^1'^:^ (C.12) 

X Smi;r-ri+r3 ( -^ ^ ^^ ) G'miV;r3 (^ ~ (2 "^ ^^) ^) ' 

Eqs. (C. 8, C.12) ahow for a self-consistent solution provided the functionals for the self-energy are 
known. Within the self-consistent Born approximation they are given by the same functionals 
as in section 4.1.2 where the same index r is added in the self-energies as well as the Green- 
functions. This diagonal structure in r is due to the fact that the scattering matrix element is 
either not time dependent (for impurity scattering) or only depends on the time difference ti — ^2 
(for phonon scattering). 

Finally, the current is given by 

In^n+l =^ j; Re{riG<+i,„(t,t,k)} 

■^ ^ (C.13) 



^T.Y.J ^^^ E ^s+.(a)riG<+i,„^,(i?, k) > 



k r 

where S'„_,_;^S'„(t) = ^^/ Jr'{a) e^^ ^* has been used. Note that for homogeneous systems 

Gm,n;riE) = Gm-nfl;riE + ncFdcd) and Iln;r(-E') = ^0;riE + ncFdcd) (C.14) 

holds which can simplify the calculation significantly. 

C.2. Sequential tunneling 

Now we want to derive the expression for sequential tunneling used in section 6.3. Like in 
Appendix B.l the lowest order in the coupling yields a current ~ Tf. For vanishing coupling 

Eqs. (C.8,C.12) give Gm{%2;r{KE) = 5mi,m2^r,oGmi{KE + mieFdcd) and T,m ""^(k^E) = 

T,m-r'^ (k, E + meF(icd)Srfi which are [up to the shift in the energy argument which eliminates the 
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Fdc-dependence of Gm,(k, E)] the same functions as applied in Appendix B.l without irradiation. 
Thus 



G<{k,E) = &^^{k,E)t<{KE)G^:^^KE) = iA{k,E)nF{E - ^„) . 

In lowest order of the coupling one finds from Eq. (C.8): 

GZli,m-A^, E) =&;^\^{k, E + {m± l)eFdcd + rhn)TiJ±r{a)GZ\K E + meFdcd) 
+ 0{Tl) 
GZ]n.A^, E) =5rfiGZ'{k, E + meF^.d) + 0{tI) 
G'^%,.,{k,E)=0{Tl) for \m-n\>2 

and in a similar way: 

G^^±i,m;r(k, E) =G^^ii(k, E + {m ± l)eFdcd)Ti J^,(a)G^^"(k, E + rh^ + meF^cd) 
+ 0{Tl) 



(C.15) 



(C.16) 

(C.17) 
(C.18) 



(C.19) 



With these relations Eq. (C.12) gives 

G'm+i,m;.(k, E) = G^*+i(k, E + {iTi + l)eF^,d + s/2hn)TiJ,{a)GZ\k, E + meFdcd - s/2hn) 
X S< (k, E + meFdcd - s/2hn)G^''{k, E + meF^^d - s/2hn) 
+ G^Vi(k, E+{m + l)eFdcd + s/2hn)t<^^(k, E + {m + l)eFdcd + s/2hn) 
X G^l^ (k, E + {m + l)eFdcd + s/2/if])ri J,G^^(k, E + meF^^d - s/2hn) 



--\TiJs 



&^_^-^{k,E + eFdcd + shn)A^ik,E)nFiE - fi„ 



+ im+i(k, E + eFdcd + sh9)nF{E + eF^cd + shn - i2m+i)G^^{k, E) 
+ 0{Tf) 



(C.20) 



where E = E + meF^cd — s/2Kl was inserted in the last line. Inserting into Eq. (C.13) and 
sorting with respect to cos(ri7t) and sin(r$7t) we obtain Eqs. (197-199) with 



IdcieFd) 



2e 



^'S/ 



d E 

—A'„%,(k,E + eFd)Amik,E) 

[npiE - Hm) - ripiE + eFd - Hm+i] 



(C.21) 



and 



K{eFd) 



4e 



T?T. 



dE 

'2^ 



Re{G::^*+i(k,^ + eFd)}Ara{k,E)np{E - ^i^) 

+ Am+i(k, E + eFd)nF{E + eFd - /i„+i)Re{G'^^"(k, E)} 
This proves the formulas applied in section 6.3 for a ^ a tunneling. 



(C.22) 
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C.3. Tunneling between different levels 



Now we want to consider tunneling between different levels. We will focus on the transitions 
from level a in well m = to the level b in well m = 1 and omit the well indices in the following. 
According to Eq. (24), the matrix element for this transition [eFdcd + hQacos{Qti)]Ri becomes 
time dependent. Neglecting the coupling one obtains again the functions Ga/bO^j^) as in the 
preceding subsection. In the lowest order Eq. (108) yields 

(iS J- -Et-E' + eFa,d\ 0;""^'(k; (,, (2) ^[eFi^d + Kla cos(f!(,)]/;^G;';y"''-(k; h,h) 



+ / ^Er'-(M.,*)G;:r>;M.) 



(C.23) 



In the same way as Eq. (C.6) we obtain 
1 



h 



dt3gr'{ti,t3)[eF^cd+hnacos{nt3)]RTGl''{t3,t2) 



-Siih 



E^-Ek + eFdcd + iO^ 



27r 



-1^1 dt,Slit,)l-l dE,e--^^^('^-':^y^Ua)- 



X e 



-isQt^ 



eFficd + Mia 



giSlta _j_ g— iSlis 



R 



ba 



2tt 



l^S\{tr)S,{t2) I dE^ e-^i*i/^ J] J,(a) — 



E'' -Ek + eFdcd + iO+ 

dE2e''^^^'^-'^y^Gl'\E2)So{t2) 
1 



2-K 
xR\^ 



E^-Ek + eFdcd + iO 



eFdcde'^^^/''-'^>'/''Gl''\Ei - shft) 



+ 



2 
hO,a 



,iiE,/h-{s-l)n)t2/hQVct^-^^ _ ^^ _ ^)^^y 



^5l(ti) / d^^e-i'-™ie-'^(*i-*2)/ft5o(i2) 



1 



271" '''^'\l ^^ "^'^'^^ E + rhn- E^ - Ek + eFdcd + iO+ 

Jr{a)eFdcd + Jr+i{a)hQa/2 + Jr-i{a)hna/2] Gl^\E) 



xR\^ 



(C.24) 



Therefore we find 

-iret /I, 771 \ /^rct 



Gl%,{\i, E) = Gr(k, E + eFAcd + rh^)Jr{a)R\^[eFAcd + r/iJ^JCf (k, E) + 0(Tf) (C.25) 
and the Keldysh relation becomes in analogy to the preceding subsection: 
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Gb,a,sO^^ E) =iJs(a)i?^[eFdcd + sh^] [Gf\K E + eF^.d + shn)Aa{]^, E)nF{E - ^o) 

+ i6(k, E + eFAcd + sh^)nF{E + eF^cd + shft - ^^Gf ^(k, E)] + 0{Tf) 

(C.26) 

Finally, the current density can be obtained by inserting the time dependent matrix element 
into Eq. (C.13) 

la^, =^Y.Y.[ ^^^ (E Js+Aa)[eF^,d + ;iOa(e'^* + e-™)/2]<G<,^,(ii;, k) e-'^4 

k r' -^ I s J 



4e 



E E / |^I^« (E Js+ria)[eF^.d +{s + r)h^]R\^G<^.^,[E, k) e-^4 



(C.27) 

For r = 0, Eq. (197) can be recovered even in the case of a linear field dependence of the matrix 
element. In contrast, for r > 1 Eqs. (198,199) only hold if the field dependence of the matrix 
element is negligible, which is appropriate if eFg,cd <^ eF^cd and hQ <^ eF^cd hold. 
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